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INTRODUCTION 

The  Turbine  Engine  Diagnostic  Development  (TEDD)  Program  began  in 
FY  1972  when  the  Naval  Air  Systems  Command  (NAVAIR,  AIR-330)  authorized 
the  Naval  Air  Propulsion  Test  Center  (NAPTC)  to  develop  and  demonstrate 
turbine  engine  dianostics  by  full-scale  engine  testing  (reference  1) . 

During  FY  1972,  NAPTC  accomplished  Phase  I of  the  program,  which  con- 
sisted of  generating  and  testing  a computer  program  which  would  track 
a TF30-P-408  engine  operating  at  sea  level  conditions  and  output 
diagnostic  messages  in  real-time  to  the  test  cell  operator.  A unique 
vibration  monitoring  system  was  used,  as  well  as  new  type  transducers 
for  oil  quality  and  speed  sensing. 

In  Phase  II,  the  computer  program  was  restructured,  using  new  engine  ^ 

operating  mode  recognition  logic  and  providing  for  thermodynamic  ; 

performance  analysis.  New  hardware  items  used  included  a radiation 
pyrometer.  Internally  mounted  accelerometer  and  real-time  graphic 

display.  The  range  of  engine  testing  was  extended  to  Include  sea  level  !! 

ram  conditions.  The  results  of  Phase  II  were  reported  in  reference  2.  | 

With  the  completion  of  Phase  II  in  FY  1973,  full-scale  engine  ] 

testing  was  discontinued.  However,  a large  quantity  of  data  was  1 

recorded  for  use  in  the  follow-on  software  evaluations.  It  was  planned 
that  follow-on  TEDD  sensor  evaluations  would  be  performed  in  conjunction 
with  other  NAPTC  engine  test  programs  on  a piggy-back  basis.  Reference  3 
authorized  NAPTC  to  continue  the  exploratory  development  effort  in  Task 
Area  WF41.433.301,  Advanced  Auxiliary/Equipment  Systems,  in  accordance 
with  Work  Unit  Plan  No.  NAPTC-624  of  June  1973.  In  TEDD  Phase  III,  this 
diagnostic  development  effort  was  continued  in  a number  of  hardware  and  1 

software  areas  where  further  development  work  was  identified.  These  j 

areas  were  determined  from  Phase  II  test  results,  as  well  as  from  an  j 

overall  knowledge  of  future  diagnostic  requirements  and  the  current  | 

state-of-the-art.  I 

j 

SUMMARY  j 

i 

Since  the  TEDD  Phase  III  effort  addressed  the  development  and 
demonstration  of  diagnostic  elements  (sensors,  software  techniques) 
rather  than  a complete  diagnostic  system,  this  summary  is  an  item  by 
item  description  of  these  efforts. 

A.  Hardware 

I 

1.  Turbine  Inlet  Temperature  Sensor  i 

Direct  measurement  of  turbine  inlet  temperature  in  modern  Navy  j 

gas  turbine  engines  is  not  attempted  due  to  the  lack  of  durability  of  ] 

sensors  located  in  these  regions  of  elevated  temperature.  If  such  a 1 

measurement  was  available,  diagnostic  capability  would  be  improved,  ' 

particulary  in  the  areas  of  gas  path  analysis,  hot  section  life,  and 

overtemperature  severity  determination,  ■ 
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b.  In  TEDD  Phase  III,  NAPTC  contracted  with  Tyco  Laboratories 
to  develop  a turbine  Inlet  gas  temperature  probe  based  on  the  concept 
of  an  incandescent  iiranersed  sapphire  light  pipe.  The  tip  of  the 
sapphire,  being  immersed  in  the  hot  gas,  becomes  incandescent  and  emits 
light  which  is  transmitted  through  the  sapphire  pipe  itself  to  a 
photosensor.  Since  the  emission  is  directly  related  to  the  temperature 
of  the  probe  tip,  it  is  a good  indication  of  that  temperature.  Sapphire 
(AI2  O3)  is  one  of  the  few  materials  that  can  be  used  in  the  high 
temperature,  cxidizirig  environment  at  the  turbine  inlet. 

c.  A prototype  probe  was  fabricated  and  tested  in  a burner 

rig  at  NAPTC.  Although  the  amount  of  testing  was  limited,  and  the 
rig  temperature  was  limited  to  700°C  (1300°F),  the  results  were 

encouraging.  The  basic  concept  that  sufficient  light  will  be  generated 
by  the  probe  tip  to  allow  accurate  measurement  of  the  temperature  was 
demonstrated,  and  some  indication  of  mechanical  integrity  was  obtained. 

Total  response  time  needs  improvement,  although  the  response  of  the 
probe  tip  itself  is  good. 

2.  Ferrograph  Oil  Monitor  - Iron  in  oil  is  one  of  the  prime 
indicators  of  problems  in  bearings  and  gears.  To  date,  no  airborne 
device  can  quantitatively  identify  the  amount  of  iron  in  oil.  The 
Ferrograph  is  a device  which  precipitates  magnetic  particles  from  the 
oil  by  size  through  the  use  of  a variable  magnetic  field.  It  then 
measures  the  optical  density  of  the  larger  and  smaller  particles. 

The  ratio  of  particle  densities  is  then  computed  and  can  be  related 
to  the  health  of  oil  wetted  components.  NAPTC  has  contracted  with 
Trans-Sonics , Inc.  to  take  the  Ferrograph  concept  and  develop  a proto- 
type real-time  engine  mountable  unit  for  test  cell  evaluation.  A 
unit  was  delivered  to  NAPTC  and  evaluated  on  a bearing  rig.  Preliminary 
analysis  of  the  rig  testing  indicated  the  instrument  is  capable  of 
detecting  incipient  bearing  failures  related  to  fatigue  spalling  and 
abnormal  wear.  Further  evaluation  of  the  unit  is  being  made  on  a 
running  engine  (TF34),  and  it  is  planned  to  develop  an  airborne  unit 
for  flight  evaluation. 

3.  Vibrating  Cylinder  Pressure  Transducer 

a.  The  vibrating  cylinder  pressure  transducer  (VCPT)  is  a 
relatively  new  type  of  pneumatic  pressure  transducer  which  utilizes 
natural  vibration  modes  of  a metal  cylinder  to  generate  a frequency 
output  as  a function  of  applied  pneumatic  pressure.  Its  frequency- 
modulated  output  makes  it  immune  to  noise  and  drift  to  which  amplitude- 
based  transducers  are  susceptible 

b.  A VCPT  manufactured  by  Hamilton  Standard  was  evaluated  at 
NAPTC  during  TEDD  Phase  III.  The  transducer  evaluated  is  an  absolute 
pressure  unit  with  a range  of  0-350,000  newtons/m2  (0-50  psia) . 

Unlike  conventional  transducers,  the  output  of  the  VCPT  is  a non- 
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linear  function  of  applied  pneumatic  pressure.  Therefore,  the 
measured  pressure  is  calculated  from  a manufacturer  supplied  polynomial 
which  relates  frequency  output  to  the  applied  pressure.  An  additional 
manufacturer  supplied  temperature  c^upensation  polynomial  is  also 
utilized  in  the  calculation. 

c.  The  evaluation  consisted  oi  laboratory  caiiorat ions  and 
installation  on  a ,15'.''  engine.  Laboratory  calibrations  before  and 
after  the  engine  test  determined  the  ,_ransducer  accuracy  to  'je  better 
than  0.2Z  fuit  scale.  The  repeatab  lit"  M'  the  transducer,  after  more 
than  200  hours  of  on-engine  use,  was  determined  to  be  well  within  the 
uncertainty  limit  jf  the  laborato'-'  , .Aa"d  1+  0.u7%  oi  re.ading). 

4.  Oil  Debris  Monitor  - The  K West  oil  debris  monitor  utilized 
during  tlie  lEOU  Phase  I and  II  tests  was  installed  on  a TFJ0-P-412A 
engine  undergoing  arterburner  lighting  tests.  This  device  consists 
of  a polyester  insulating  screen  interwoven  with  stainless  steel  con- 
ductors, When  conductive  particles  are  trapped  on  the  surface  of  the 
screen,  an  electrical  "short"  occurs,  causing  a decrease  in  total 
resistance  at  the  sensor  terminals.  The  unit  was  installed  in  the 
scavenge  oil  line  for  main  bearings  no.  4,  4 1/2,  5 and  6 on  the 
TF30-P-412A  engine  for  a total  of  65  hours  test  time.  There  was  no 
change  in  the  output  reading  during  the  test,  indicating  the  integrity 
of  the  oil  system  components.  This  agrees  with  the  results  of  SOAP 
samples  for  this  engine.  Thus,  no  unsatisfactory  operation  nas  been 
noted  with  the  K West  unit  during  any  of  the  TEDD  phases,  but  it  has 
not  been  subjected  to  an  engine  problem  situation  where  its  effective- 
ness could  be  evaluated. 

5.  Hot  Section  Analyzer 

a.  A Howell  Instrument  Computer  Recorder  was  insta'led  on  a 
TF30-P-412A  engine  which  was  undergoing  a 150  hour  accelerated  cyclic 
endurance  test  at  NAPTC.  This  unit,  also  referred  to  as  a Time 
Temperature  Recorder  and  Integrator  (TTRl),  has  inputs  of  turbine 
inlet  temperature  (indirectly  measured)  and  high  rotor  speed,  and 
provides  a cumulative  Hot  Section  Factor  Count  based  on  the  Integrated 
time-temperature  exposure  and  the  hot  section  metallurgical  character- 
istics of  the  engine.  It  also  displays  Hot  Section  Factor  Time 
(1015°C,  1900“F  and  above),  time  above  il7/°C  (2150°F),  total  engine 
time,  and  maximum  temperature,  and  has  three  over temperature  flags 
and  one  overspeed  flag.  The  unit  performed  .sacj.sf actorlly  throughout 
the  duration  of  the  test. 

b.  This  type  of  diagnostic  an.ilysls  has  the  potential  of 
impacting  maintenance  philosophy,  that  is  to  perform  maintenance  "on- 
conditlon"  rather  than  at  the  current  fixed  time  intervals.  In  ordei 
to  accomplish  this,  a correlation  of  Hot  Section  Factor  Counts  and 
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hot  section  deterioration  must  be  determined  from  a large  sampling  of 
engines  in  fleet  service.  This  unit,  which  is  limited  because  of 
the  number  of  parameters  measured,  may  not  obtain  sufficient  infor- 
mation to  generate  a reliable  correlation  for  hot  section  condition. 
However,  this  approach,  expanded  to  include  low  cycle  fatigue,  thermal 
shock  and  other  hot  section  distress  mechanisms,  offers  the  potential 
for  significant  improvement  in  engine  maintenance  philosophy. 

6 . Vibration  Analyzer 

a.  he  vibration  analyzer  was  built  by  General  Electric 
Company,  Binghamton,  New  York,  for  the  TF34--GE-2  and  subsequently 
adapted  for  computer  control.  There  was  no  problem  in  using  it  on 
the  TFlU-P-AOB  s"'nce  the  system  is  easily  converted  to  any  engine  by 
changing  ratios  and  limits, 

b.  The  analyzer  contains  signal  processing  electronics  to 
evaluate  four  classes  of  malfunctions  in  an  engine. 

(1)  Bearings 

(2)  Mass  unbalance 

(3)  Gears  - local  defect 

(4)  Gears  - gross  defect 

c.  Digital  Comb  Filter  - In  the  latter  three  processes,  signal 
extraction  is  done  by  the  digital  comb  filter.  The  digital  comb  filter 
is  a time-averaging  device  of  256  discrete  points.  The  time-averaging 
feature  tends  to  cancel  noise  on  the  signal.  The  comb  filter  is 
synchronized  with  a rotating  member  so  that  each  of  256  points 
examines  the  same  point  on  a rotating  member.  The  filter  has  responses 
at  its  tuning  frequency  and  integer  multiples  (harmonics)  of  it.  Any 
signals  or  noises  which  are  not  exactly  Integer  multiples  of  the 
tuning  frequency  will  be  rejected.  The  significance  of  this  type  of 
filtering  can  be  seen  when  a vibration  signal  from  a gearbox  is  applied 
to  the  filter.  This  signal  will  consist  of  the  sum  of  all  gearbox 
shaft  and  meshing  vibrations.  By  tuning  the  filter  to  a gearshaft 
frequency,  only  the  shaft  and  its  gear  meshing  vibration  are  passed 
through  the  filter,  since  its  frequency  (number  of  teeth  x shaft 
frequency)  is  a shaft  harmonic. 

d.  Bearing  Malfunction  Processor  - The  bearing  malfunction 
discriminant  processor  evaluates  the  Impact  Index  of  the  acceleration 
signal  from  a bearing  housing.  The  Impact  Index  is  a normalized 
dimensionless  quantity  whose  value  is  indicative  of  the  incipient 
bearing  malfunction.  The  Impact  Index  value  is  one-half  the  ratio  of 
the  peak  signal  level  to  its  average  level.  An  Impact  Index  for  a 
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normal  bearing  will  range  from  two  to  three  and  during  spal]  initiation, 
from  three  to  four.  As  the  spall  increases  in  size  (but  is  still 
relatively  small),  the  Impact  Index  may  increase  to  eight  t^  ten. 

Beyond  this,  with  Increasing  spall  size,  the  Impact  Index  will  decrease 
due  CO  t!.e  increase  in  average  acceleration.  The  Impact  Index 
liscrimir.  ri.t  for  normal  bearings  is  not  a sanction  of  engine  spared; 
j-.e.,  it  Will  remain  between  three  and  four  over  the  entire  engine 
operating  range.  Foi-  a bearing  with  malfunction,  the  Impact  index  may 
increase  by  up  to  30  percent  with  increasing  engine  speed.  Full  scale 
on  the  indicator  is  an  Impact  Index  of  10. 

e.  Ma.- s b'ob.ilance  Processor  - " ic  mass  ur:balance  discriminant 
processor  rs  a narrow  band  tracking  filter  which  selectively  i liters 
vibration  energy  associated  with  the  mass  unbalance  of  a rotor  system. 

Tht  output  Ls  displayed  as  a disnS.acement  on  a meter  calibrated  from 

0 to  100  percent  full  scale,  where  lalx  scale  Is  equal  to  2.54  x 10~^m 
(10  mils)  double  amplitude.  The  signal  is  first  filtered  to  accept 
frequencies  in  the  10  Hz  to  AOO  riz  band.  It  is  then  filtered  by  the 
digital  f omb  filter,  which  is  tuned  to  eiti.er  once  per  fan  or  core 
frequency.  The  output  of  the  comb  filter  now  represents  the  acceleration 
associated  with  the  rotor  mass  unbalance.  A double  integration  then 
yields  the  displacement  associated  with  the  rotor  mass  unbalance.  To 
obtain  a signal  suitable  for  display,  the  displacement  signal  is 
average  detected  and  displayed  on  a meter  calibrated  to  read  peak-to- 
peak  values. 

f.  Local  Gear  Defect  Processor  - The  local  gear  defect 
discriminant  processor  evaluates  the  Impact  Index  associated  with  a 
given  gear  mesh  vibration  signature.  Local  defects  on  a gear  consist 
of  spalled,  deformed,  or  cracked  teeth. 

(1)  A normal  gear  in  mesh  will  generate  a sinusoidal 
vibration  at  the  gear  meshing  frequency.  '/Then  a local  defect  is  present 
on  a tooth  of  the  gear,  a transient  vibration  will  be  generated  each 
time  that  tooth  meshes.  The  level  of  the  transient  will  bt  considerably 
higher  than  the  normal  meshes.  This  type  of  signal  can  easily  be 
discriminated  by  an  Impact  Index  measurement.  Since,  in  practice,  a 
single  gear  mesh  vibration  Is  mixed  with  other  mesh  vibrations  and 
noise,  this  mesh  must  be  extracted  from  the  total  signal  before  its 
Impact  Index  can  be  evaluated. 

(2)  The  gear  discriminant  is  implemented  by  tuning  the 
digital  comb  filter  to  the  shaft  frequency  of  the  gear  being  analyzed. 

At  this  frequency,  only  harmonics  will  be  passed  through  the  comb 
filter.  The  harmonic  of  Interest  will  be  the  nth.  harmonic  corres- 
ponding to  the  gear  under  analysis  which  has  "n"  number  of  teeth. 

After  the  gear  mesh  has  been  isolated,  the  signal  is  evaluated  for 
its  Impact  Index  and  displayed. 
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g.  Gross  Gear  Defect  Processoj  - The  gross  gear  defect  dis- 
criminant processor  evaluates  the  modulation  of  gear  mesh  vibration 
due  to  its  shaft  frequency.  A gross  gear  defect  includes  misaligned, 
loose,  or  eccentric  gears  on  a gearshaft. 

(1)  A normal  gear  in  mesh  will  generate  a sinusoidal 
vibration  due  to  its  shaft  frequency.  When  a gear  has  a gross  defect, 
the  gear  meshing  vibration  will  be  modulated  by  the  gearshaft  frequency. 
Since  this  signal  Is  normally  mixed  with  other  gear  mesh  vibrations 

and  niise,  the  characteristic  gear  mesh  vibration  signature  must  be 
extracc?d  from  the  total  vibi-ation  signal. 

(2)  The  discriminant  measurement  is  implemented  by 
tuning  the  digital  corah  filter  to  the  gearshaft  frequency  of  Interest. 
This  will  allow  Integer  multiples  of  the  shaft  frequency  to  pass 
through  the  filter  while  rejecting  all  other  signals.  The  modulation 
will  be  passed  through  the  filter.  The  measurement  which  is  then 
made  on  the  signal  is  the  ratio  of  the  modulation  amplitude  to  the 
carrier  fgear  mesh)  amplitude  yielding  the  modulation  index,  which 
ranges  from  0 to  100.  This  Is  a dimensionless  parameter.  To  normalize 
the  ratio,  an  automatic  gain  control  (AGC)  is  used.  The  AGC  is  a 
variable  gain  amplifier  that  is  electronically  controlled  to  boost  low 
level  signals  up  to  levels  that  may  be  detected. 

(3)  The  AC  portion,  which  contains  the  modulation,  is  low 
pass  filtered  to  allow  only  the  modulating  frequency  to  pass.  The  low 
pass  filter  is  programmable  so  that  the  filter  cutoff  may  be  set  for 
various  gears  and  various  engine  power  settings.  The  output  of  the  low 
pass  filter  is  then  peak  detected  and  displayed  on  the  output  meter. 

(4)  The  section  of  the  analyzer  of  interest  for  this 
phase  of  the  TEDD  program  was  the  gross  gear  defect  processor.  The 
analyzer  gave  unsatisfactory  readings  for  the  TEDD  engine,  and  experi- 
ments were  conducted  to  determine  why.  Part  of  the  problem  was  due  to 
the  method  of  speed  measurement. 

(5)  To  provide  a speed  signal  for  the  gear  box  experiments, 
the  signal  is  taken  from  a tach  generator  mounted  on  the  gear  box.  The 
signal  has  modulation  on  It  due  to  gear  back  lash  and  play  in  the  gear 
mesh,  which  looks  like  speed  changes  to  the  analyzer.  A percent 
modulation  can  be  measured  by  comparing  the  period  changes  that  occur 

to  a steady  signal  of  the  same  frequency.  In  addition  to  the  problems 
from  the  modulated  tachometer  signal,  additional  filtering  was 
required  to  smooth  the  digital  comb  filter  output. 

7.  Speed  Sensor  - In  conjunction  with  the  vibration  analyzer  effort 
in  Phase  III,  a study  of  various  methods  of  speed  sensing  was  conducted. 
The  gear  box  mounted  tach  used  on  the  TF30  proved  to  be  insufficient  for 
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for  vibration  auaiysis  due  to  gear  backlasn  and  play,  creauang  signal 
modulation.  Tlie  primary  emphafe.s  of  thxs  study  was  the  determination 
of  signa.1  modulation  as  a function  of  pickup  location.  Ihe  locations 
in”  .-.stlgatci.  were  the  gear  bot,  tf'.e  nose  oor.  j ana  a direct  shat  t 
picsip.  A blade  passage  technique  was  alt..  lUvestigatid,  and  this 
techniqLe  otters  acceptable  moduiatio..  i acteristics  tor  use  with 
the  vibratior.  .analyzer. 

Kadi  ’tion  Pvrouiecer  - fhe  Solar  radi.ation  pyrometers  used 
during  TEI  i fl.ase  If  engine  testing  exhibits  stgn.-.l  characteristics 
which  renuired  further  invet Ligation.  This  Pnase  111  sti.ay  addressed 
particuiarl'  the  ucise  in  the  nnajo-  cutprt  ir^-ult.  Dincussiori  with 
other  users  of  rauiatlon  pv,.ometers  revt-ai.-d  tbf  relat  iorsh’ p between 
signal  noi.sc  and  nrobe  location.  Locat'.ng  the  nrobe  such  that  the 
flame,  c.-nnot  be  view^'d  (i.e.,  bebina  the  o.Lade  cow  locking  forward  or 
looking  radially)  resulted  in  reduced  noise  levels.  Another  method 
of  reducing  the  noise  which  was  investigated  is  the  use  oi  the 
vibration  anaj  ,-zer  fitter  circuitry  to  remove  the  noise.  The  filtering 
technique  also  shows  promise  lor  impro 'ing  the  pyrometer  output. 

9.  Internal  Accelerometer  - One  of  the  reasons  vibration  monitoring 
techniques  are  quite  complex  is  the  communication  patfi  between  the 
vibration  sources  and  the  vibration  pickup.  An  accelerometer  which  is 
not  closely  coupled  to  the  bearing  (or  other  vibration  source)  which 
is  to  be  diagnosed  will  pickup  vibration  from  other  sources.  This 
requires  the  analyzer  to  "sort  out"  the  vibrations,  making  the 
technique  complex.  Therefore,  reduction  in  analysis  complexity  v.an 
be  realized  if  accelerometers  are  directly  mounted  to  the  bearing  to 
be  monitored.  There  are  some  common  problems  associated  with  this 
task  that  have  been  under  design  consideration  for  several  years. 

These  problems  are  the  environment.il  restrictions  on  the  acceJ erometer , 
signal  lead  routes,  and  bearing  housing  accessibility,  which  places 
restrictions  on  the  design  of  the  engine.  The  problem  of  limited  life 
for  the  accelerometer  makes  it  necessary  that  it  be  accessible  without 
engine  disassembly.  One  approach  is  to  place  the  accelerometer  at 
the  bearing  through  a "dipstick".  This  requires  a suitable  engine 
design,  but  places  the  accelerometer  ana  lead  in  a position  of  access- 
ibility. This  approach  is  well  suited  to  rigidly  mounted  bearings; 
soft-mounted  bearings  (such  as  Intershatt  bearings  and  other  bearings 
not  rigidly  mounted  to  the  engine  frame)  are  not  conducive  to  close 
coupled  accessible  accelerometer  installation. 

B.  Software 


1.  Data  Smoothing  Techniques  - During  the  TEDD  Phase  II  test 
program,  a limited  moving  average  smoothing  technique  was  used  to  improve 
the  stability  of  the  mode  detection  routine  and  reduce  the  scatter  in 
the  performance  analysis  results.  The  moving  average  technique  means 
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the  parameter  used  in  the  diagnostic  logic  is  the  average  of  the  last  ) 

n values  of  that  parameter  (n  in  TEDD  Phase  II  was  adjustable  between 

1 and  4).  This  technique  (for  n > 1)  will  reduce  the  noise  level  of  ] 

the  smoothed  parameter  and  will  reduce  the  rate  of  response.  It  was  ’ 

recognized  that  further  investigation  was  required  in  the  area  of  i 

data  smoothing,  and  especially  the  consideration  of  exponential  ? 

smoothing.  Exponential  smoothing  is  a .net  ■ >d  of  estimating  each  new 

parameter  value  in  a time  series  ot  data  by  letting  the  new  smoothed  | 

value  be  equal  to  a fraction  (1  - a)  oi  the  previous  smoothed  value,  j 

plus  a fraction  (a)  of  ihe  new  parameter  value.  u is  defined  as  the  j 

weight  given  to  the  new’  parameter  value.  A study  of  data  smoothing  I 

was  made  during  TEDD  Phase  III,  using  data  recorded  during  Phase  II  ^ 

and  also  from  another  NAPTC  engine  test  program.  As  a result  of  this  j 

study,  it  was  shown  that  the  exponential  smoothing  technique  is  better  i 

Suited  to  TEDD  applications  than  the  moving  average  technique. 

2.  Multiple  Fault  Gas  Path  Analysis  - The  NAPTC  steady-state 
performance  analysis  used  in  TEDD  Phase  II  was  designed  to  detect  single 
gas  path  performance  faults.  However,  both  single  and  multiple  fault 
conditions  were  tested  and  recorded.  This  actual  test  data,  acquired 
at  knowm  fault  level  conditions,  provided  a means  of  evaluating  the 
Hamilton  Standard  "Parameter  Interrelationships"  approach  to  multiple 
fault  steady-state  performance  analysis.  During  TEDD  Phase  III,  NAPTC 
contracted  (Contract  No,  N00140-74-C-0582)  with  the  Hamilton  Standard 
Division  of  United  Technologies  Corporation  to  apply  and  demonstrate 
their  methods,  using  the  TF30-P-408  engine  data  recorded  during  Phase  II 
testing.  The  results  of  this  program  verify  the  capability  of  this 
approach  to  identify  and  quantify  single  and  multiple  gas  path  per- 
formance faults.  A further  study  is  being  conducted  by  NAPTC  in-house 
to  determine  the  sensitivity  of  this  approach  to  "time  series"  data, 
i.e.,  continuous  real-time  steady-state  performance  analysis. 

3.  Diagnostics  on  Variable  Geometry  Engines  - The  need  for  the 
capability  to  simulate  the  performance  of  a gas  turbine  engine  with 
variable  geometry  arose  from  the  consideration  of  gas  path  analysis 

as  one  element  of  a diagnostic  package  for  a 1980-85  aircraft.  During 
the  period  from  January  1974  through  March  1975,  a variable  geometry 
engine  cycle  deck  was  developed  at  NAPTC.  The  first  application  of  the 
data  from  this  cycle  deck  was  to  provide  variable  geometry  data  to 
the  Hamilton  Standard  Division  of  United  Technology  Corporation  for 

study  using  the  parameter  interrelationships  method  of  gas  path  ‘ 

analysis.  The  purpose  of  this  study  was  to  identify  unique  problem  1 

areas  associated  with  diagnosing  variable  geometry  engines.  2 

Additionally,  this  study  identified  critical  instrumentation  accuracy  fl 

requirements.  This  work  was  conducted  under  Contract  No.  N00140-  1 

75-C-0049,  from  1 April  to  1 October  1975.  3 


8 


NAPTC-PE-88 


IHDD  Seniattivlty  Study  - The  ability  of  diagnostic  logic  to 
perform  Its  > iteni-ad  task  is  limited  by  the  accuracy  and/or  repeat- 
abili.y  of  the  incoming  data.  The  characteristics  of  candidate 
sensort  must  be  viewed  with  respect  to  their  application,  that  is, 
can  the  sene-  r nerform  adequately  for  the  logic  using  that  parameter. 
Sensitivity  studlrS  establish  the  relationship  between  sensor  per- 
formance and.  e diagnostic  logic  requirements  for  that  sensor. 

>)'ir1r>c<  " 'e  Tifi,  1 sensitivity  study  was  performed  utilizing  the 
gso  path  ; :eci-i.ique  developed  during  TEDD  Phase  II.  The 

results  ' cn'.s  m-noy  ■’ar  be  used  to  determine  the  sensor  character- 
istic-J  requiiod  i-,r  a .tevel  of  confidence  in  the  gas  path 

analysis  results.  Another  aspect  of  the.  study  considered  the  sensor 
chaia.;r.fc.  j.st  ' ts  i - iiulrei  for  hot  section  life  rnniput.ation, 

CnNCLuSlutiS 

1.  Lntttai  tig  testing  of  a prototype  sapphire  light  pipe  turbine  inlet 

temperatuci  sensor  has  demonstrated  the  feasibilitv  ,.if  the  basic  concept 
CO  705'’C.  and  has  shown  encouraging  results  with  respect  to 

sensitivity,  strength  anJ  potential  time  response. 

2.  Analysis  of  data  acquired  on  a bearing  rig  indicates  that  the  engine 
moui^table  real  time  Ferrograph  is  capable  of  detecting  abnormal  wear  and 
bearing  failure  in  each  case  tested. 

3.  ITie  vibrating  cylinder  pressure  transducer  operated  satisfactorily 
and  is  superior  to  the  conventional  diaphram-type  transducer  in  terms  of 
accuracy,  stability,  repeatability  and  ease  of  operation.  Its  major 
drawbacks  are  its  relatively  slow  response  and  its  sensitivity  to  varia- 
tions in  gas  cot; position. 

4.  TTie  K West,  oli  debris  monitor  has  been  installed  on  three  TF30  engine 
tests  and  has  performed  satisfactorily.  However,  it  was  not  subjected  to 
an  engine  oil  system  fault  situation  where  Its  effectiveness  could  he 
evaluated . 

5.  The  Howell  Hot  Section  Analyzer  (Time  Temperature  Recorder  and  Integ- 
rator) operated  satisfactorily  on  a TF30-P-412A  engine,  but  more  extensive 
operational  experience  would  be  necessary  as  an  end  item  to  verify  the 
potential  maintenance  benefits  for  any  particular  engine.  It  should  be 
noted  that  the  TTR&I  does  not  record  LCF  cycles  or  LCF  rupture  inter- 
actions nor  has  a correlation  between  TTR&I  output  and  hot  section  distress 
been  developed. 

6.  Installation  of  a filter  before  the  automatic  gain  control  success- 
fully eliminated  noise  caused  by  the  digital  comb  filter  of  the  General 
Electric  vibration  analyzer;  however,  at  high  engine  speeds,  gear  gross 
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defect  experiments  were  still  over  limits  due  to  frequency  modulation  on 
the  tachometer  signal,  indicating  tliat  blade  passage  speed  data  would  be 
satisfactory  for  use  with  the  GE  vibration  analyzer. 

7.  Noise  in  the  radiation  pyrometer  output  can  be  reduced  by  placing 
the  probe  where  it  cannot  view  the  flame. 

8.  Circuitry  In  the  GE  vibration  analyzer  could  be  used  to  process  and 
analyze  pyrometer  signals. 

9.  The  problems  of  mounting  internal  accelerometers  can  be  overcome  if 
they  are  addressed  during  initial  engine  design.  Where  the  bearing  is 
soft  mounted,  there  is  a communication  problem  between  bearing  and 
accelerometer  if  the  accelerometer  is  not  coupled  directly  to  the  bearing. 

10.  The  exponential  data  smoothing  technique  is  better  suited  for 
diagnostic  system  application  than  the  moving  average  technique.  To  apply 
exponential  smoothing  in  a diagnostic  system,  it  is  necessary  to  have  pre- 
defined separate  a's  for  steady  state  and  transient  modes,  depending  on 
the  response  and  accuracy  requirements  for  each  parameter. 

11.  The  capability  of  the  Hamilton  Standard  multiple  fault  diagnostic 
approach  to  identify  and  quantify  single  and  multiple  gas  path  performance 
faults  has  been  verified. 

12.  The  variable  geometry  cycle  deck  developed  for  diagnostics  studies 
on  advanced  engines  is  a very  flexible  deck  with  many  potential  uses  in 
addition  to  diagnostics,  such  as  the  study  of  the  effects  of  advanced 
technology  on  cycle  performance,  and  the  evaluation  of  performance  quotes 
in  an  engine  development  proposal. 

13.  When  a non-dimensional  technique,  is  used  for  diagnosis,  such  as  in 
the  TEDD  Steady  State  Performance  Analysis,  the  repeatability  of  the  data 
is  more  important  than  the  absolute  accuracy  of  the  data. 

14.  The  TEDD  turbine  life  computations  are  extremely  sensitive  to  turbine 
blade  temperature  error.  Calculated  erosion  life  will  be  in  error  by  10% 
for  a 5.6°C  (10°F)  error  in  gas  temperature  while  creep  life  shows  35% 
error  for  a 5,6°C  (10°F)  error  in  gas  temperature. 

15.  High  and  low  rotor  speeds  (N2  and  Ni)  are  the  most  sensitive  para- 
meter in  the  TEDD  Steady  State  Performance  analysis  at  the  sea  level  static 
condition  studied. 

RECOMMENDATIONS 

1.  Further  development,  test  and  evaluation  of  a sapphire  light  pipe 
turbine  inlet  temperature  sensor  should  be  pursued  to  determine  the  sen- 
sor's accuracy  and  response  characteristics,  and  rig  testing  should  be 
extended  up  to  the  limits  of  the  probe. 
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2,  Tl'ie  real  Lime  Ferrograph  should  be  lurther  evaluated  on  running 
engines  In  a test  cell,  and  work  should  continue  to  develop  an  airborne 
uni  t . 

j.  A n-.i.iioer  of  possible  failure  mechanisms  shoi.lu  be  investigated  for 
engine  hot  section  life  diagnostics,  incl'iding  creep,  low  cyi le  fatigue, 
thermal  shock  and  erosion,  because  current  analytical  approacnes  do  not 
provide  sufficient  accuracy  for  successful  on-condltion  maintanance. 

4.  Additional  vibration  analyzer  analvsis  should  be  dene  wiih  data  that 
has  a ta«'homc-ter  signal  with  less  than  0.8%  modulation,  to  resolve  th< 
question  of  i req  i"nce  modulation  versus  lack  of  cc:imui  icatici.,  to  provide 
valid  gear-gross  experiment  results. 

5.  For  future  vibration  analyzer  use,  a blade  passage  speed  pickup 
should  he  used  for  the  high  speed  rotor  and  gear  box  experiments. 

6.  Pyrometer  development  should  be  continued  for  use  in  future  engine 
diagnostic  systems,  witli  careful  consideration  given  to  proper  placement 
of  the  probe  in  new  engine  designs  so  as  to  minimize  inaccuracies  due  to 
reflected  radiation. 

/.  If  Internally  mounted  accelerometers  are  used  where  the  engine  bearing 
is  soft  mounted,  they  should  be  mounted  directly  on  the  bearing. 

8.  Exponential  data  smoothing  should  be  considered  for  incorporation  in 
all  future  iNAP'iC  diagnostic  and  general  transient  data  reduction  procedures. 

9.  Exponential  data  smoothing  should  be  considered  for  use  in  future 
diagnostic  systems  because  of  its  ability  to  reduce  signal  noise  while 
maintaining  adequate  response  characteristics. 

10.  A study  should  be  made  to  determine  the  sensitivity  of  the  Hamilton 
ittandard  multiple  fault  diagnostic  approach  to  real  time  data. 

/ 

< 11.  Investigation  should  be  made  of  the  Hamilton  Standard  multiple  fault 
diagnostic  approach  as  applied  to  an  advanced  variable  geometry  gas 
turbine  engine  to  determine  the  impact  of  variable  geometry  on  the 
analysis  complexity,  storage  requirements,  and  Instrumentation  accuracy/ 
repeatability  characteristics  requirements. 

DISCUSSION 

A.  Hardware 

1 . Turbine  Inlet  Temperature  Sensor 

The  ability  to  directly  measure  turbine  inlet  temperature  would 
be  a very  valuable  diagnostic  tool  as  well  as  being  useful  as  a control 
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parameter.  In  TEDD  Phase  III  NAPTC  contracted  with  Tyco  Laboratories  to 
develop  a turbine  inlet  gas  temperature  probe  based  on  the  concept  of 
an  incai;  Ic  u-eiit  Immersed  sapphire  light  pip^.  The  tip  of  the  sapphire, 
being  immersed  in  the  hot  gas,  becomes  incandescent  and  emits  light  which 
is  transmit tt  i tnrotigh  the  sapphire  pipe  itself  to  a photosensor.  Since 
tne  emission  is  directly  related  to  the  temperature  of  the  probe  tip,  it 
is  a good  indication  f that  temperature.  Sapphire  (AI2  O3)  is  one  of 
the  few  iiuteri.iJs  tiiat  can  be  used  in  the  high  temperature,  oxidizing 
• :ivironi.;e  It  at  t no  t ai  i ne  inlet.  i'he  high  strength  and  hardness  of  the 
s ippiiire  should  i>ermt  t the  probe  to  operate  at  temperatures  up  to  1371°C 
250()“K)  . 

A pro.  t probe  iliouj"  1)  was  t.chric.nted  and  tested  in  a burner 

rig  at  NAi  l'(  . .h..i;  :‘i  the  amount  of  testing  was  limited,  and  the  rig 

t.mnerato  vi,  limited  to  (1300 the  results  were  encouraging. 

. ne  b.isle  .on  t-pt  tliit  ‘-uf  f ieient  light  will  he  generated  by  the  probe 
lip  to  allow  ai  i.-re  me.is  irement  of  the  temperature  was  demonstrated,  and 
some  Ind  ic.ti  ii;u  </f  mech  ■‘r.ical  integrity  was  obtained.  Total  response  time 
needs  imorriV'.  rn  ,it , altiioupb  the  response  of  the  probe  tip  itself  is  good. 
This  development  el f<irt  is  further  described  in  reference  A. 

2 . Ee rrograt  h Oil  Monitor 

('ll  .inalysis  has  long  been  recognized  as  an  effective  method  of 
determining  the  healtli  ot  engine  oil  wetted  parts  as  shown  by  the  Spec- 
trometric  Oil  Analysis  Program  (SOAP).  This  approach  requires  periodic 
sampling  oi  the  oil,  the  delivery  of  the  sample  to  an  analyzer,  and  the 
subsequent  leporting  of  the  analysis  results  to  the  maintenance  personnel 
responsible,  for  that  engine.  This  approach  has  at  least  two  drawbacks: 

(1)  samples  are  taken  as  a function  of  engine  hours  rather  than  an 
indicated  need  for  analysis  and,  (2)  the  time  elapsed  between  sampling 
and  reporting  of  the  analysis  results  may  exceed  the  time  available  to 
perform  effective  preventive  maintenance.  These  drawbacks  are  inherent  in 
any  non-real  time,  non-aJrborne  analysis  approach.  Chip  detectors  are  a 
real  time,  airborne  approach  currently  being  used  which  are  limited  by 
the  large  particle  size  required  to  indicate  an  impending  (?  already 
present)  failure.  Another  real  time,  airborne  approach  has  been  formulated 
by  Trans-S'onics , Inc.  utilizing  their  Ferrograph  concept.  The  Ferrograph 
(Figure  2)  is  a device  which  precipitates  magnetic  particles  from  an  oil 
sample  and  de()osits  them  by  size  through  the  use  of  a variable  density 
magnetic  field.  It  has  been  determined  by  Trans-Sonics  that  the  progres- 
sion of  many  gear  and  bearing  failure  mechanisms  is  directly  related  to 
the  rate  of  Increase  of  the  presence  of  large^  particles  in  the  oil 
(reference  5).  The  Ferrograph  consequentlv  optically  measures  the  density 
of  the  deposited  large  and  small  particles  and  is  thereby  capable  of  in- 
dicating the  presence  of  a failure  mechanism,  that  is  an  incipient  failure. 


^I.arge  In  this  Inst.mce  means  between  2 and  10  microns  in  major  dimension. 
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NAPTC  contracted  with  Trans-Sonics,  Inc.  under  Contract  NOOlAO-73-1065 
to  take  the  Ferrcgraph  concept  and  develop  a prototype  real  time  engine 
mountable  systems  for  test  cell  evaluation.  Reference  5 is  the  Final 
Report  concerning  this  contract. 

For  convenience  of  operation  in  a test  cell,  the  Peal  Time  (RT) 
Ferrograph  Svstem  consists  of  three  units;  (1)  RT  Fe'-”ograph  ro  be 
mounted  on  an  engine  above  the  oil  tank;  (2)  a pump  which  is  normally 
mounted  below  the  ank;  and,  (3)  a Coetro' /] ndic.  tor  designed  to  be  mount  i 
on  the  wall  of  the  control  room  so  that  the  indications  may  bn  viewed  and 
recorded  during  engine  operation.  The  RT  Fcrrc.ijraph  op  -rates  by  passing  a 
fixed  volume  of  ha  engine's  lubricating  oil  through  a precipitator  tube 
located  in  the  field  of  an  electromagnet.  As  tf o oil  flews  slowly  through 
the  tube  in  *-he  electromagnet's  field,  the  wear  -'.art-cies  are  deposited 
on  the  inside  surface  of  the  tube.  The  large  particles  precipitate  first 
and  the  particle  size  grades  continuously  as  a functioc,  o^  distance  along 
the  tube. 

In  order  to  indicate  the  density  of  particles  ..f  a selected  size,  a 
light  beam  is  passed  across  the  tube  and  the  attenuation  cf  the  light  is 
measured.  The  light  is  conducted  from  a common  lamp  assemblv  to  three- 
positions  on  the  tube  by  means  cf  a trifurcated  fiber  optic  bundle.  ihe 
first  position  is  a reference  and  is  external  to  the  magnetic  field. 

Light  crossing  the  tube  is  conducted  by  a fiber  optic  bundle  to  a reference 
photo  resistor.  If  the  color  of  the  oil  changes  from  light  to  dark,  the 
lamp  Intensity  is  automatically  increased  so  that  the  intensity  of  light 
reaching  the  photo  resistor  remains  constant,  independent  of  color.  Light 
from  the  same  lamp  is  also  conducted  to  two  other  positions  on  the 
precipitation  tube  by  branches  of  the  trifurcated  fiber  optic  bundle. 

The  second  position,  close  to  the  point  where  the  oil  enters  the  mag- 
netic field,  measures  the  deposit  of  the  larger  particles.  These  particli 
typically  range  from  2 up  to  10  microns.  The  light  which  passes  up 
through  the  tube  is  blocked  by  the  presence  of  the  opaque  metal  wear 
particles.  The  fractional  area  of  the  light  beam  which  is  blocked  by 
the  particles  determines  the  precentage  change  of  light  intensity  which 
is  measured  by  one  of  the  sensing  photo  resistors  in  the  photo  resistor 
bank. 


The  sensing  photo  resistors  are  connected  to  Integrated  operational 
amplifiers  located  in  the  wall-mounted  Control /Indicator . The  circuits 
are  designed  to  generate  a voltage  linearly  proportional  to  the  percent 
area  of  the  light  beam  covered  by  the  particles.  For  medium  to  low 
densities,  the  particles  are  deposited  in  a single  layer  and,  therefore, 
the  output  voltage  Is  linearly  porportional  to  the  area  of  the  particles 
in  the  field  of  view.  A third  light  beam  is  located  where  smaller 
particles,  on  the  order  of  0.5  microns  are  deposited. 
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At  the  start  of  an  operating  cycle,  the  pump  is  energized  to  circulate 
oil  from  the  engine  through  the  oil  passages  of  the  Ferrograph,  including 
a small  reservoir  of  approximately  lOcc  volume.  During  this  time,  an 
automatic  zero  circuit  adjusts  each  of  the  readout  circuits  to  read  zero 
as  indicated  on  the  digital  displays.  At  the  end  of  the  flushing  cycle, 

B^hich  lasts  approximately  one  minute  as  determined  by  an  internal  timing 
circuit,  the  pump  turns  ol^,  the  vent  line  drains,  and  the  oil  in  the 
lOcc  reservoir  begins  to  run  down  through  the  precipitation  tube.  The 
speed  with  wnich  the  oil  flows  through  the  tube  is  determined  by  a 
restrictor  and  by  the  viscosity  of  the  oil.  If  the  oil  viscosity  is  high, 
the  particles  migrate  through  the  oil  slowly  under  the  influence  of  the 
magnetic  field,  but  the  oil  flows  through  the  restrictor  slowlv  in  exactly 
the  same  ratio  so  that  the  trajectory  of  the  particles  is  nearly  constant 
and  independent  of  temperature.  As  the  oil  heats  up,  the  viscosity 
decreases  so  that  the  oil  flows  more  rapidly.  However,  the  particles  also 
move  faster  so  that  the  trajectory  of  the  particles  remains  constant. 

Therefore,  the  particles  are  precipitated  in  approximately  the  same  location  j 

on  the  precipitation  tube  independent  of  temperature.  1 

i 

Wlien  all  of  the  oil  in  the  reservoir  runs  out,  air  passes  the  j 

reference  poistion  causing  a sudden  increase  in  the  amount  of  light  sensed  I 

by  the  reference  photo  resistor.  This  serves  as  a signal  to  turn  the  | 

flushing  pump  on  and  the  system  then  automatically  recycles  itself.  At  1 

normal  operating  oil  temperatures  of  121°C  (250°F),  the  system  will  j 

recycle  approximately  once  every  10  minutes.  At  lower  temperatures  the  j 

time  increases  in  proportion  to  the  viscosity.  ^ 

i 

The  percent  area  covered  by  the  particles  at  the  end  of  each  cycle  j 

is  proportional  to  the  density  of  particles  having  a selected  size.  During  | 

the  flushing  cycle  the  particles  which  have  been  precipitated  are  returned 
to  the  engine's  oil  tank  from  which  they  came.  ( 

The  ability  of  the  engine  mountable  Ferrograph  to  detect  bearing  and 
gear  failures  was  evaluated  by  rig  testing.  Ten  bearing  tests  were  con- 
ducted. Bearing  failures  were  generated  by  loading  the  bearings  such  as  ' 

to  achieve  failures  in  (1)  contact  fatigue  spalling;  (2)  wear  leading  to 
contact  fatigue  spalling;  and  (3)  excessive  wear.  Tests  were  also 
conducted  to  simulate  various  levels  of  rubbing  wear  experienced  in  gear 
tooth  contact  using  a test  rig  called  the  Geared  Roller  Test  Machine  (GRTM) . 

Data  analysis  indicated  that  the  Ferrograph  is  capable  of  detecting 

abnormal  wear  in  each  case  tested.  ■ 

After  rig  testing,  the  unit  was  mounted  on  the  TF34  engine  to  evaluate  J 

its  ability  to  function  in  a test  cell  environment.  Certain  deficiencies  ; 

in  the  unit  were  revealed,  the  major  ones  relating  to  false  readings 
obtained  as  a result  of  oil  aeration  and  engine  vibrations.  Trans-Sonlcs 
made  fixes  to  the  unit  to  correct  these  deficiencies.  Further  testing 
on  the  TF34  was  aborted  because  an  engine  oil  leak  caused  oil  seepage  into 
the  unit  and  impaired  its  ability  to  make  accurate  readings.  As  a result,  ^ 

no  reliable  data  was  taken  during  the  engine  test.  J 
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3.  Vibrating  Cylinder  Pressure  Tranducer 
a.  Description 

The  vibrating  cylinder  pressure  transducer  (VCPT)  is  a 
relatively  new  type  of  transducer  having  some  inherent  advantages  over 
other  types  (reference  6).  This  diaphragmiess  pneumatic  pressure  tra.is- 
ducer  utilizes  natural  vibration  modes  of  a metal  cylinder  to  generate 
a frequency  output  as  a function  of  applied  pneumatic  pressure.  Its 
natural  frequency-modulated  output  makes  it  immune  to  noise  and  drift 
to  which  amplitude  based  transducers  are  susceptible. 

The  operating  principles  can  be  seen  from  Figure  3,  where  a sim- 
plified sketch  of  the  VCPT  is  shown.  The  vibrating  cylinder  is  energize' 
by  the  driving  coll  electromagnet.  The  driving  coil,  in  turn,  Is  driven 
by  the  electronic  amplifier.  The  amplifier  has  a positive  feedback  loot, 
causing  it  to  be  inherently  unstable  and  oscillatory.  When  the  VCPT 
amplifier  is  first  turned  on,  electronic  noise  in  the  amplifier  is 
amplified  and  sent  to  the  drive  coil,  causing  minute,  random  ex"itatiori 
of  the  vibrating  cylinder.  The  vibrating  cylinder  will,  in  turn,  oscxll' 
at  its  natural  frequency.  This  oscillation  is  picked-up  by  the  sensing  coil 
and  transmitted  to  the  amplifier,  where  further  amplification  occurs. 

Thus,  the  osciallation  amplitude  rapidly  builds  up  until  the  limiting 
value  (set  by  amplifier  electronics)  is  reached.  The  output  of  the 
amplifier  can  be  passed  through  a pulse  shaping  circuit  and  monitored  on 
a counter.  The  natural  frequency  of  the  vibrating  cylinder  is  a functior, 
of  its  wall  tension.  As  the  applied  pressure  changes,  the  wall  tension 
varies,  causing  the  natural  frequency  to  change.  In  this  way,  the  frequen, 
output  of  the  transducer  is  a function  of  applied  pressure. 

Two  other  factors  also  contribute  to  the  wall  tension  of  th.e  cylinder; 
the  temperature  and  the  gaseous  composition  of  the  applied  pressure 
(reference  7) . Temperature  variations  change  the  elastic  modulus  and 
linear  dimensions  of  the  metal  cylinder,  causing  a change  in  cylinder 
frequency.  The  effects  of  temperature  can  be  minimized  by  using  low  temp- 
erature coefficient  type  metal  alloys.  In  addition,  this  effect  can  be  ! 

eliminated  by  monitoring  the  temperature  of  the  cylinder  via  suitable  1 

sensing  elements  imbedded  inside  the  transducer  to  provide  proper  compen-  i 

sation  to  output  data.  The  effect  of  gas  composition  is  due  to  the  ‘ 

molecular  weight  of  the  gas  species.  As  the  cylinder  vibrates  near  its  < 

natural  frequency,  the  gas  molecules  will  move  with  the  cylinder  wall, 
thereby  adding  mass  to  the  cylinder.  This  effect  will  vary  with  dificT 
gas  species,  being  more  significant  with  heavier  molecules.  This  effect 
can  be  minimized  by  making  the  cylinder  wall  thick.  However,  this  neces- 
sitates a design  compromise  as  the  Increased  thickness  adversely  affects 
the  sensitivity  of  the  transducer  (reference  7) . 

Unlike  conventional  transducers,  the  frequency  output  of  the  VCPT  is 
a non-linear  fuctlon  of  applied  pneumatic  pressure.  The  exact  form  of 
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the  function  is  a polynomial  curve.  The  coefficients  of  the  polynomial 
are  determined  for  each  transducer  through  individual  calibrations  and 
applying  computer  curve  fitting  techniques  to  the  new  calibration  data. 

The  temperature  effect  can  be  compensated  by  varying  the  calibration 
temperature  and  obtaining  a family  of  temperature  correction  curves. 

b.  Test  Results:  A model  PT-050S-1D  transducer,  manufactured 

by  Hamilton  Standard  was  evaluated  at  NAPTC  during  TEDD  Phase  III.  This 
transducer  is  an  absolute  pressure  unit  with  a factory  calibrated  range 
of  0-350,000  ncwtons/m2  (0-30  psla).  All  pressures  were  calculated  by  a 
polynomial  equation  and  temperature  correction  curve  supplied  by  the 
manufacturer.  The  VCPT  was  first  checked  against  a laboratory  pressure 
standard  (^0.07%  of  reading).  After  calibration,  the  transducer  was 
mounted  on  a J-52  engine  for  endurance  evaluations.  The  VCPT  was  mounted 
on  the  compressor  section  of  the  engine  and  measured  turbine  discharge 
pressure  through  an  in-line  air  filter.  During  more  than  200  hours  of 
engine  testing,  the  VCPT  output  was  monitored  and  checked  against  other 
standard  engine  Instrumentation  measuring  the  same  pressure  station. 

Finally,  after  the  engine  tests,  the  VCPT  was  again  calibrated  against 
the  same  laboratory  standard  to  determine  any  observable  shifts  in  trans- 
ducer characteristics.  Throughout  the  entire  test  program,  no  adjust- 
ments were  made  on  the  VCPT. 

During  laboratory  calibration  checks  before  (pre-cal)  and  after 
(post-cal)  the  engine  test,  the  transducer  accuracy  was  determined  to  be 
better  than  0.2%  full  scale.  Shown  in  Figure  4 is  a plot  of  percent  full 
scale  difference  between  the  laboratory  standard  and  the  VCPT,  versus 
pressure.  Referring  to  the  Figure,  it  can  be  seen  that  the  VCPT  accuracy 
based  on  manufacturer  supplied  coefficients  is  within  that  of  the  labora- 
tory standard  for  pressures  up  to  approximately  240,000  newtons/m2  (35  psia). 
Although  the  VCPT  readings  were  consistently  lower,  the  lower  reading  of 
the  VCPT  is  probably  due  to  the  differences  in  gas  composition  between  the 
manufacturer’s  calibration  air  supply  and  that  of  NAPTC.  This  effect  be- 
comes more  significant  at  higher  pressures  because  more  gas  molecules 
fill-up  the  cylinder  volume,  thereby  adding  significant  "effective"  mass 
to  the  cylinder  wall.  However,  it  should  be  noted  from  Figure  4 that  the 
repeatability  of  the  transducer,  even  after  more  than  200  hours  of 
endurance  testing,  was  well  within  the  uncertainty  limit  of  the  laboratory 
standard  (+0.07%  of  reading).  Thus,  if  polynomial  coefficients  were 
used  for  the  calibration  gas,  an  accuracy  of  better  than  + 0.07%  would 
be  obtained. 

During  engine  tests,  the  VCPT  output  was  checked  against  an  engine 
control  room  gage  (+1%  full  scale)  and  the  NAPTC  computer  data  acquisi- 
tion system  (4^.3%  full  scale).  All  (calculated)  VCPT  pressures  agreed 
to  well  within  the  accuracy  limit  of  the  two  instrumentation  systems. 

In  summary,  the  following  advantages  and  possible  drawbacks  relative 
to  conventional  diaphragm-type  transducers  can  be  concluded  from  NAPTC 
evaluations.  First  the  advantages: 
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(1)  Ease  oi  operation:  There  were  no  calibration  adjustments 

necessary  throughout  the  engine  test  program. 

{,2)  Accuracy  and  stability:  Throughout  the  test,  the  (cal- 

culated; pressuie  data  from  the  VCPT  agreed  to  well  within  the  accuracy 
range  of  the  stanuara  engine  Instrumentation  system.  The  transducer 
performed  equally  weti  at  the  end  ot  the  program  as  in  the  beginning. 

Its  evaluaieU  a"curacy  was  better  than  0.2%  full  scale  and  is  capable  of 
better  than  0.i;7%  ot  reading  if  gas  composition  differences  are  elimi- 
nated ot  compensated. 

(3)  Kepea tai.'i  1 ity : The  rtpeatability  of  the  transducer  was 

witliLn  the  uncertainties  of  the  laboratory  standard. 

Possible  drawbacks  are  as  follows: 

(1)  Data  reduction  complexity:  The  need  for  polynomial  eval- 

uations and  temperature  compensation  required  the  use  of  a computer  for 
real-time  pressure  readouts.  However,  linearization  and  temperature 
compensation  circuits  have  already  been  manufactured  to  provide  direct 
read-out  pressure;  this  should  remove  this  drawback  of  the  transducer. 

(2)  Response  Time:  The  finite  "fill"  time  necessary  for  the 

cylinder  volume  in  the  transducer  inherently  limits  its  response  time. 
For  example,  for  the  350,000  newton/m^  (50  psia)  VCPT  evaluated,  the 
manufacturer  has  quoted  a rise  time  of  approximately  10  milliseconds. 

(3)  Sensitivity  to  gas  composition:  For  extreme  accuracy  appli- 

cations claimed  possible  with  the  VCPT  (better  than  0.015%  of  full 
scale),  the  gas  composition  of  the  unknown  pressure  must  be  similar 

to  the  calibration  pressure  source.  For  example,  for  an  aircraft  flying 
from  sea  level  to  high  altitude  at  25°C  (77°F)  and  measuring  a pressure 
of  350,000  newton/m2  (50  psia),  the  accuracy  of  the  VCPT  could  vary  by 
as  much  as  0.037%,  using  the  same  calibration.  This  is  due  to  the 
difference  in  relative  humidity  between  sea  level  and  high  altitude. 

4.  Oil  Debris  Monitor 


As  previously  discussed,  the  progression  of  oil  system  component 
(gears  and  bearings)  failures  generate  large  metal  particles  in  the  oil. 
In  some  diagnostic  system  applications  it  will  be  most  cost  effective  to 
utilize  a relatively  small,  lightweight,  simplified  real-time  airborne 
oil  system  monitor.  The  K West  Model  0252  debris  monitor  is  a candidate 
for  this  type  of  application. 

This  monitor  utilizes  a sensing  grid  woven  in  a unique  pattern. 
Transparent  polyester  filaments  wnich  serve  as  Insulators  are  Interwoven 
with  stainless  steel  conductor.  These  elements  are  locked  together 
with  a third  small  diameter,  stainless  steel  wire  which  greatly  improves 
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the  contact  surface  when  debris  is  impinged  on  the  grid  by  the  oil  flow. 

A properletary  method  of  interconnecting  each  conductive  strand  of  the 
screen  allows  the  detection  of  conduct'' ve  debris  as  it  shorts  out 
adjacent  wire  pairs.  A solid  cone  at  the  downstream  end  of  the  screen 
deflects  oil  through  the  screen.  As  debris  collects  along  the  surface 
of  the  screen,  the  flow  will  gradually  divert  itself  to  the  remaining 
open  area  and  thus,  randomly  distribute  the  debris  over  the  entire  area 
of  the  sensor.  Any  debris  that  is  electrically  conductive  (both  ferrous 
and  non-ferrous)  will  be  detected.  The  effective  reduction  of  electrical 
resistance  as  buildup  occurs  is  read  out  on  the  data  system.  Figure  5 
shows  the  equivalent  circuit. 

This  unit  has  been  installed  on  three  TF30-P-408  engine  tests,  TEDD 
Phases  I and  II,  and  on  a TF30-P-41ZIA  engine  undergoing  afterburner 
lighting  tests.  The  unit  was  installed  in  the  scavenge  oil  line  for 
main  bearings  no.  4,  4 1/2,  5 and  6 on  the  TF30-P-412A  engine  for  a 
total  of  65  hours  test  time.  There  was  no  change  in  the  output  reading 
during  the  test.  Analysis  of  material  found  on  the  screen  after 
completion  of  the  test  showed  only  a minute  quantity  of  debris  (2.0 
milligrams)  (7  x 10“5  ounces),  consisting  mostly  of  fine  grit  clay  with 
a few  carbon  particles  and  one  long  magnetic  particle.  Thus,  no 
unsatisfactory  operation  has  been  noted  with  the  K West  unit  during  any 
of  the  TEDD  phases,  but  it  has  not  been  subjected  to  an  engine  problem 
situation  where  its  effectiveness  could  be  evaluated. 

5.  Hot  Section  Analyzer 

An  important  element  in  any  diagnostic  system  is  the  detection 
of  operational  limit  exceedances,  especially  overtemperatures  and  over- 
speeds. When  those  parameters  necessary  to  check  these  limit  exceedances 
are  available,  additional  information  about  hot  section  component  status 
is  available.  When  engine  maintenance  is  to  be  performed  "on-condltion" 
rather  than  at  fixed  time  intervals,  it  is  necessary  to  utilize  this 
information  to  determine  the  need  for  hot  section  inspection  or  parts 
replacement.  Howell  Instruments  has  developed  a Computer  Recorder  which 
attempts  to  interpret  this  Information  and  provide  maintenance  guidance. 

A Howell  Instrument  H158A-15  Computer  Recorder  (Figure  6)  was  in- 
stalled on  a TF30-P-412A  engine  which  was  undergoing  a 150  hour  accele- 
rated cyclic  endurance  test  at  NAPTC.  This  unit,  also  referred  to  as  a 
Time  Temperature  Recorder  and  Integrator  (TTRI),  has  inputs  of  turbine 
inlet  temperature  (Indirectly  measured)  and  high  rotor  speed,  and  provides 
a cumulative  Hot  Section  Factor  Count  based  on  the  integrated  time-terapera- 
ture  exposure  and  the  hot  section  metallurgical  characteristics  of  the 
engine.  It  also  displays  Hot  Section  Factor  Time  (1015°C  (1900°!)  and 
above),  time  above  1177“C  (2150“F),  total  engine  time,  and  maximum  temp- 
erature, and  has  three  overtemperature  flags  and  one  overspeed  flag. 
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The  TF30-P-412A  engine,  with  113.i  hours  since  new,  was  received 
from  the  Naval  Air  Rework  Facility,  Norfolk,  Virginia  (NARF  NORVA) 
following  a comp.^ete  engine  repair.  ITie  total  engine  test  time  at  NAPTC 
was  218.5  hours,  of  which  150  hours  were  cyclic  endurance  running.  Five 
different  cycles  were  used  during  the  endurance  test.  Test  cycle 
profiles  were  based  on  the  mociified  Qual  if  1 .:ation  test  of  an  F401-PW-400 
engine.  To  :dd  to  the  severity  of  the  test,  the  turbine  inlet  tempera- 
ture was  adjusted  to  obtain  the  maximum  allowable  value  of  1177°C 
(2150°F).  To  determine  the  ainount  of  distress  in  the  engine  hot  section, 
three  periodic  hot  section  inspections  were  made,  and  at  the  end  of 
the  test,  a complete  engine  teardown  and  analytical  inspection  was  per- 
formed at  NARF  NOKVA,  This  test  program  was  reported  on  by  reference  (8). 

The  Howell  unit  was  install'd  after  the  first  10. V hours  of  testing 
and  remained  in  use  for  the  duration  of  the  test.  Readings  were  recorded 
periodically  and  the  Instrument  was  zeroed  out  after  each  reading,  except 
for  engine  time,  which  has  no  external  zero  adjustment.  At  the  conclusion 
of  the  test,  the  following  total  readings  were  recorded: 

Maximum  Temperature  1231°C  (2248°F) 

Engine  Hours  234.1 

Hot  Section  Factor  Count  241,704  units 

Hot  Section  Factor  Time  77.5  hours 

Time  Over  1177°C  (2150°F)  21.4  hours 

Over temperature  Flag  A (1177°C,  10  second  delay)  had  been 
tripped  on  41  of  the  53  readings  taken. 

Overtemperature  Flag  B (1215“C,  10  second  delay).  Overtemperature 
Flag  C (1250°C,  10  second  delay),  and  the  Overspeed  Flag 
(103.7  % N2)  were  not  tripped. 

It  is  seen  that  there  was  a very  large  amount  of  overtemperature 
operation,  mostly  due  to  the  high  engine  trim  setting,  but  there  were 
no  events  extreme  enough  to  trip  the  B Flag. 

The  total  engine  time  recorded  (234.1  hours)  was  higher  than  the 
actual  engine  operating  time,  during  which  the  unit  was  installed 
(207.6  hours).  This  problem  was  unique  to  this  installation  and  test 
program.  The  Howell  unit  counts  engine  time  when  the  temperature  reaches 
a predetermined  level,  in  this  case  200°C  (392“F)  turbine  inlet  tempera- 
ture, provided  power  is  on  to  the  unit.  During  the  test,  the  TF30-P-412A 
turbine  inlet  temperature  was  still  between  29C'’C  (555°F)  and  315°C 
(600°F)  after  completion  of  rundown,  and  power  renuilned  on.  Since  there 
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were  approximately  500  starts  and  shutdowns  made,  the  extra  time 
amounted  to  a significant  amount.  This  is  not  a problem  in  the  field, 
where  there  is  no  power  to  the  unit  when  the  engine  is  shut  down,  nor 
did  this  invalidate  any  of  the  other  data  obtained  during  this  test. 

Figure  7 shows  a plot  of  the  Hot  Section  Factor  (HSF)  count  versus 
engine  time  (from  Howell  Unit).  The  relative  temperature  exposure 
severity  of  different  phases  of  the  test  can  be  seen  from  the  different 
rates  of  increase  in  HSF  count.  The  times  of  hot  section  inspections 
are  indicated,  along  with  the  parts  replaced.  Reference  8 provides 
more  details  on  engine  inspection  and  teardown  results.  The  signifi- 
cance of  a HSF  count  must  be  determined  by  correlating  it  with  the  need 
for  maintenance  action.  This  is  done  by  comparing  the  count  level  with 
actual  engine  hot  section  condition  over  a large  enough  sample  of  engines 
This  has  not  been  established  for  the  TF30-P-412A  engine. 

In  addition  to  being  an  indicator  of  hot  section  condition,  this 
type  of  unit  has  been  found  to  provide  other  maintenance  benefits  when 
used  on  the  F-8/J57,  revealing  the  existence  of  bad  trim  boxes, 
incorrect  trim  procedure,  short  throttle  rods,  incorrectly  wired 
thermocouple  harnesses,  etc.,  and  documenting  overtemperature  experi- 
ences (reference  9). 

Whatever  the  overall  maintenance  benefits  of  this  type  of  unit,  it 
is  necessarily  limited,  since  it  is  a single  parameter  system.  Future 
hot  section  diagnostics  should  take  more  factors  into  account,  such  as 
low  cycle  fatigue,  thermal  shock,  and  performance  degradation.  It  is 
noted  that  Howell  has  developed  a more  sophisticated  system,  called. 

Jet  Engine  Monitor  (JEM),  but  this  has  not  been  fully  evaluated  yet. 

6.  Vibration  Analyzer 

Vibration  analysis  is  an  asset  to  a diagnostic  system.  A vib- 
ration analyzer  coupled  with  well  placed  vibration  transducers  can 
provide  an  indication  of  mechanical  health  not  available  from  other 
sources  as  well  as  provide  second  source  information  to  confirm  other 
analysis.  The  vibration  analyzer  used  in  the  TEDD  program  is  capable 
of  examining  all  areas  of  an  engine,  from  the  gearbox  to  the  main 
rotor  bearings. 

During  the  running  of  Phase  I and  Phase  II  of  the  Turbine  Engine 
Diagnostic  Development  program,  the  General  Electric  engine  vibration 
analyzer  was  used  and  found  to  give  unsatisfactory  readings  for  gear- 
gross  experiments  (reference  2) . To  determine  the  reason  for  these 
erroneous  readings,  the  data  from  the  TEDD  program  v/as  rerun  and  a 
detailed  analysis  was  performed  on  the  vibration  analyzer  in  Phase  III. 

The  technique  used  to  do  gear-gross  analysis  is  shown  in  the  block 
diagram  of  Figure  8.  The  tachometer  signal  is  conditioned  by  multi- 
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plier-dividers  co  provide  a clock  signal  for  the  digital  comb  filter 
and  the  synchronous  detector.  Lhe  digital  comb  filter,  a 256  section 
re.,u  time  enhancement  filter,  extracts  the  gear  mesh  frequency  of 
interest  and  its  associated  harmonics.  If  there  are  gross  gear  defects, 
such  as  a worn  shaft  or  bearings  or  unbalance,  the  gear  mesh  frequency 
will  be  amplitude  moiiu'.ated  by  the  once  per  revolution  frequency  of  the 
gear.  In  real  time,  tbLs  signal  will  appear  as  shown  in  Figure  8b.  The 
synchronous  detector  is  used  to  extract  the  modulating  frequency.  In 
the  synchronous  detector,  the  gear  mesh  frequency  or  carrier  frequency  is 
rein  .ved  and  only /the  modulation  Is  left.  The  amplitude  of  the  modulation 
is  an  indication  of  gross  defecc.s. 

luitlaily,  tne  pronlera  with  gi ar-gross  results  was  thought  to  be  an 
unstable  tacnometer  signal  caused  by  instantaneous  speed  changes  in  the 
gear  i.rai  n due  to  backlash  and  pi  iv  in  the  gear  mesh. 

This  torsional  piav  results  in  a frequency  modulated  tachometer 
signal  that  cTeates  a problem  in  the  analyzer  by  causing  rapid  phase 
shif.s  through  the  multiplier-divider  circuits.  These  circuits  have  a 
transient  response  sufficient  to  track  an  accelerating  engine;  however, 
they  require  a finite  amount  of  time  to  settle  to  a phase  locked  con- 
dition. In  order  to  do  gear-gross  analysis,  five  multiplier  phase- 
locked  loops  a|:e  used  to  condition  the  tachometer,  and  it  was  thought 
that  the  rapid ' phase  shifts  were  propagated  through  the  multipliers, 
causing  the  defection  circuit  to  saturate. 

To  analyze  the  synchronous  detector  circuit,  a signal  from  a function 
generator  was  used  for  both  tachometer  and  sensor  inputs  to  insure  that 
the  signal  was  stable.  The  multiplier  circuits  were  set  to  one  over  one 
ratio  and  the  signal  was  traced  through  the  various  stages  of  condition- 
ing. It  was  found  that  the  detector  showed  modulation  on  the  signal, 
while  there  was  none  on  the  input.  It  was  found,  at  this  point,  that 
there  were  some  switching  incongruities  from  the  digital  comb  filter 
that  were  being  amplified  by  the  automatic  gain  control  used  in  the 
synchronous  detector. 

The  nature  of  the  digital  comb  filter  is  such  that  the  switching 
spikes  are  evident  in  the  output  hut  can  be  smoothed  out  with  proper 
filtering.  The  digital  comb  filter  switches  the  input  signal  to  a 
series  of  capacitors  in  sequence  and  at  a rate  determined  by  the 
multiplier-divider  settings.  A simplified  diagram  is  shown  in  Figure 
9a  and  the  resultant  output  in  Figure  9b.  The  transitions  between 
capacitors  is  the  noise  that  was  being  passed  to  the  synchronous 
filters,  generating  erroneous  gear-gross  readings. 

Wlien  this  noise  was  filtered  out  just  before  the  automatic  gain 
control,  this  produced  satisfactory  results  on  the  synthetic  signal. 

The  taped  data  from  the  TEDD  program  was  then  played  into  the  analyzer. 
The  results  are  presented  in  Table  I,  along  with  the  original  readings 
from  the  Phase  11  engine  test.  The  first  set  of  readings  were  taken 
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at  lower  engine  speed  (Idle)  with  good  correlation  obtained  for  gear 
local  and  much  better  results  for  gear-gross  experiments.  At  the 
higher  engine  speed  (Intermediate),  however,  gear-gross  experiments  were 
mostly  over  limits.  This  problem  is  the  result  of  much  more  frequency 
modulation  on  the  tachometer  signal,  and/or  poor  communication  between 
the  object  being  analyzed  and  the  accelerometer. 

In  addition  to  the  speed  signal  changing  phase,  the  gear  passage 
frequencies  will  also  be  frequency  modulated  and  perhaps  add  to  the 
problems  encountered  ii.  the  gear  gross  defects  experiments.  This 
situation  can  be  analyzed  only  with  an  unmodulated  tachometer  signal. 
Additional  vibration  analyzer  analysis  should  be  done  with  data  that 
iias  a tachometer  signal  with  less  than  .8%  modulation  to  resolve  the 
question  of  frequency  modulation  versus  lack  of  communication  to  provide 
valid  gear-gross  experiment  results. 

7 . Speed  Sensor 

In  an  attempt  to  improve  the  performance  of  the  vibration 
analyzer  on  overall  gear  defect  experiments,  an  investigation  was  done 
on  various  speed  sensing  techniques.  It  is  believed  that  part  of  the 
problem  with  the  gear-gross  defects  is  due  to  an  unsteady  tachometer 
signal  caused  by  torsional  effects  and  gear  meshing  tolerances.  It 
is  also  possible  that  the  various  loads  on  the  gear  box  are  not  steadily 
applied  and  may  affect  the  rotation  of  the  tachometer  pad.  There  are 
several  places  on  an  engine  where  speed  data  can  be  collected.  These 
include  a gear  box  pad,  the  rotor  shaft  itself  in  a nose  cone  pickup, 
or  from  a compressor  or  turbine  stage.  Below  is  a presentation  of 
data  recorded  and  evaluated  from  tach  pickups  on  the  gear  box,  from  a 
gear  reduction  from  the  rotor  shaft  and  directly  from  the  shaft.  Also 
presented  is  the  variation  of  a calibration  signal,  recorded  and  re- 
produced, which  can  be  used  as  an  Indication  of  the  error  Introduced 
by  tape. 
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Ic  (_an  ba  si'  *n  tnaL  tnere  is  a^ipi  oximately  2%  jitter  on  the  gear 
box  pickup  at  aighei  rpni.  This  amouat  is  intoierable  for  vibration 
anaiyzei  exper iiuei'.t a . it  was  noted  that  variations  to  approximately 
1%  could  be  tolerated  by  the  vibration  analyzer. 

Bkice  passage  J<ua  ..as  not  analyzed  quaatilative.iy;  however,  it  is 
modulate ..  i(  3s  tliaii  ei  t seen  on  data  t rom  the  gear  reduction  in  the 
front  ae-  of  the  ir-  .d,  aad  would  prove  to  be  satisfactory  for  the 
anai . zer . 

It  is  ; e coioraendeo  that  ’oc  tui'.rt  ibra’.  ion  analyzer  use,  a blade 
passage  pickup  be  use.'  lor  the  t;i.d;  sp.-ed  rotor  and  gear  box  experiments. 

o.  : adiatlon  Pyrometer 


file  bo  tar  radiation  pyrometers  used  on  TEDU  Phase  II  produced 
results  that  required  'urther  analysis.  The  readings  from  the  pyro- 
meters for  peak  and  average  peak  showed  higher  values  than  calculated 
gas  temperature.  An  (.scilloscope  was  used  to  look  at  the  analog 
signal  troiu  the  pyroi..uter,  and  it  was  noted  that,  in  addition  to  the 
blade  passage  signal,  spurious  spikes  were  seen  by  the  pyrometer. 

Aside  from  the  noise  problems,  the  aperture  probe  showed  a decreasing 
blade  temperature  versus  time  for  comparable  engine  settings,  indicating 
a possible  sooting  condition.  When  the  pyrometer  was  removed,  a post 
calibration  was  run  and  a degradation  was  noted.  Attempts  to  clean 
the  probe  ir.iproved  its  sensitivity  to  some  extent.  Upon  shipment  to 
Solar,  it  was  found  to  be  inoperative.  Solar  disassembled  the  probe 
and  found  that  the  fiber  optics  had  broken.  There  was  also  an  oily 
deposit  in  the  purge  air  tube  that  Indicated  that  possibly  the  purge 
had  been  dirty.  This  did  not  seem  to  be  a problem  with  the  lens  probe. 
These  pr<,biems,  particularly  the  noise  in  tne  analog  output  signal,  are 
the  ones  Investigated  during  this  phase  of  TEDD. 

In  discussing  pyrometry  with  other  manufacturers  (at  first  primarily 
with  Detroit  Diesel  Allison,  Division  of  General  Motors,  since  they 
had  one  of  tlie  early  probes  developed  bv  Solar),  it  was  found  that 
others  had  seen  the  type  oi  signal  output  observed  during  the  TEDD 
program.  in  investigations  at  oth *r  facilities,  using  different  in- 
staiiation  locations,  it  was  found  that  the  noise  could  be  reduced, 
or  eJlmlnated,  by  shieiding  the  probe  from  the  flame.  This  has  been 
done  by  looking  at  the  back  side  of  the  blade  or  radially  to  the  blade 
and  choosing  angles  such  tiiat  the  probe  cannot  view  the  flame,  either 
directly  or  through  reflections.  Tills  approT-h  has  resulted  in  the 
best  reduction  of  spurious  noise,  partiiuiarly  on  engines  running  at 
current  turbine  inlet  temperatures. 

Another  approach  to  measuring  tlie  average  and  peak  average  blade 
timperature  was  to  enliance  the  signal  tlir.nigh  techniques  used  for 
vibration  analysis.  Data  was  obtained  from  Allison  that  looked  very 
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similar  to  that  seen  in  the  TEDD  program.  It  was  played  into  the 
vibration  analyzer  and  processed  through  the  digital  comb  filter  to 
extract  the  spurious  signals.  Figure  10  shows  scope  traces  of  both 
the  unconditioned  signal  and  the  conditioned  signal.  It  can  be  seen 
that,  in  general,  the  digital  comb  filter  does  a reasonably  good  job 
of  extracting  the  blade  passage  frequency.  There  are  some  cases, 
however,  where  the  input  electronics  were  saturated  and  at  these  points, 
there  is  no  blade  passage  frequency  to  extract. 

The  data  analyzed  by  the  vibration  analyzer  was  not  calibrated  to 
provide  temperature  information,  and  it  was  coupled  through  the 
analyzer  such  that  an  overall  temperature  could  not  be  recovered. 
However,  some  of  the  analyzer's  detection  circuits  could  be  used  to 
automatically  flag  hot  blades.  The  impact  index,  the  ratio  of  peak  to 
average  signal  level,  could  be  used  to  detect  blades,  or  a blade  that 
was  significantly  hotter  than  the  average.  Separate  circuitry  would 
be  required  to  provide  the  overall  average  blade  temperature  level. 

A second  technique  that  is  being  developed  by  Solar  and  Allison  is 
to  detect  the  incoraming  spikes  and  reject  them  above  a prefixed  level, 
usually  a not-to-exceed  gas  temperature.  This  seems  to  work  fairly 
well,  but  was  not  analyzed  at  NAPTC.  The  disadvantage  of  this  system 
is  that  during  high  temperature  running,  the  signal  could  conceivably 
be  rejected  more  often  than  it  would  be  measured. 

9 . Internal  Accelerometer 

In  order  to  monitor  bearing  health  more  accurately  than  has  been 
done  in  the  past,  it  has  been  recommended  that  internal  accelerometers 
should  be  mounted  on  the  bearings  (reference  2).  There  are  some  common 
problems  associated  with  this  task  that  have  been  under  design  con- 
sideration for  several  years.  These  problems  are  the  environmental 
restrictions  on  the  accelerometer,  signal  lead  routes,  and  bearing 
housing  accessibility,  which  places  restrictions  on  the  design. 

Accelerometers  are  relatively  new  devices,  coming  into  full  develop- 
ment in  the  last  five  to  ten  years;  and  have  been  considered  delicate 
transducers.  The  piezoelectric  material  used  as  the  transduction 
element  is  sensitive  to  temperature  and  has  an  upper  limit  dependent 
on  the  type  of  material  used.  The  various  housing  and  lead  coupling 
techniques  used  in  accelerometers  are  not  rugged  in  construction.  Due 
to  the  characteristic  of  the  piezoelectric  material,  very  low  level 
signals  are  produced;  thus,  careful  consideration  must  be  given  to 
the  manner  in  which  the  signal  leads  are  run. 

Current  state-of-the-art  accelerometers  are  fairly  rugged  in  that 
they  can  withstand  temperatures  to  500°C  (932°F),  and  the  leads  are 
integral  to  the  case  so  that  connections  on  the  accelerometer  are  not 
required.  An  accelerometer  of  this  type  was  used  for  the  internal 
accelerometer  in  the  TEDD  program. 
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The  second  problem  for  internal  accelerometer  placement  is  how  to 
get  the  signal  to  the  outside  of  the  engine.  If  it  is  assumed  that 
the  engine  will  be  disassembled  for  repair  of  other  than  the  accelero- 
meter, the  signal  lead  should  not  cross  engine  boundaries  that  would 
be  separated.  In  some  engine  installations  it  is  possible  to  mount 
the  accelerometer  in  a cooled  location,  and  then  it  becomes  necessary 
to  trace  the  leads  through  extremely  hot  sectio-is.  This  is  not  as 
serious  a problem  as  is  the  -’C'-eleTineter ' s installation;  however,  it 
must  he  considered. 

The  p;vipiem  of  limited  lift  for  the  accelerometer  makes  it  neces- 
sar>  that  it  be  acce:  ible  from  the  engine  case.  Several  designs  have 

been  advocated  to  permit  this  accessibility.  me  concept  is  to  place 
the  accelerometer  at  the  bearing  through  a "dipstick”  approach;  that 
is,  a rod,  with  the  signal  lead  integral  with  it  and  the  accelerometer 
at  its  end,  is  screwed  into  the  mount.  This  requires  a suitable  engine 
design,  but  places  the  accelerometer  and  lead  in  a position  for 
accessibility. 

This  approach,  though  relatively  straight  forward,  runs  into 
problems  when  the  actual  bearing  design  is  considered.  In  dual  rotor 
engines  and  engines  of  advanced  design,  some  of  the  hearings  are 
designed  to  be  soft  mounted.  Thus,  mounting  is  affected  by  a canti- 
levered design  that  permits  the  bearing  and  bearing  housing  to  float 
in  the  engine.  The  soft  mounted  bearing  is  not  conducive  to  an 
accessible  accelerometer  installation. 

This  type  of  bearing  mount  was  used  on  the  TF30-P-408  rear 
bearing  (TEDD  Phase  II).  The  accelerometer  was  hard  mounted  to  the 
bearing  housing  and  was  only  accessible  by  removing  the  exhaust 
section  and  by  entry  into  the  sump  oil  housing. 

In  general,  these  problems  were  the  ones  considered  by  Pratt  & 
Whitney  when  installing  internal  accelerometers  on  the  F401  engine. 

The  bearings  in  which  they  were  Interested  were  the  if2  and  #3  bearings. 
The  //3  bearing  is  soft  mounted  and  a dipstick  approach  was  used  for 
accelerometer  accessibility.  It  is  their  assessment  that  the  infor- 
mation obtained  was  not  pertinent  to  the  bearing  problems  due  to  lack 
of  coramxinication  between  bearing  and  accelerometer  and,  therefore,  it 
has  been  recommended  that  the  accelerometer  not  be  used  in  F401 
production  engines. 

B.  Software 


Data  Smoothing  Techniques 


In  order  to  perform  accurate  and  valid  diagnostic  logic,  the 
parameters  used  in  the  logic  must  truly  represent  the  state  of  the 
engine.  Evaluation  of  the  "truthfulness"  of  a sensed  parameter  must 
consider  the  accuracy,  repeatability  and  time  response  of  the  sensed 
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parameter  relative  to  the  actual  state  of  the  parameter.  Accuracy, 
repeatability  and  time  response  are  performance  characteristics  of  a 
sensor  which  nriy  become  involved  in  a system  effectiveness  trade-off; 
that  is,  sensors  with  the  best  performance  may  be  large,  expensive,  or 
delicate  such  that  their  usefulness  in  an  airborne  diagnostic  system  is 
limited.  Therefore,  it  will  be  necessary  at  times  to  sacrifice  sensor 
performance  for  system  suitability.  However,  it  may  be  possible  to 
"buy  back"  the  sacrificed  sensor  performance  by  applying  software  tech- 
niques designed  to  optimize  the  per fcrmance . It  is  obvious  that  nothing 
can  be  done  to  improve  tne  accuracy  of  a given  sensor,  but  software 
techniques  which  pertorm  data  validity  tests  (l.e.,  identify  failed 
sensors  and  signal  driit)  can  be  used  to  maintain  the  highest  possible 
accuracy.  Likewise,  it  is  impossible  to  improve  the  time  response  of  a 
given  sensor,  bur  software  techniques  can  be  used  to  minimize  the  time 
response  distortion  Inherant  In  some  diagnostic  logic  applications. 
Improvement  in  parameter  repeatability  can  be  directly  addressed  by 
application  of  data  smoothing  software  techniques,  and  the  following 
discussion  addresses  the  TEDD  Phase  III  effort  on  this  subject.  A 
thorough,  mathematical  discussion  of  this  effort  is  presented  in 
Appendix  A. 

In  order  to  evaluate  the  effectiveness  of  data  smoothing  techniques, 
the  software  performance  goals  or  criteria  must  be  defined.  The  following 
criteria  were  developed  and  used  in  this  study: 

a.  Storage  Required:  This  is  concerned  with  how  much  data 

storage  is  required  to  perform  the  various  smoothing  techniques  and  also 
how  much  storage  the  technique  itself  will  occupy.  Obviously,  the  less 
storage  required  the  better. 

b.  Ability  to  Adjust  Rates  of  Response:  The  nautre  of  TEDD 

data  is  such  that  it  is  not  known  whether  the  next  data  sample  will  have 
a similar  magnitude  as  the  last,  or  represent  a radical  change.  This 
characteristic  dictates  that  the  techniques  must  be  able  to  assign  weights 
to  the  data  depending  on  certain  engine  parameter  rates  of  change.  How 
well  the  techniques  can  do  this  is  an  important  criteria. 

c.  Simplicity  of  Computation:  Since  the  ultimate  application  of 

TEDD  programming  is  a real  time  situation  in  which  the  computations  are 
repeated  with  each  new  data  scan,  and  since  applied  smoothing  will  be 
essentially  an  iterative  processs,  the  amount  of  time  spent  smoothing  is 
critical.  Those  techniques  that  employ  fewer  m;ithematical  steps  and  can 
work  effectively  on  less  data  will  be  of  high  Interest.  Not  only  will 
less  memory  be  used,  but  they  will  also  be  faster. 

d.  Response : This  criteria  applies  mainly  to  transient  con- 

ditions. Due  to  mode  detection  techniques,  little  distortion  of  rates 
of  change  of  certain  key  parameters  can  be  tolerated.  The  tradeoff  here 
is  to  get  the  true  rate  of  change  of  the  parameter  by  smoothing  out  noise 
but  not  to  the  extent  of  distorting  the  true  event. 


26 


NAP TC- PE- 8 8 


e.  Accuracy  and  Noise  Reduction:  This  criteria  is  used  to 

compare  how  well  the  smoothed  data  compares  to  the  model  of  the  true 
event.  Standard  deviation  defined  as  the  measure  of  how  much  the  data 
varies  from  the  mean  of  the  sample,  will  be  used  to  judge  this  criteria. 

Two  techniques  of  data  smoothing  were  e.’aiuated  in  this  study, 
namely  moving  average  anu  exponential  averaging.  The  moving  average 
technique  is  one  of  the  methods  chosen  for  evaluation  since  it  was  used 
previously  dui Ing  lET.J  Phase  II  and,  although  it  had  some  drawbacks,  was 
fairly  successful  and  accurate  enough  to  warrant  further  investigation. 
Simply  s ated,  the  moving  average  technique  calculates  the  average  of  the 
lasi  N values  of  a given  p..;ameter.  it  minimizes  the  sum  of  the  squares 
of  the  differences  between  the  N most  recent  observations  and  the  best 
guess  for  that  average.  The  simplicity  of  computation  for  this  method  is 
an  obvious  plus  in  a TLDD  operating  environment.  The  value  that  is  picked 
for  N,  the  number  of  data  points  to  be  included  in  each  average,  will  not 
be  the  same  for  all  data.  A s.^ail  value  of  N is  usefui  for  eliminating 
small  variations  in  the  data,  but  is  relativel\  ineffective  for  large 
variations.  Conversely,  if  too  large  a value  of  N is  selected,  there  is 
the  problem  of  excessive  smoothing,  thereby  ignoring  some  of  the  true  event. 

An  inherent  problem  in  applying  this  method  is  the  lack  of  data  for 
initial  conditions  since  there  are  no  previou-  data  points  to  use  in  com- 
puting the  initial  average.  Extra  weight  will  be  given  to  the  initial 
data  points  until  N readings  are  obtained. 

Exponential  smoothing  is  a method  of  estimating  each  new  data  point 
in  a time  series  by  letting  the  new  smoothed  value  be  equal  to  a fraction 
(1  - a)  of  the  previous  smoothed  value,  plus  a fraction  (a)  of  the  new 
data  point,  a is  defined  as  the  weight  given  to  the  new  data  point. 

This  method  of  smoothing  is  simple  to  compute  and  has  the  advantage  of 
small  amounts  of  data  storage.  Only  two  words  of  storage  are  required 
for  each  parameter  instead  of  the  N words  required  by  the  moving  average 
technique.  Another  advantage  is  its  flexibility  because  only  a needs  to 
be  varied  to  change  the  rate  response.  When  a is  small,  the  method  gives 
little  weight  to  the  new  data,  which  is  good  for  smoothing  random  noise 
in  a steady  stage  condition.  A large  a will  weight  the  new  data  heavily 
so  it  responds  to  rapid  changes  in  the  time  series,  which  is  necessary 
for  tracking  transients. 

The  data  used  in  this  study  was  stored  on  magnetic  tape  during  TEDD 
Phase  II  TF30-P-408  engine  testing.  TF30-P-412  data  from  a NAPTC  Fleet 
Service  Test  Program  was  also  used.  The  parameters  which  were  used  in 
the  study  are  power  lever  angle  (PLA) , high  pressure  spool  rotor  speed 
(N2)  and  burner  pressure  (Pb).  The  values  of  N for  the  moving  average 
technique  were  1,  2,  3,  ...  10,  and  for  a in  the  exponential  averaging 
technique,  the  values  were  .1,  .2,  .3,  ...  .9.  The  detailed,  mathematical 
results  of  the  study  appear  in  Appendix  A.  The  following  summarizes  the 
results. 


27 


NAPTC-PE-88 

When  comparing  the  two  methods  of  smoothing  with  the  criteria 
developed  earlier,  it  is  obvious  the  exponential  smoothing  technique  is 
better  suited  to  TEDD  applications  than  the  moving  average.  The  fol- 
lowing chart  shows  the  comparison. 


Criteria 

Storage  Required 

Ability  to  adjust 
rate  of  response 

Simplicity  of 
computation 

Response 

Accuracy  and 
noise  reduction 


Moving  Average 

N X (#  of  parameters) 

Programming  more 
complicated 

More  instructions 
to  implement 

O.K.  for  proper  N 
O.K.  for  proper  N 


Exponential 

Smoothing 

2 words  for  each 
parameter 

easy  to  implement 

Fewer  instructions; 
therefore  faster 

Always  better 
for  tested  range 

Always  better 
for  tested  range 


To  apply  the  exponential  smoothing  to  TEDD  data,  it  is  necessary  to 
have  pre-defined  separate  a's  for  steady  state  and  transients,  depending 
on  response  and  accuracy  desired  for  each  parameter.  Upon  deciding  which 
mode  is  current,  the  appropriate  a is  applied  to  the  data.  Programming 
details  such  as  initial  conditions  cannot  be  overlooked  when  switching 
a during  a transient. 

2.  Multiple  Fault  Gas  Path  Analysis 

In  a diagnostic  system  the  gas  path  analysis  element  (GPA)  is  the 
portion  of  the  logic  which  addresses  the  aerothermal  performance  of  the 
engine  components.  GPA  quantitatively  determine  changes  in  some  or  all 
of  the  following:  Compressor  and  turbine  efficiencies,  changes  in 

turbine  and  exhaust  nozzle  areas,  bleed  and  power  extractions,  and 
changes  in  compressor  pumping  capacity.  When  one  of  the  goals  of  a 
diagnostic  system  is  to  improve  maintenance  through  the  use  of  "on- 
condition"  maintenance  rather  than  maintenance  at  fixed  time  intervals, 

GPA  provides  the  "condition"  of  the  engine  components  which  is  Information 
not  available  through  visual  or  other  inspection  techniques. 

A number  of  GPA  techniques  exist  covering  a wide  range  of  capability 
and  complexity.  Once  a GPA  technique  has  been  verified  to  perform  as 
intended,  it  is  possible  to  evaluate  the  technique  in  terms  of  capability 
and  complexity.  The  capability  of  a technique  Includes  its  flexibility, 
qualitative  performance  and  quantitative  performance.  Flexibility  addres- 
ses the  ability  of  the  technique  to  detect  the  types  of  gas  path  faults 
that  the  particular  engine  to  be  diagnosed  exhibits.  Flexibility  also 
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considers  the  measured  parameter  requirements  (are  they  measurable?)  as 
well  as  the  ability  of  the  technique  to  detect  erroneous  input  para- 
meters and  to  function  (fully  or  partially)  without  that  parameter. 
Qualitative  performance  addresses  the  GPA  technique's  ability  to 
identify  the  physical  fault  (dirty  compressor,  eroded  turbine,  etc.) 
responsible  for  the  detected  change  in  gas  path  performance.  An  import- 
ant consideration  is  the  technique's  ability  to  identify  multiple 
simultaneous  as  well  as  single  faults.  Quantitative  performance  is  the 
GPA  technique's  ability  to  calculate  the  level  of  component  deteriora- 
tion and  Lue  precision  of  chat  calcul.^tion.  Tlie  complexity  of  a GPA 
technique  involves  at  least  three  elements;  background  or  development 
requirements,  computacional  requirements,  and  computer  storage  re  ,uire- 
nents.  The  background  or  development  requirements  effect  the  lead  time 
requii ements  for  system  developmei^t  and/or  modification  as  well  as  costs. 
Computational  complexity  considers  the  type  and  number  of  calculations 
required  by  the  techniqtie.  Minimizing  the  computational  requirements 
optimizes  the  speed  of  the  technique,  which  is  particularly  critical  when 
GPA  is  to  be  performed  airborne.  Computer  storage  requirements  addresses 
the  technique's  input  parameter  storage,  storage  of  the  computation  coding, 
storage  of  constants  and  coefficients  associated  with  the  technique,  and 
output  data  storage.  For  an  airborne  application,  the  computer  storage 
requirements  are  likely  to  be  a limiting  factor. 

In  TEDD  Phase  11  and  111,  two  approaches  to  GPA  have  been  examined. 

In  Phase  II,  an  in-house  GPA  technique  was  developed  using  a fault 
signature  recognition  approach,  described  in  Reference  2.  In  Phase  III, 
a GPA  technique  developed  by  Hamilton  Standard  Division  of  United  Tech- 
nologies Corporation  known  as  Parameter  Interrelationships  was  evaluated 
using  Phase  II  TF30-P-A08  data.  The  final  report  by  Hamilton  Standard 
(Reference  10)  completely  describes  the  Parameter  Interrelationships 
methodology  and  the  results  using  TEDD  Phase  II  data.  Appendix  B of 
this  report  summarizes  the  methodology  and  results.  The  following  chart 
compares  the  two  methods. 


CRITERIA 
A.  CapAblllty 

(1)  FlcxihlUty 


(2)  'jualltacive 
Perforaant  e 


TEDD  PHASP  II  HFTHCI 


PARAMETER  INTERRELATIONSHI PS 


Fault  selection  llmlred  only  bv  the  flex- 
ibility of  the  cycle  deck  used.  Requires 
only  normal  oessurable  parameter*-.  Detect 
Inatruotentatlon  failure^  svsten  inoperative^ 


Fault  ^election  Halted  by  para- 
meter selectlon-lnvolvea  trade- 
off. Instrumentation  failure 
not  ad'^reaaed  In  thla  phase. 


Single  tault  onlv.  Provide:*  ranking  of  Multiple  fault  capabllltv. 

probable  faults. 


O)  'Juantltatlve 
Performance 

B.  Complexity 

fl)  Background 

<2)  I'onputatlon 

(^)  'torage 


Loaputea  nlngle  fault  level  and  Confidence  Coaputea  a fault  level  for  each 

Factor  poaalhle  fault. 


Extensive  cycle  deck  exercise  followed 
bv  lengthy  data  nanlpulatlon. 

Multiple  matrix  operations.  Extensive 
calculat Iona. 

Thousands  of  words  of  coefficient  stor- 
age - lengthv  program  coding. 


First  time  development  extensive. 
Completed  and  reusable. 

t>ne  matrix  aulclplicatlon. 


Less  'ban  a thousand  words  - atoT' 
age  and  progrsenlng  combined. 
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Whatever  form  it  may  take,  any  engine  diagnostic  system  must  rely 
on  discernable  changes  in  observable  parameters  in  order  to  detect 
physical  faults.  Physical  faults  in  an  engine  are  in  general  some  com- 
bination of  problems  such  as  erosion,  corrosion,  foreign  object  damage 
(F.O.D.),  worn  seals,  etc.  Although  the  physical  faults  result  in 
changes  in  observable  parameters  such  as  rotational  speeds,  temperatures, 
pressures,  etc.,  it  is  more  fundamental  to  consider  the  physical  faults 
as  causing  changes  in  engine  component  performance  parameters  such  as 
component  efficiencies,  turbine  and  exhaust  nozzle  areas,  compressor 
pumping  capacities,  etc.,  end  that  these  changes,  in  turn,  result  in  dis- 
cernible changes  in  the  measurable  engine  parameters.  If  proper  equa- 
tions can  be  estabiishea  relating  the  dependent  engine  variables  (measur- 
able quantities)  to  the  generally  nonmeasurable  independent  component 
performance  variables,  then  an  inverted  set  of  relationships  can  be 
obtained  relating  the  independent  variables  to  the  dependent  variables. 

This  could  then  form  the  basis  of  an  engine  diagnostic  system  which  would 
have  the  capability  of  detecting  multiple  faults  as  well  as  the  magnitude 
and  direction  of  the  fault. 

This  is  the  basic  technique  employed  by  Hamilton  Standard  in  their 
engine  diagnostic  system.  The  major  obstacle  in  employing  such  a diag- 
nostic scheme  is  in  establishing  the  proper  relationships  that  exist 
between  the  dependent  and  independent  engine  variables.  Two  different 
approaches  were  used  to  establish  these  relationships.  Before  discussing 
these  approaches,  it  is  appropriate,  at  this  time,  to  define  two  terms 
which  will  be  used  frequently.  The  first  term  is  "delta"  (A).  In 
everything  that  follows,  the  term  "delta"  will  be  taken  to  mean  the  percent 
of  point  difference  in  a particular  variable  relative  to  its  nominal  or 
baseline  value.  All  A's  are  dimensionless  quantities  whose  value  for  any 
variable  does  not  depend  on  the  units  used  to  measure  the  variable.  The 
second  term  is  "fault  coefficient".  This  is  just  the  percent  change  in 
a dependent  variable  per  unit  percent  chat  ge  in  an  Independent  variable. 

Note  that  since  the  deltas  are  dimensionless,  the  fault  coefficient  is 
also  dimensionless. 

Two  methods  were  used  in  this  study  to  calculate  the  required  fault 
coefficients.  The  first  method  used  an  elaborate  gas  path  analysis  of 
the  engine  developed  by  Hamilton  Standard  and  is  primarily  theoretical  in 
nature.  The  second  method  is  empirical  and  used  available  baseline  test 
data  as  well  as  engine  test  data  obtained  with  known  changes  made  to  the 
baseline  configuration. 

As  shown  in  reference  11  a general  Influence  coefficient  matrix  can 
be  derived  for  any  particular  gas  turbine  cycle.  This  matrix  represents 
a set  of  linearized  equations  which  Interrelates  the  various  dependent 
and  independent  engine  performance  parameters.  Since  the  book  was  written, 
extensive  refinements  and  modifications  were  made  to  Improve  the  precision 
of  the  coefficients  and  to  extend  the  approach  to  twin  spool  turbofan  and  i 

other  type  engines.  Each  of  the  coefficients  is  derived  by  proper  manl-  I 
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pulation  of  the  differential  forms  of  the  basic  thermodynamic  equations 
for  each  engine  typeo  Factors  such  as  component  efficiency  variations; 
specific  shapes  of  the  component  performance  maps;  engine  power  balance; 
conservation  of  mass,  energy,  and  momentum;  variable  specific  heat  effects; 
nozzle  unchoking  effects,  etc.,  are  taken  into  account.  For  the  actual 
TFiO-P-408  engine  used  in  the  test,  such  details  as  component  efficiency 
variations,  specific  shapes  of  component  performance  maps,  etc.,  were  not 
available.  Judgement  numbers  based  on  experience  and  the  model  specifica- 
tion for  this  engine  were  used  to  obtain  the  theoretical  fault  coefficients. 

An  alternate  approach  using  actual  engine  data  supplied  by  NAPTC  is 
available  tor  calculativig  tne  fault  coel i icients . Ihis  method  is 
empirical  and  makes  use  oi  nomina'  baseline  engine  data  as  well  as  steady 
state  engine  data  with  deliberate  perturbations  of  known  amounts  intro- 
duced individually.  From  test  data  where  all  independent  variables  are 
held  fixed  except  one,  the  resulting  delta  of  any  measurable  dependent 
variable  is  a result  of  the  one  known  perturbation.  In  empirically 
determining  the  lault  coefficients,  it  has  been  tacitly  assumed  that 
independent  component  performance  parameters  such  as  component  efficiencies, 
turbine  nozzle  area,  etc.,  did  not  independently  change  during  the  test. 

In  other  words,  the  only  perturbations  were  exhaust  nozzle  area  changes 
and  varying  the  three  compressor  bleeds.  it  follows,  therefore,  that 
only  the  fault  coefficients  with  respect  to  these  four  disturbances  can 
be  computed  empirically. 

Once  the  numerical  values  for  the  fault  coefficients  have  been  deter- 
mined, applying  them  to  engine  diagnostics  is  relatively  straightforward. 

In  general,  if  n faults  are  to  be  computed,  (n  + 1)  measurements  are 
required.  One  measurement  serves  as  the  abscissa  against  which  n deviations 
from  baseline  values  are  obtained.  If  the  column  vector  x represents  the 
set  of  n sought  after  deltas  (such  as  changes  in  component  efficiencies, 
bleed  flows,  geometries,  etc.)  and  the  column  vector  mp  is  the  set  of  n 
delta  measurements,  then  they  are  related  by  an  nxn  ra^Trlx  of  fault 
coefficients  FI  as: 

Eq.  1 mj^  = Fj_x 

A subscript  1 is  attached  to  m and  F to  emphasize  that,  in  general, 
there  may  be  more  than  one  set  of  sensor  measurements  which  can  be  used 
to  diagnose  a given  set  of  faults.  If  the  measurement  deltas  bear  a 
meaningful  relationship  to  the  sought  after  deltas,  then  Fj  in  Eq.  1 
will  have  an  Inverse  and  we  can  solve  Eq.  1 for  to  obtain: 

Eq.  2 _x  = (Fj^  ~ 1)  Hi  = Hi 

The  evaluation  of  Parameter  Interrelationships  approach  using  the 
Phase  II  data  shows  that  the  approach  is  capable  of  identifying  the 
presence  of  faults,  singly  and  multiply.  The  computed  fault  levels 
generally  agree  quite  well  with  the  known  implanted  fault  level.  Com- 
plete results  appear  in  Appendix  B. 
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3.  PiaRnostlcs  on  Variable  Geometry  Engines 

The  need  for  the  capability  to  simulate  the  performance  of  a gas 
turbine  engine  with  variable  geometry  arose  from  the  consideration  of 
gas  path  analysis  as  one  element  of  a diagnostic  package  for  a 1980-85 
aircraft.  During  the  period  from  January  1974  through  March  1975,  a 
variable  geometry  engine  cycle  deck  was  developed  at  NAPTC.  The  specific 
gas  turbine  cycle  simulated  by  this  deck  is  a two-spool,  T-fan2,  mixed- 
flow  turbofan  with  afterburning.  In  addition,  this  cycle  deck  allows  for 
variable  area  high  and  low  pressure  turbine  (HPT  and  LPT)  nozzles  and 
variable  exhaust  nozzle  (A8)  area. 

Basically,  this  cycle  deck  consists  of  the  thermodynamic  equations 
of  a two-spool,  mixed-flow  turbofan.  These  equations  account  for  flow 
continuity,  energy  balance,  speed  matching,  and  pressure  balance.  In 

addition,  the  program  is  required  to  satisfy  the  requested  HPT,  LPT,  and  | 

A-8  areas.  The  thermodynamic  properties  of  air  are  obtained  from  GASTAB, 
which  is  a computerized  gas  table. 

This  cycle  deck  has  been  written  as  a design  or  Parametric  deck,  so  j 

that  any  engine  in  the  two-spool,  T-fan,  mixed-flow  turbofan  catagory  I 

can  be  simulated.  To  simulate  a particular  engine,  the  engine  design 

point  information  must  be  provided  to  the  program,  operating  in  the  ! 

Design  Mode.  j 

i 

Operating  in  its  Normal  Mode,  the  program  accepts  inputs,  such  as  i 

altitude  and  Mach  number,  horsepower  and  bleed  air  extraction,  changes  in  j 

component  performance,  changes  in  HPT,  LPT,  and  exhaust  nozzle  areas,  and  : 

rating  code  (ZRC) , which  is  the  power  level  input.  The  program  can  be  i 

operated  to  generate  a single  data  point  or  to  generate  a "power  series"  ’ 

at  the  given  input  conditions.  J 

( 

! 

The  data  generated  by  this  cycle  deck  is  output  to  the  line  printer  j 

and/or  magnetic  tape.  Several  output  formats  are  available  to  provide  1 

flexibility. 

Some  of  the  features  of  this  program  which  make  it  flexible  are  as  i 

follows.  Each  compressor  and  turbine  is  individually  represented  by  per-  ; 

formance  maps.  The  compressor  maps  are;  (a)  referred  speed  as  a function 
of  pressure  ratio  and  referred  flow  and,  (d)  efficiency  as  a function  of 
pressure  ratio  and  referred  flow.  The  turbine  maps  are:  (a)  flow  para- 

meter as  a function  of  pressure  ratio  and  referred  speed  and,  (b)  effi- 
ciency as  a function  of  velocity  ratio  and  corrected  work.  These  are 

stored  as  percentages  of  design  so  that  the  same  maps  can  be  used  for  ^ 

many  engines.  Pressure  drops  are  not  considered  to  be  constant.  The  fan  I 

duct,  burner,  diffuser,  mixer,  and  afterburner  friction  pressure  drops  are  4 


^T-fan  implies  fan  and  low  pressure  compressor  on  low  spool 
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computed  as  a function  of  the  referred  flow  and  the  design  value  of  the 
pressure  drop,  providing  more  realistic  calculations.  Also,  the  after- 
burner momentum  pressure  loss  is  computed  as  a function  of  afterburner 
Mach  number  and  temperature  rise.  Wherever  possible,  the  program  coding 
was  organized  into  building  blocks,  that  is  a subroutine  for  a compressor, 
a burner,  a turbine  vane,  a turbine  blade,  etc.  This  has  been  done  to 
facilitate  construe  tion  of  other  engine  cycles  without  duplicate  effort. 
Another  option  within  this  program  involves  the  number  of  turbine  stages. 
Both  the  high  and  low  pressure  turbines  tnve  the  option  of  being  single 
or  two-stage  components  independently. 

Once  data  has  been  written  cn  ro’';netic  tape,  there  are  a number  of 
ways  to  retrieve  desired  infonaatlon,  i)ne  option  is  to  list  all  runs  on 
tape.  It  is  also  possible  to  output  a chosen  run  or  set  of  runs  in  any 
of  the  data  output  formats. 

Another  tape  feature  which  is  available  is  called  Tape  Edit.  This 
allows  output  of  up  to  eight  stored  values  (inputs  and/or  outputs)  for 
all  or  any  selected  group  of  runs.  The  "editing"  is  accomplished  by 
disqualifying  runs  which  do  not  satisfy  the  editing  criteria,  and  any 
stored  value  may  be  used  as  an  editing  criteria.  Each  editing  criteria 
must  be  assigned  a minium  and  maximum  value. 

The  building  block  approach  to  the  coding  of  this  deck  allows  for 
the  creation  of  other  gas  turbine  cycles  with  a minimum  effort.  All 
of  the  programming  required  for  input/output  and  program  control  will 
remain  essentially  unchanged.  The  peripheral  capabilities  are  independent 
of  the  cycle  and  are  therefore,  available  for  use  on  any  cycle. 

The  cycles  which  are  planned  to  be  developed  in  the  immediate  future 
are  the  (a)  two-spool,  unmixed  turbofan  with  duct  burning  and/or  after- 
burning, (b)  two-spool,  mixed  turbofan  with  afterburning,  (c)  two-spool, 
regenerative  turboshaft  with  a free  power  turbine  and,  (d)  high  bypass 
ratio  turbofan. 

0 ' 

Variable  geometry  data  from  this  cycle  deck  was  supplied  to  Hamilton 
Standard  Division,  United  Technologies  Corporation  for  study  using  the 
parameter  interrelationships  method  of  gas  path  analysis. 

This  data  represents  an  engine  whose  geometry  schedules  are  designed 
for  9,144  meters  (30,000  ft.),  Mach  0.8  Cruise  condition  because  the 
cruis«.’  condition  is  an  appropriate  choice  for  steady  state  gas  path 
analysis.  The  geometry  schedules  were  selected  to  provide  Installed 
performance  gains.  This  Is  accomplished  by  maintaining  the  inlet  design 
airflow  while  reducing  thrust  v ^hereby  reducing  or  eliminating  spillage 
drag)  and  by  opening  the  exhaust  nozzle  to  reduce  boattail  losses  (aft 
end  losses).  Additionally,  maintaining  overall  compression  ratio  tends  to 
minimize  uninstalled  TSFC.  A complete  description  of  this  data  and  a 
more  thorough  discussion  of  the  cycle  deck  itself  appears  as  Appendix  C. 
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4.  TEDD  Sensitivity  Study 

In  developing  an  engine  diagnostic  system,  sensitivity  studies 
are  necessary  to  determine  which  engine  parameters  are  most  responsive 
to  changes  in  engine  condition  caused  by  various  faults;  and  therefore, 
should  be  used  in  the  diagnosis  of  these  faults.  In  addition,  the 
measurement  accuracy  of  these  parameters,  required  to  detect  given 
minimum  fault  levels,  must  be  determined.  Sensitivity  studies  can  take 
various  forms  and  the  results  expressed  in  different  ways.  They  are 
needed  to  establish  monitoring  parameter  requirements  for  engine  thermo- 
dynamic performance  analysis,  but  are  also  done  on  various  engine  sub- 
systems. Some  examples  of  sensitivity  studies  are  given  in  Reference  12 
which  describes  the  development  of  a methodology  to  identify  thermodynamic 
performance  monitoring  parameters  and  their  sensitivities,  and  also  gives 
the  results  of  two  computer  studies  run  on  engine  lube  system  models  to 
determine  the  relative  sensitivity  of  lube  system  parameters  to  failure/ 
deterioration  of  various  lube  system  components.  The  Hamilton  Standard 
approach  to  multiple  fault  diagnostics  also  provides  a direct  measure  of 
system  sensitivity  (Appendix  A of  Reference  11). 

When  the  TEDD  Phase  II  steady-stage  performance  analysis  program  was 
developed,  the  relative  sensitivity  of  each  of  23  elements  (pairs  of 
corrected  performance  parameters)  was  identified,  and  this  information  was 
used  in  setting  up  the  Directional  Analysis  and  Quantitative  Analysis 
(calculated  fault  level  and  Confidence  Factor)  parts  of  the  program 
(Reference  2).  During  TEDD  Phase  III,  a system  sensitivity  study  was 
conducted  on  the  TEDD  program  to  determine  the  effects  of  parameter  errors 
on  the  diagnostic  capabilities  of  the  system  and  to  get  an  indication  of 
the  measurement  accuracies  required.  The  results  of  this  study  are 
presented  in  Appendix  D. 
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Vibration  Experiments  Listing  and  Results 

Note:  All  readings  under  50  percent  are  considered  valid. 

Experiment  Class:  Gear/Pump /Accessory  Drive  - Local  Defects 

Function  SW  - Local 


Experiment 

Sensor 

Reading  Percent 
Idle  Intermediate 

No. 

Component 

No. 

Eng 

Tape 

Eng 

Tape 

1 

Tower  Shaft  Gears  (9),  (20) 

1 

13 

8 

10 

20 

2 

Gears  (21',,  (22) 

1 

16 

6 

18 

20 

3 

Gears  (27A),  (27B),  (35) 

1 

17 

6 

20 

20 

4 

UHP  Drive  Gear  (37) 

1 

21 

8 

22 

20 

5 

I ' Drive  Gear  (23) 

2 

20 

10 

22 

20 

6 

osj  35  Tooth 

2 

20 

10 

26 

18 

/ 

Starter  Drive  (28) 

1 

15 

8 

21 

20 

8 

uears  (24A) , (24B),  (24C) 

2 

21 

8 

17 

25 

9 

De-air  Drive  Gear  (25) 

2 

21 

10 

19 

24 

10 

Fuel  Pump  Drive  (29) 

1 

19 

10 

19 

20 

11 

Fuel  Pump  (41) 

1 

21 

6 

20 

20 

12 

Fuel  Pump  (42) 

1 

16 

8 

25 

20 

13 

Gears  ( 30A) , (30B) 

1 

17 

8 

16 

20 

14 

Main  Oil  Pump  Drive  (31) 

1 

22 

6 

23 

20 

15 

N2  Tach  Drive  (32) 

1 

21 

6 

20 

20 

16 

17  - 24 

Fuel  Boost  Drive  (26) 
Unassigned 

2 

19 

10 

19 

19 

Experiment  Class:  Cear/Pump/Support  Bearing  - Overall  Defects 

Function  SW  - Gross 


25 

Shaft 

A - C (9)  and  SB 

1 

99 

22 

99 

24 

26 

Shaft 

A - G (20)  and  SB 

1 

52 

20 

77 

24 

27 

Sha  f t 

B - G (21)  and  SB 

1 

34 

20 

90 

24 

28 

Shaft 

B - G (22)  and  SB 

1 

31 

40 

33 

65 

29 

Shat  t 

C - G (23)  and  SB 

2 

99 

40 

71 

80 

30 

Snaft 

C - CSD  2''  and  SB 

2 

77 

99 

44 

80 

31 

Shaf  t 

W - CSD  35 

2 

59 

99 

99 

70 

32 

Shatt 

D - G (24A)  and  SB 

2 

99 

45 

99 

95 

33 

Shaft 

D - G (24B) 

2 

99 

60 

65 

99 

34 

Shaft 

D - G (24C)  and  SB 

2 

99 

99 

99 

99 

35 

Shaft 

E - G (25)  and  SB 

2 

99 

20 

99 

95 

36 

Shaft 

G - G (27A)  and  SB 

1 

14 

40 

31 

80 

37 

Shaft 

H - G (27B) 

1 

28 

65 

23 

99 

38 

Shaft 

H - G (35)  and  SB 

1 

65 

95 

85 

99 

39 

Shaft 

J - G (37)  and  SB 

1 

99 

90 

99 

99 

40 

Shaft 

J,  K - UHP,  LPH 

1 

98 

99 

99 

99 

41 

Shaft 

L - G (28)  and  SB 

1 

22 

30 

50 

80 

42 

Shaft 

M - G (29)  and  SB 

1 

77 

40 

99 

99 

43 

Shaft 

M - PPG  (40) 

1 

89 

90 

84 

99 

44 

Shaft 

Ml  - FPG  (40) 

1 

99 

99 

90 

99 

45 

Shaft 

P - FPB  (41) 

1 

99 

99 

99 

99 

46 

Shaft 

Q - FPG  (42) 

1 

99 

99 

84 

99 

47 

Shaft 

R - G (30A) 

1 

32 

40 

19 

99 

35 
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FIGURE  3:  VIBRATING  CYLINDER  PRESSURE  TRANSDUCER  SCHEMATIC 
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FIGURE  5: 


K WEST  DEBRIS  MONITOR  EQUIVALENT  CIRCUIT 
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FIGURE  8:  OVERALL  GEAR  DEFECT  ANALYSIS 
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FIGURE  9:  DIGITAL  COMB  FILTER 
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APPENDIX  A j 

TEDD  Data  Smoothing  Study  j 

A . Introdn ction  I 

The  purpose  of  this  study  was  to  analyze  data  smoothing  techniques 
in  order  to  improve  on  the  smoothing  calculations  used  in  TEDD  Phase  II 
software,  particularly  in  the  logic  of  the  mode-detect  section  of  the 
program. 

During  Phase  II  a limited  moving  average  technique  was  employed. 

It  was  limited  because  of  core  size  restrictions,  and  was  only  partially 
successful  for  the  following  reasons:  (1)  lack  of  sufficient  data  to 

define  high  rate  of  change  engine  events,  and  (2)  inflexibility  to 
change  the  number  cf  points  included  in  the  moving  average  during  real 
time  processing. 

Specifically,  one  severe  problem  was  encountered  in  the  mode  de- 
tection logic.  Due  to  noise  fluctuations,  the  speed  derivatives  taken 
with  respect  to  time  were  unstable.  To  have  good  mode  detection,  it 
was  therefore  necessary  to  use  derivatives  that  were  taken  over  a longer 
time  span.  However,  because  of  these  longer  derivative  response  times, 
the  mode  detection  logic  could  not  properly  recognize  rapid  transients, 
such  as  Bodies,  that  have  two  rapid  changes  in  slope  polarity. 

The  data  used  for  this  analysis  was  stored  on  magnetic  tape  during 
TEDD  Phase  II  TF30-P-408  engine  testing,  TF30-P-412  data  from  a NAPTC 
Fleet  Service  Test  Program  was  also  used. 

This  report  describes  the  data  used  and  the  different  techniques 
used  to  smooth  it.  Criteria  are  developed  that  judge  these  techniques, 
and  some  viable  methods  of  smoothing  applicable  to  diagnostic  data  are 
presented. 

B.  TRANSIENT  DATA  - GENERAL  DISCUSSION 

TEDD  transient  data,  especially  that  used  for  mode  detection,  have 
characteristics  that  make  any  fixed  degree  of  smoothing  difficult  to 
apply.  During  Phase  II  of  TEDD,  short  time  duration  Bodies  (rapid 
deceleration  followed  immediately  by  a rapid  acceleration)  were  encount- 
ered which  occurred  in  less  than  three  seconds.  The  high  rates  of 
change  in  the  speed  parameters,  combined  with  a sample  rate  of  about  5 
samples  per  second,  eliminated  the  possibility  of  smoothing  these 
parameters,  because  any  attempt  to  smooth  the  noise  would  distort  the 
rates  of  change  so  vital  for  transient  mode  detection.  The  other 
extremes  encountered  were  slow  transients  and  steady  state  conditions 
where  there  were  very  little  changes  in  parameter  magnitudes.  In  such 
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cases,  data  was  sampled  about  every  0.5  to  2.5  seconds  and  a high  degree 
of  smoothing  could  be  done  without  any  concern  about  distorting  the 
rates  of  change.  However,  other  criteria  such  as  accuracy  entered  into 
consideration. 

In  this  report,  moving  averages  and  exponential  smoothing  techniques 
are  applied  to  the  data.  Due  to  Che  uniqueness  of  the  application  of 
TEDD  data,  criteria  to  judge  the  results  of  data  smoothing  must  be 
developed. 

C.  CRITERIA 

Storage  Required 

This  is  concerned  with  how  much  data  storage  is  required  to 
perform  Che  various  smoothing  techniques  and  also  how  much  storage  the 
technique  itself  will  occupy.  Obviously,  the  less  storage  required  the 
better. 

2.  Ability  to  Adjust  Rates  of  Response 

The  nature  of  TEDD  data  is  such  that  it  is  not  known  whether  the 
next  data  sample  will  have  a similar  magnitude  as  the  last,  or  repre- 
sent a radical  change.  This  characteristic  dictates  that  the  techniques 
must  be  able  to  assign  weights  to  the  data  depending  on  certain  engine 
parameter  rates  of  change.  How  well  the  techniques  can  do  this  is  an 
important  criteria. 

3.  Simplicity  of  Computation 

Since  the  ultimate  application  of  TEDD  programming  is  a real  time 
situation  in  which  the  computations  are  repeated  with  each  new  data  scan, 
and  since  applied  smoothing  will  be  essentially  an  iterative  process, 
the  amount  of  time  spent  smoothing  is  critical.  Those  techniques  that 
employ  fewer  mathematical  steps  and  can  work  effectively  on  less  data 
will  be  of  high  interest.  Not  only  will  less  memory  be  used,  but  they 
will  also  be  faster. 

A.  Response 

This  criteria  applies  mainly  to  transient  conditions.  Due  to 
mode  detection  techniques,  little  distortion  of  rates  of  change  of 
certain  key  parameters  can  be  tolerated.  The  tradeoff  here  is  to  get 
the  true  rate  of  change  of  the  parameter  by  smoothing  out  noise  but 
not  to  the  extent  of  distorting  the  true  event. 

5.  Accuracy  and  Noise  Reduction 

This  criteria  is  used  to  compare  how  well  the  smoothed  data 
compares  to  the  model  of  the  true  event.  Standard  deviation  defined 
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as  the  measure  of  how  much  the  data  varies  from  the  mean  of  the 
sample,  will  be  used  to  judge  this  criteria. 

D . TECH MQUES  APPLIED 

L . Moving  Average  Techniques 

This  Is  me  of  the  methods  chosen  for  evaluation  since  it  was 
used  previously  daring  TEDD  Phase  II  and,  although  it  had  some  draw- 
backs, was  fairly  successful  and  accurate  enough  to  warrant  further 
investigation.  It  minimizes  the  sum  of  the  squares  of  the  differences 
between  the  most  recent  N observations  and  the  best  guess  for  that 
a'/erage.  The  simplicity  of  computation  for  this  method  is  an  obvious 
plus  in  a TEDD  operating  environment. 

The  moving  average  technique  operates  on  a time  series  by  computing 
a series  of  averages,  based  on  a given  number  of  points  in  the  original 
series.  For  example,  if  Xi  is  the  1th  data  point  in  a time  series  of 
data  points,  and  N is  the  number  of  discrete  points  to  be  averaged,  the 
corresponding  point  in  the  average  series  would  be  computed  as  shown 
below: 


N 


In  the  same  way  the  following  point  X^+j^  would  be: 

j = 1+1 

j = i-N+2 


N 

The  value  that  is  picked  for  N,  the  number  of  data  points  to  be  included 
in  each  average,  will  not  be  the  same  for  all  data.  A small  value  of 
N is  useful  for  eliminating  small  variations  in  the  data,  but  is 
relatively  ineffective  for  large  variations.  Conversely,  if  too  large 
a value  of  N is  selected,  there  is  the  problem  of  excessive  smoothing, 
thereby  ignoring  some  of  the  true  event. 
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An  inherent  problem  in  applying  this  method  is  the  lack  of  data 
for  Initial  conditions  since  there  are  no  previous  data  points  to  use 
in  computing  the  initia]  average.  Extra  weight  will  be  given  to  the 
initial  data  points  until  N readings  are  obtained. 

2 . Exponential  Averaging  Technique 

Exponential  smoothing  is  a method  of  estimating  each  new  data 
point  in  a time  series  by  letting  the  new  smoothed  value  be  equal  to 
a fraction  of  the  previous  smoothed  value,  plus  a fraction  of  the  new 
data  point.  This  can  be  written  as: 

St(x)  = aXt  + (1-a)  St_2^(x) 

St(x)  = new  smoothed  data  point 

a = smoothing  constant;  the  weight 
given  to  the  new  data  point 

St_j^(x)  = previous  smoothed  data  point 
X(-  = new  data  point 

This  method  of  smoothing  is  simple  to  compute  and  has  the  advantage  of 
small  amounts  of  data  storage.  Only  two  words  of  storage  are  required 
for  each  parameter  instead  of  the  N words  required  by  the  moving 
average  technique.  Another  advantage  is  its  flexibility  because  only 
a needs  to  be  varied,  to  change  the  rate  of  response.  When  a is  small, 
the  method  gives  litt'le  weight  to  the  new  data,  which  is  good  for 
smoothing  random  noise  in  a steady  stage  condition,.  A large  a will 
weight  the  new  data  heavily  so  it  responds  to  rapid  changes  in  the  time 
series,  which  is  necessary  for  tracking  transients. 

To  illustrate  the  weighting  of  data  in  both  the  moving  average  and 
exponential  smoothing  techniques,  see  Figure  A-1  where  N = 4 for  the 
moving  average  and  a = 0.6.  In  the  exponential  smoothing  technique, 
the  weight  of  data  for  a given  age  may  be  found  by  substituting  for  the 
previous  smoothed  value  the  equation  that  defines  in  the  same  terms  a 
still  earlier  smoothed  value  (Reference  1).  Since  a = 0.6,  the  current 
data  point  has  a weight  of  0.6  and  the  previous  data  points  are  weighted 
0.24,  0.096,  0.0384,  respectively. 

St(x)  = aXt  + (l-a)St-l(x)  substituting  for  the  St-l(x)  term 
St(x)  = aXt  + (1-a)  [aXt_2^  (l“0i)St;_2(x)  ] and  continuing  to  substitute 

k = t-1 

for  N earlier  values  of  St(x)  St(x)  = a ^ (l-a)*'^Xj._j^  + (l-a)tXQ 

k = 0 
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The  weight  given  to  the  X^-k  term  with  a = 0.6;  for  K = 1,  2,  3 are:  ^ 

.6(1-. 6)1  = .6(.4)  = 0.24  i 

i 

.6a-. 6)2  = ,6(.16)  = 0.096  j 

.6(1-. 6)3  = .6(0.064)  =*  0.0384 

When  appj-ing  the  smoothing  techniques  to  transient  data  it  was 
obvious  chat  no  fixeo  degree  of  smoothing  would  satisfy  the  criteria 
already  discussed.  As  an  example  of  this  in  Figure  A-2,  a moving 
average  with  N = 2 is  shown  along  with  an  exponentially  smoothed  timed 
series  with  a = 0.2.  it  is  seen  that  the  moving  average  follows  the 

real  data  with  almost  no  change  in  amplitude,  but  with  a very  slight  ; 

time  lag.  The  exponentially  smoothed  trace,  with  little  weight  given 
to  each  new  reading,  smoothes  the  steady  state  interval  (from  5.0  to  j 

6.8  seconds)  rather  well.  During  the  transient,  however,  there  is  ' 

serious  distortion. 

If  the  other  extremes  (N  = 8,  a = 0.8)  are  plotted  for  the  same  j 

data,  it  is  seen  that  the  smoothing  results  are  almost  reversed. 

Figure  A-3  shows  that  the  exponential  method  tracks  the  transient  almost 

exactly  and  does  not  have  as  much  time  lag  as  did  the  moving  average  in  ; 

Figure  A-2.  Also,  as  expected,  the  moving  average  smoothed  the  steady 

state  data  well  but  distorted  the  transient.  S 

Figures  A-2  and  A-3  illustrate  the  '’classical"  problem  of  data 
smoothing  for  TEDD  applications.  The  main  purpose  of  the  TEDD  computer  ' 

program  is  to  diagnose  engine  problems.  Therefore,  it  is  desirable  ] 

to  spend  only  a small  percentage  of  time  reading  and  smoothing  data  j 

and  to  devote  a larger  percentage  to  doing  actual  diagnostics.  The  ij 

conflict  arises  when  trying  to  get  good  smoothing  results.  A high  j 

sample  rate  gives  better  smoothing,  but  takes  more  time  to  process.  |i 

A low  sample  rate  takes  less  time  to  process,  but  gives  poorer  results  ) 

in  smoothing.  Figure  A-4  states  this  problem  graphically.  i 

E.  TECHNIQUES  AN/ALYZED  | 

When  analyzing  the  data,  it  was  convenient  Co  consider  the  steady  : 

state  and  transient  data  separately.  This  is  a natural  division  in  I 

that  it  gives  the  two  extremes  of  data  characteristics  to  be  encountered.  I 

I 

The  steady  state  case  wJll  be  considered  first.  Accuracy  and  i 

how  sample  rate  affects  accuracy  are  the  main  criteria.  Standard 

deviation  will  be  used  as  the  Indicator  of  accuracy.  The  standard  j; 

deviation  is  defined  as  the  measure  of  how  much  the  data  varies  from  I 

the  mean  of  a given  sample.  The  expression  is  given  below.  I 
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1=1 


N-1 

Where  Xj.  is  the  ith  data  point  in  the  sample,  X is  the  mean  of  the  sample. 

Three  parameters  (PLA,  N2,  and  PB)  will  be  analyzed.  PLA  and  N2 
are  key  parameters  for  mode  detection  and  PB  is  representative  of  a 
pressure  whose  signal  noise  has  to  be  reckoned  with.  The  data  was 
processed  at  sample  rates  of  100,  50,  10  samples/second.  Moving 
averages  with  N ranging  from  2 to  10  were  computed.  Exponential 
smoothing  was  done  for  a = 0.9  to  a = 0.1.  The  ranges  for  N and  a were 
selected  to  show  trends  and  do  not  necessarily  show  the  optimum 
smoothing  that  is  possible.  The  results  of  the  processing  are  shown  in 
Tables  I,  II,  and  III. 

The  tables  show  that  each  parameter  responds  to  degrees  of  smoothing 
in  different  ways.  Not  only  does  the  magnitude  of  the  a's  vary  from 
parameter  to  parameter,  but  for  increasing  values  of  N the  standard 
deviation  does  not  always  improve.  Because  of  various  frequency  com- 
ponents present  in  the  different  parameters,  higher  sampling  rates  did 
not  always  give  lower  standard  deviations.  However,  the  general  trends 
shown  in  the  tables  indicate  that  as  N went  from  2 to  10  and  a from  1.0 
to  0.9,  the  standard  deviation  of  the  parameters  improved;  likewise  for 
increasing  sampling  rates.  It  is  obvious  from  the  tables  that  for  any 
degree  of  smoothing,  the  smoothed  data  has  a lower  standard  deviation 
than  the  unsmoothed  data.  For  the  values  of  N = 10  and  a = 0.1,  the 
exponential  technique  was  consistently  better. 

Figure  A-5  is  a graphic  illustration  of  the  effectiveness  of 
exponential  smoothing.  The  unsmoothed  data  has  a a of  38.665  which  is 
interpreted  as  the  signal  varying  ±38.665  rpm  about  the  mean  of  the 
sample.  The  smoothed  trace  has  an  a of  6.734.  Overall,  the  percentage 
of  relative  error  improves  from  0.591%  to  0.102%. 

AQ  X 100  = % relative  error 

Q 

77.33  X 100  = 0.591% 

13078 

13.468  X 100  = 0.102% 

13078 
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TABLE  I 


EFFECT 

OF  VARYING  N. 

a,  AND 

SAMPLE  RATE  ON 

STANDARD 

DEVIATION  OF 

STEADY 

STAGE  PLA  DATA 

PLA, 

DEGREES 

N 

100 

samples/sec 

a 

50 

samples/sec 

0 

10  samples/i 
0 

2 

.0737 

.0436 

.0602 

3 

.0532 

.0311 

.0500 

4 

.0350 

.0308 

.0413 

5 

.0245 

.0233 

.0323 

6 

.0260 

.0278 

.0253 

7 

.0284 

.0269 

.0199 

8 

.0274 

.0252 

.0269 

9 

.0241 

.0285 

.0298 

10 

,0236 

.0280 

.0359 

Unsmoothed 

data  a 

.0947 

.0918 

.0779 

a 

.1 

.0138 

.0109 

.0098 

.2 

.0226 

.0178 

.0197 

.3 

.0316 

.0244 

.0293 

.4 

.0140 

.0314 

.0383 

.5 

.0507 

.0389 

.0465 

.6 

.0604 

.0472 

.0539 

.7 

.0698 

.0566 

.0656 

.8 

.0788 

.0669 

.0667 

L 
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TABLE  II 


EFFECT 

OF  VARYING  N,  a 

, AND  SAMPLE  RATE  ON 

STANDARD 

DEVIATION  OF 

STEADY  STATE  N2  DATA 

N2,  RPM 

100 

N 

samples/ sec 
a 

50  samples/sec 
a 

10  samples/sec 
a 

2 

20.666 

25.903 

21.823 

3 

14.66 

19.959 

17.957 

4 

12.881 

17.667 

15.454 

5 

11.643 

16.068 

19.202 

6 

11.332 

16.370 

19.141 

7 

13.053 

16.111 

22.219 

8 

12.888 

16.103 

22.373 

9 

13.301 

16.302 

24.927 

10 

13.946 

16.189 

24.777 

Unsmoothed  data  o 

38.665 

36.927 

30.715 

a 

.1 

6.734 

10.479 

9.5036 

.2 

7.605 

12.223 

11.469 

.3 

10.156 

14.687 

13.358 

.4 

13.261 

17.412 

15.573 

.5 

16.682 

20.304 

17.894 

.6 

20.380 

23.336 

20.264 

.7 

24.372 

26.502 

22.677 

.8 

28.702 

29.812 

25.177 

.9 

33.437 

33.280 

27.833 

( 

I 

) 


I 


i 


1 


\ 


A-8 


NAPTC-PE-88 


TABLE  III 


EFFECT  OF  VAPYING  N, 

g,  AND 

SAMPLE  RATE  ON 

STANDARD 

DEVIATION  OF 

STEADY 

STATE  Pb  DATA 

Pb,  Psia 

N 

100  samples/sec 
0 

50 

samples/sec 

0 

10  samples/sec 
0 

2 

.6322 

.3610 

.5392 

3 

.4653 

.2582 

.4926 

4 

.3127 

.2705 

.4530 

5 

.2239 

.2224 

.3940 

6 

.2297 

.2242 

.3436 

7 

.2454 

.2174 

.3358 

8 

.2432 

.2179 

.3641 

9 

.2240 

.2537 

.3976 

10 

.2034 

.2593 

.4116 

Unsmoothed  data 

0 .7805 

.7656 

.6649 

g 

.1 

.1352 

.1059 

.1240 

.2 

.2029 

.1474 

.1997 

.3 

.2770 

.1989 

.2800 

.4 

.3551 

.2561 

.3562 

.5 

.4348 

.3192 

.4236 

.6 

.5137 

.3892 

.4820 

.7 

.5893 

.4673 

.5332 

.8 

.6597 

.5550 

.5793 

.9 

.7235 

.6538 

.6225 
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In  Figure  A-5,  there  is  a good  Illustration  of  response.  The 
first  data  point  in  the  time  series  is  13,400  rpm.  With  a = 0.1  there 
is  little  weight  given  to  each  new  data  point.  The  lag  of  the  smoothed 
data  is  slightly  less  than  0.4  seconds.  This  lag  would  have  to  be 
taken  into  account  when  this  technique  is  applied  to  data.  A higher 
alpha  could  be  applied  until  the  initial  lag  time  is  reached. 

The  transient  data  smoothing  analysis  Illustrates  the  criteria  of 
response  and  noise  reduction.  Again,  as  in  the  steady  state  case,  PLA, 

N2  and  PE  are  the  parameters  to  be  examined.  One  criterion  for  trans- 
ient data  will  be  time  lag,  defined  as  the  time  in  seconds,  the  smoothed 
values  lag  the  actual  data.  Also  how  smooth  the  speed  and  PLA  time 
derivatives  are  will  determine  how  well  the  mode  detection  routine  will 
function.  Effects  of  sample  rates  will  not  be  analyzed  because  the 
minimum  sample  rate  will  be  dictated  by  the  worst  case  parameter.  In 
this  situation  PB  was  the  worst  case  parameter,  having  a negative  slope 
of  about  56  Psla  in  0.08  seconds.  It  is  considered  that  a minimum  of 
three  equally  spaced  data  points  is  required  to  define  a given  slope. 

The  transient  event  analyzed  was  a Bodie  combined  with  a stall. 

Tables  IV  through  VI  show  time  lags  for  moving  averages  with  N ranging 
from  2 to  10  and  exponential  smoothing  with  a's  ranging  from  0.1  to  0.9. 

The  time  lags  recorded  in  Tables  IV  through  VI  occur  at  the  most  abrupt 
change  in  slope  during  the  transient.  Again,  a compromise  must  be  made 
between  response  and  noise  reduction. 

The  steady  state  analysis  showed  that  for  TEDD  application,  the 
exponential  smoothing  did  a better  job  of  noise  reduction.  For  this 

reason,  the  moving  average  technique  will  not  be  considered  further.  I 

The  mode  detection  routine  used  in  TEDD  Phase  II  looked  for  a PLA  ; 

change  of  one  degree  in  one  second  to  help  determine  transient  and 
steady  state  modes.  The  PLA  signal  was  repeatable  to  within  ±0.5 
degrees.  To  apply  smoothing,  it  would  be  desired  to  smooth  the  data 
to  be  repeatable  to  within  ±0.1  degrees. 

It  is  estimated  that  a lag  of  0.25  second  from  the  true  event  to 
its  detection,  for  current  engines,  is  tolerable.  From  Table  II  it  1 

is  determined  that  an  a of  0.5  will  meet  the  response  criteria  and  in  j 

Table  I an  a of  0.5  has  a standard  deviation  of  0.0465  which  is  1 

acceptable  for  repeatability.  Figure  A-6  shows  the  results  of  j 

smoothing  with  a = 0.5  and  Illustrates  that  the  above  criteria  can  be  3 

met. 

PB,  the  parameter  that  dictates  sample  rate  for  transient  analysis, 
has  the  fastest  rate  of  change  of  all  the  parameters  examined  for 
diagnostics.  In  this  case,  response  is  the  most  important  criterion 
and  any  noise  reduction  we  can  get  will  be  acceptable.  As  mentioned 
before,  during  a stall  PB  dropped  as  fast  as  56  Pisa  in  0.08  seconds. 

To  be  able  to  detect  this  slope  before  the  event  is  over,  a time  lag 
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TABLE  IV 


EFFECT  OF  VARYING  N AND 

OL  ON  TIME  LAG 

OF  TRANSIENT  PLA 

DATA 

i 

PLA 

Samples/sec  = 

50 

N 

Time  Lab  (sec) 

a 

Time  Lay 

10 

.16 

.9 

.02 

9 

.14 

.8 

.02 

8 

.12 

.7 

.025 

7 

.10 

.6 

.04 

6 

.07 

.5 

.06 

5 

.065 

.4 

.10 

4 

.05 

.3 

.18 

3 

.03 

.2 

.31 

2 

.01 

.1 

.82 
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TABLE  VI 


EFFECT  OF  VARYING  N AND 

a ON  TIME  LAG 

OF  TRANSIENT  Pb  DATA 

pb 

Samples/sec  = 

50 

N 

Time  Lag  (sec) 

g 

Time  Lag 

10 

.08 

.9 

.01 

9 

.06 

.8 

.01 

8 

.06 

.7 

.01 

7 

.05 

.6 

.02 

6 

.04 

.5 

.03 

5 

.03 

.4 

.05 

4 

.03 

.3 

.05 

3 

.02 

.2 

.05 

2 

.02 

.1 

.06 
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of  0.02  seconds  would  be  acceptable.  Table  VI  shows  that  an  a of  0.6 
will  give  this  response.  In  Figure  A-7,  the  raw  data  is  plotted  with 
the  smoothed  data  and  illustrates  that  the  above  criterion  can  be  met. 
With  an  a of  0.6,  standard  deviation  of  steady  state  before  the  trans- 
ient is  0.39. 

Another  important  parameter  for  mode  detection  is  N2.  It  was  used 
in  the  mode  detection  logic  in  Phase  II  to  help  determine  steady  state 
or  transient  modes,  and  then  further  to  indicate  fast  or  slow  transient. 
As  stated  earlier,  a change  in  N2  over  a three  second  period  was  sought. 
If  the  absolute  value  of  the  delta  was  greater  than  160  rpm,  a trans- 
ient mode  was  indicated.  Please  note  this  was  not  the  only  criterion 
used  to  indicate  transient  or  steady  state  modes.  The  three  second 
interval  was  used  because  of  the  noise  in  the  N2  trace.  It  was  deter- 
mined that  for  the  given  signal  noise  level  and  sample  rate,  the  above 
was  the  best  criteria.  If  the  sample  rate  were  faster  and  exponential 
smoothing  methods  have  been  available,  the  mode  detection  would  have 
been  much  improved,  especially  in  the  transient  section.  Again,  looking 
at  the  Bodie-Stall  sequence  in  Figure  A-8,  it  is  seen  that  the  rate  of 
change  of  N2  to  be  detected  is  approximately  -625  rpm/sec.  To  detect 
this  slope  and  still  remain  as  close  in  time  to  the  true  event  as 
practical,  an  a of  0.5  is  used.  The  time  lag  is  0.03  seconds  with  the 
results  shown  in  Figure  A-8.  From  Table  II,  the  standard  deviation  is 
20.3.  This  is  a rather  large  deviation  about  the  mean  and  probably 
would  not  give  good  slope  definition.  A standard  deviation  of  10.479 
has  very  good  noise  reduction. 

F.  CONCLUSION 

When  comparing  the  two  methods  of  smoothing  with  the  criteria 
developed  earlier,  it  is  obvious  the  exponential  smoothing  technique  is 
better  suited  to  TEDD  applications  than  the  moving  average.  Table  VII 
shows  the  comparison.  To  apply  the  exponential  smoothing  to  TEDD  data, 
it  is  necessary  to  have  pre-deflned  separate  a's  for  steady  state  and 
transients,  depending  on  response  and  accuracy  desired  for  each  para- 
meter. Upon  deciding  which  mode  is  current,  the  appropriate  a is 
applied  to  the  data.  In  Figure  A-10,  an  a of  0.1  is  applied  to  the 
steady  state  section  and  an  a of  0.3  is  applied  to  the  transient  section. 
Programming  details  such  as  Initial  conditions  cannot  be  overlooked  when 
switching  a during  a transient. 
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CRITERIA 


Storage  Required 


Ability  of  adjust 
rate  of  response 

Simplicity  of 
computation 


Response 


Accuracy  and 
nose  reduction 


TABLE  VII 
MOVING  AVERAGE 
N X (#  of  parameters) 

Programming  more 
complicated 

More  instructions 
to  implement 


O.K.  for  proper  N 


O.K.  for  proper  N 
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EXPONENTIAL 

SMOOTHING 


2 words  for  each 
parameter 

Easy  to  implement 


Fewer  instructions; 
therefore  faster 


Always  better  for 
tested  range 

Always  better  for 
tested  range 
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FIGURE  A-4:  DATA  SMOOTHING  TRADE-OFFS  AND  OPTIMIZATION 


Mlninmm 

A GOODNESS  OF 
SMOOTHING 

B TIME  AVAIL- 
ABLE 


SAMPLE  RATE 


CURVE  A:  For  a high  sample  rate  we  get  a high  goodness  of  smoothing 

and  definition  of  event.  Conversely,  for  a low  sample  rate, 
a poor  goodness  of  smoothing  and  lack  of  definition. 

CURVE  B:  For  a high  sample  rate  there  is  little  time  for  diagnostic 

analysis  and  for  a low  sample  rate  there  is  more  time 
available  for  diagnostics. 

CURVE  C:  Algebraic  combination  of  A and  B,  the  minimum  of  which  indi- 

cates the  optimum  combination  of  sample  rate,  diagnostic 
time  available,  and  goodness  of  smoothing. 
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APPENDIX  B 

Multiple  Fault  Gas  Path  Analysis 
Applied  to  TF30-P-408  EnRlne  Data 
(SUMMARY  OF  HAMILTON- STANDARD  REPORT  HSER  6587) 

A.  INTRODUCTION 

Under  Phase  II  of  the  Turbine  Engine  Diagnostics  Development  Program, 
steady  state  data  from  the  TF30-P-408  engine  was  obtained,  with  and  with- 
out deliberately  Introduced  faults.  This  data  was  supplied  to  the  : 

Hamilton- Standard  Division  of  United  Technologies  Corporation,  under  I 

Contract  Number  N00140-7A-C-0582  for  the  purpose  of  applying  the  Parameter 
Interrelationships  method  of  gas  path  analysis.  The  results  of  this 
study  provide  a means  for  evaluating  this  approach  to  multiple  fault  gas 
path  analysis.  This  document  serves  only  as  a summary  of  the  effort 
conducted  by  Hamilton-Standard.  A thorough  description  of  the  Parameter 
Interrelationships  method,  as  well  as  a complete  set  of  results,  are  con- 
tained in  Reference  B-1. 

B.  DESCRIPTION  OF  DATA  SUPPLIED 

The  data  supplied  to  Hamilton-Standard  for  this  study  are  of  three 
types;  steady  state  engine  data  with  and  without  known  faults,  engine 
base  line  calibration  data,  and  TF30  cycle  deck  data  at  sea  level,  Mach 
number  >0.  Each  of  the  engine  parameters  was  acquired  200  times  through 
the  transient  data  system  over  a period  of  six  seconds  and  statistically 
averaged  to  a single  value.  Fifty  channels  of  data  comprise  a "snapshot" 
and  97  "snapshots"  were  provided  for  analysis.  Fifty-three  of  these 
"snapshots"  are  data  acquired  at  sea  level  static  with  one  or  two  delib- 
erately introduced  faults.  The  remaining  44  "snapshots"  are  data  acquired 
at  sea  level  with  inlet  ram  conditions  simulating  Flight  Mach  numbers  from 
.25  to  .75,  with  and  without  Introduced  faults.  Five  difference  faults 
were  introduced:  customer  bleed  air  extraction  from  the  ninth  (BL3)  and 

sixteenth  (BL4)  stage,  surge  bleed  air  extraction  (BL3.5)  and  two  abnormal 
exhaust  nozzle  areas  (AJN+) (AJN-) . The  engine  power  settings  ranged  be- 
tween 1.6  and  2.2  in  engine  pressure  ratio  (PT7/PT2). 

C.  DESCRIPTION  OF  ANALYTICAL  TECHNIQUES 

Whatever  form  it  may  take,  any  engine  diagnostic  system  must  rely  on 
discernable  changes  in  observable  parameters  in  order  to  detect  physical 
faults.  Physical  faults  in  an  engine  are  in  general  some  combination  of 
problems  such  as  erosion,  corrosion,  foreign  object  damage  (F.O.D.),  worn 
seals,  etc.  Although  the  physical  faults  result  in  changes  in  obaervable 
parameters  such  as  rotational  speeds,  temperatures,  pressures,  etc.,  it  is 
more  fundamental  to  consider  the  physical  faults  as  causing  changes  in 
engine  component  performance  parameters  such  as  component  efficiencies, 
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turbine  and  exhaust  nozzle  areas,  compressor  pumping  capacities,  etc.,  and 
that  these  changes,  In  turn,  result  in  discernible  changes  in  the  measur- 
able engine  parameters.  If  proper  equations  can  be  established  relating 
the  dependent  engine  variables  (measurable  quantities)  to  the  generally 
nonmeasurable  independent  component  performance  variables,  then  an  inverted 
set  of  relationships  can  be  obtained  relating  the  independent  variables  to 
the  dependent  variables.  This  could  then  form  the  basis  of  an  engine 
diagnostic  system  which  would  have  the  capability  of  detecting  multiple 
faults  as  well  as  the  magnitude  ana  direcf'cn  of  the  fault. 

This  is  the  basic  technique  employed  by  Hamilton-Standard  in  their  en- 
gine diagnostic  system.  The  major  obstacle  in  employing  such  a diagnostic 
scheme  is  in  establishing  the  proper  relationships  that  exist  between  the 
dependent  and  independent  engine  variables.  Two  different  approaches  were 
used  to  establish  these  relationships.  Before  discussing  these  approaches, 
it  is  appropriate  to  define  two  te^'ras.  The  first  term  is  "delta"  (A).  The 
term  "delta"  means  the  percent  of  point  difference  in  a particular  variable 
relative  to  its  nominal  or  base  line  value.  For  example,  if  yj^  and  are 
the  measured  values  of  two  engine  parameters  y and  x,  then  Ay  at  x = 
constant  and  Ax  at  y = constant  would  be  calculated  as  illustrated  in 
Figure  B-1.  Note  that  by  the  above  definition,  all  A's  are  dimensionless 
quantities  whose  value  for  any  variable  does  not  depend  on  the  units  used 
to  measure  the  variable.  The  second  term  is  "fault  coefficient."  This  is 
just  the  percent  change  in  a dependent  variable  per  unit  percent  change  in 
an  independent  variable.  For  example,  if  in  Figure  B-1,  Ay  at  x = constant 
was  the  result  of  a change  in  exhaust  nozzle  area  of  an  amount  AAjj^  only , 
then  the  fault  coefficient  of  y with  respect  to  Ajjj  at  x = constant  would 
be  simply  Ay/AAjj^.  Note  that  since  the  deltas  are  dimensionless,  the  fault 
coefficient  is  also  dimensionless.  Two  methods  are  used  in  this  study  to 
calculate  the  required  fault  coefficients.  The  first  method  uses  an 
elaborate  gas  path  analysis  of  the  engine  developed  by  Hamilton-Standard 
and  is  primarily  theoretical  in  nature.  The  second  method  is  empirical 
and  uses  available  base  line  test  data  as  well  as  engine  test  data  ob- 
tained with  known  changes  made  to  the  base  line  configuration. 

As  shown  in  Reference  B-2,  a general  Influence  coefficient  matrix  can 
be  derived  for  any  particular  gas  turbine  cycle.  This  matrix  represents 
a set  of  linearized  equations  which  Interrelates  the  various  dependent 
and  Independent  engine  performance  parameters.  Each  of  the  coefficient 
is  derived  by  proper  manipulation  of  the  differential  forms  of  the  basic 
thermodynamic  equations  for  each  engine  type.  Factors  such  as  component 
efficiency  variations;  specific  shapes  of  the  component  performance  maps; 
engine  power  balance;  conservation  of  mass,  energy  and  momentum;  variable 
specific  heat  effects;  nozzle  unchoking  effects,  etc.,  are  taken  into 
account.  For  the  actual  TF30-P-408  engine  used  in  the  test,  such  details 
as  component  efficiency  variations,  specific  shapes  of  the  component 
performance  maps,  etc.,  were  not  .available . Judgment  numbers  based  on 
experience  and  the  model  specification  for  this  engine  were  used  to  obtain 
the  theoretical  fault  coefficients.  The  theoretical  fault  coefficients 
were  computed  at  low,  medium  and  high  power  conditions  at  sea  level  static 
and  a linear  least  square  line  was  fitted  to  the  data. 
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An  alternate  approach  using  actual  engine  data  is  available  for  cal- 
culating the  fault  coefficients.  This  method  is  empirical  and  makes  use 
of  nominal  base  line  engine  data  as  well  as  steady-state  engine  data  with 
deliberate  perturbations  of  known  amounts  Introduced  individually.  From 
test  data  where  all  independent  variables  are  held  fixed  except  one,  the 
resulting  delta  of  any  measurable  dependent  variable  is  a results  of  the 
one  known  perturbation.  In  empirically  determining  the  fault  coefficients, 
it  has  been  tacitly  assumed  that  independent  component  performance  para- 
meters such  as  component  efficiencies,  turbine  nozzle  area,  etc.,  did  not 
independently  change  during  the  test.  In  other  words,  the  only  pertur- 
bations were  exhaust  nozzle  area  changes  and  varying  the  three  compressor 
bleeds.  It  follows,  therefore,  that  only  the  fault  coefficients  with  re- 
spect to  these  four  disturbances  can  be  computed  empirically. 

Once  the  numerical  values  for  the  fault  coefficients  have  been  deter- 
mined, applying  them  to  engine  diagnostics  is  relatively  straightforward. 

In  general,  if  n faults  are  to  be  computed,  n+1  measurements  are  required. 
One  measurement  serves  as  the  abscissa  against  which  n deviations  from 
base  line  values  are  obtained.  If  the  column  vector  jc  represents  the  set 
of  n sought  after  deltas  (such  as  changes  in  component  efficiencies, 
bleed  flows,  geometries,  etc.)  and  the  column  vector  m^^  is  the  set  of  n 
delta  measurements,  then  they  are  related  by  an  nxn  maTrlx  of  fault  co- 
efficients as: 

mi  = Fi 

A subscript  i is  attached  to  m and  F to  emphasize  that,  in  general, 
there  may  be  more  than  one  set  of  sensor  measurements  which  can  be  used  to 
diagnose  a given  set  of  faults.  If  the  measurement  deltas  bear  a mean- 
ingful relationship  to  the  sought  after  deltas,  then  Fj^  will  have  an 
inverse  and  we  can  solve  the  above  equation  for  jc  to  obtain: 

X = (Fl)"^  Bi 

This  equation  represents  the  diagnostic  routine  used  in  this  study. 

Eleven  different  diagnostic  routines  (eleven  different  sets  of  sensors 
and  sought  for  faults)  were  considered  in  this  study  and  a comparison  was 
made  of  their  ability  to  diagnose  a known  set  of  changes  made  to  the  base 
engine  configuration.  Figure  B-2  summarizes  the  eleven  systems  considered 
in  this  study.  There  eleven  systems  by  no  means  represent  all  possible 
systems.  However,  these  eleven  systems  are  sufficient  to  demonstrate  the 
validity  and  versatility  of  Hamilton-Standard' s approach  to  engine  diag- 
nostics . 

Referring  to  Figure  B-2  Systems  1 and  2 use  two  different  sets  of 
sensor  measurements  to  look  for  the  same  set  of  faults.  The  sought  for 
faults  for  Systems  1 and  2 are  precisely  those  which  were  varied  during 
the  test.  Because  of  the  very  limited  amount  of  data  available  with  a 
surge  bleed  (WgL  3.5)  change  introduced,  the  remaining  systems  replace 
the  3.5  calculation  with  one  or  more  additional  independent  parameter 
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changes.  The  magnitudes  of  these  additional  parameter  changes  were  assumed 
to  be  zero.  This  assumption  is  reasonable  since  the  data  was  taken  over  a 
time  period  too  short  for  the  engine  components  to  degrade  significantly. 

It  is  necessary,  of  course,  to  use  the  theoretically  determined  fault 
coefficients  to  relate  the  measurement  deltas  to  these  additional  indepen- 
dent parameter  deltas.  Systems  3,  4,  and  5 have  replaced  the  AWgL  3.5 
calculation  with  AA5.  Systems  6 and  7 have  increased  the  number  of  measure- 
ment deltas  (and  sought  for  deltas)  by  1.  Systems  8-11  have  increased  the 
number  of  measurement  deltas  by  2.  Note  that  Systems  1 and  4 and  Systems 
2 and  3 use  the  same  set  of  measurement  deltas  to  look  for  different  sets 
of  independent  parameter  deltas. 

D.  ANALYSIS  RESULTS 


A Fortran  IV  Computer  Program  for  use  on  the  IBM  370/155  Digital  Computer 
was  written  to  evaluate  the  various  diagnostic  systems  indicated  in  Figure 
B-2.  For  any  of  the  systems,  the  program  establishes  the  matrix  of  fault 
coefficients  using  N1C2  as  the  abscissa;  calculates  the  set  of  measurement 
deltas  at  N1C2  = constant  from  the  actual  snapshot  measurements  and  the 
base  line  performance  curves;  performs  the  required  matrix  inversion  and 
calculates  the  independent  parameter  deltas.  This  procedure  is  repeated 
using  N2C2  and  PT7/PT2  as  the  abscissa.  Theoretically,  the  calculated  in- 
dependent deltas  should  be  the  same  regardless  of  which  parameter  is  used 
as  the  abscissa.  Results  were  obtained  using  both  the  theoretically 
determined  fault  coefficients  as  well  as  the  empirical  fault  coefficients. 
Results  obtained  with  the  latter  are  identified  by  the  term  "mixed." 

1.  Sea  Level  Static 


In  order  to  assess  the  relative  merits  of  the  various  diagnostic 
system  shown  in  Figure  B-2,  several  individual  snapshots  were  passed  through 
all  eleven  systems  and  a comparison  made  of  their  ability  to  predict  the 
engine  configuration. 

Figure  B-3  shows  what  the  predicted  independent  deltas  would  be  for 
each  of  the  eleven  systems  using  the  data  of  Snapshot  114.  For  Snapshot  114, 
the  engine  pressure  ratio  was  approximately  2.0.  The  actual  parameter  deltas 
are  indicated  by  the  solid  horizontal  lines.  As  shown  in  Figure  B-3,  the 
only  change  to  the  base  engine  configuration  was  a +5%  change  in  the  ex- 
haust nozzle  area  Ajfj.  Results  using  both  the  theoretical  as  well  as  the 
mixed  fault  coefficients  are  shown  in  Figure  B-3.  Note  that  for  each 
system,  three  points  are  shown.  The  left  most  point  for  each  system 
gives  the  results  using  N1C2  as  the  abscissa  in  evaluating  the  fault  co- 
efficients and  calculating  the  measurement  deltas.  The  middle  point  for 
each  system  gives  the  results  using  N2C2  as  the  abscissa  and  the  right 
most  point  for  each  system  gives  the  results  using  PT7/PT2  as  the  abscissa. 

All  systems  attempt  to  diagnose  the  changes  in  Ajj,j,  WBL3  and  WBL4. 

Systems  1-6  do  about  equally  well  in  predicting  AAjij.  Depending  on  the 


I 


B-4 


NAPTC-PE-88 


abscissa,  they  calculate  AAjjj  to  be  4,4  to  4.6%,  which  is  quite  close  to 
the  actual  value  of  5%.  System  7 does  a little  better  with  a AAjjj  calcu- 
lation of  about  4,7%.  Figure  B-3  shows  that  System  8 calculates  AAjn  most 
accurately.  Since  the  only  perturbation  deliberately  introduced  was  a +5% 
change  in  Ajjj^  the  actual  levels  of  the  other  independent  deltas  is  taken 
to  be  zero.  Systems  1-7  do  a reasonably  accurate  calculation  of  the  in- 
dependent deltas.  The  results  obtained  using  the  mixed  fault  coefficients 
are  somewhat  better  overall  than  the  results  using  purely  theoretical  co- 
efficients. Systems  8-11  show  significant  differences  between  the  actual 
and  calculated  Independent  parameter  changes.  Systems  10  and  11  for 
example  predict  that  the  3.0  bleed  was  approximately  2.9%  and  2%,  respec- 
tively, even  though  for  this  snapshot  though  for  this  snapshot  the  actual 
bleed  was  zero.  The  calculated  deltas  for  WBL3,  WBL4,  13  and  for 
System  8 using  the  mixed  fault  coefficients  showed  such  a vide  variance 
depending  on  the  abscissa  usea  that  these  points  were  not  plotted  in 
Figure  B-3.  Thus,  although  System  8 does  the  best  job  of  predicting  AAjjj, 
it  predicts  the  other  deltas  so  poorly  that,  based  on  the  results  for  this 
single  snapshot,  it  would  have  to  be  rejected  as  an  acceptable  diagnostic 
system. 

Figure  B-4  shows  the  diagnostic  results  using  the  data  for  Snapshot  115. 
For  this  case,  two  perturbations  were  introduced  simultaneously;  a 5%  in- 
crease, Ajj^  ^nd  a 4.0  bleed  of  approximately  2.2%,  Systems  1-7  using  the 
mixed  matrices  predict  these  changes  quite  well.  For  all  systems,  the 
predicted  AAjjj  using  the  theoretical  fault  coefficients  is  about  0.5  to  1% 
lower  than  the  AAjj^  predicted  using  the  mixed  matrices.  System  8 again  stands 
out  as  an  undeslraole  diagnostic  system. 

Figure  B-5  shows  the  diagnostic  system  results  for  Snapshot  116.  This 
is  a multiple  fault  condition  where  the  actual  change  to  the  base  engine 
configuration  was  a 5%  increase  in  Aj|^  and  a 3.0  bleed  of  about  2.6%. 

Systems  1-7  do  a good  job  of  predicting  these  changes.  Systems  8 and  9 
are  judged  unacceptable  since  they  have  difficulty  distinguishing  between 
the  3.0  and  4.0  bleeds.  Using  the  theoretical  fault  coefficients. 

Systems  10  and  11  predict  the  3.0  bleed  significantly  higher  than  the 
actual  value.  System  10  predicts  a fan  efficiency  change  of  -2.6%. 

System  11  predicts  a fan  efficiency  change  of  -2%.  In  both  cases,  the 
actual  '’an  efficiency  change  is  assumed  to  be  zero.  For  these  reasons. 

Systems  10  and  11  are  also  judged  unacceptable  as  diagnostic  systems. 

Figures  B3-5  show  one  way  of  presenting  the  diagnostic  results. 

figures  show  various  snapshots  passed  through  all  systems.  Such 
■ ; "li  Is  useful  In  qulcxly  isolating  unacceptable  systems.  An 

.•  .-  »sv  1 presenting  the  diagnostic  results  is  to  pass  all  snap- 
■,4'  * ■•Ingle  system.  Rather  than  considering  all  seven  of  the 

. . 'stle  avstems,  four  which  were  judged  best  In  overall  per- 

1 -r*  ir  fh*<i  manner.  Figure  B-6  shows  the  diagnostic 

• s*e«  'or  ;v«;tem  1,  The  abscissa  used  In  the 

s ••  »««ure  ratio  (PI7/PT2).  The  actual  fault 
r • le  greater  visibility  to  the 
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results,  various  snapshots  were  divided  into  six  groups.  The  first  group 
had  just  a +5%  Ajj.  fault  Introduced.  This  group  consisted  of  seven  snap- 
shots. The  second  group  has  eleven  snapshots  and  are  those  cases  that 
had  just  a 3.0  bleed.  The  third  group  had  just  a 4.0  bleed  fault.  The 
total  number  of  snapshots  in  this  group  is  14.  Tlie  fourth  grouping  are 
those  snapshots  which  have  both  a 5%  Aj>^  increase  and  a 4.0  bleed.  The 
fifth  grouping  are  those  snapshots  which  have  both  a 5%  Aj^j  increase  and 
a 3.0  bleed.  The  total  number  of  snapshots  in  each  of  the  latter  two 
groups  is  four.  The  last  group  are  those  snapshots  for  which  no  fault  was 
indicated.  It  was  assumed  that  these  cases  were  run  with  the  base  engine 
configuration  and  therefore  all  independent  deltas  were  taken  to  be  zero. 
The  total  number  of  snapshots  in  this  last  group  is  seven.  Figure  B-6 
shows  Snapshot  97  to  be  a bad  data  point.  Close  examination  of  the  data 
confirmed  the  error.  Figures  B-7,  B-8  and  B-9  show  the  diagnostic  results 
for  all  snapshots  for  Systems  2,  5 and  7,  respectively.  These  figures 
also  indicate  that  Snapshot  97  is  a bad  data  point. 

A total  of  53  snapshots  were  run  at  the  sea  level  static  condition. 

Six  of  the  53  snapshot  cases  were  not  included  in  Figures  B-6  through  B-9. 
TWO  cases  were  with  just  a 5%  Aj^  decrease  and  two  cases  were  with  a 5Z 
Ajjj  decrease  and  bleed.  For  Systems  1 and  2 which  look  for  the  surge 
bleed  change,  the  calculated  deltas  for  the  two  surge  bleed  cases  were 
in  very  close  agreement  with  the  assumed  actual  values  using  the  empirical 
coefficients  but  showd  a very  wide  variance  using  the  theoretical  coeffi- 
cients. The  wide  variation  between  the  assumed  and  calculated  fault  levels 
was  obtained  for  all  the  calculated  faults.  The  reason  for  this  is  that 
the  surge  bleed  is  such  a severe  disturbance  to  the  base  line  configura- 
tion and  produces  such  large  measurement  deltas  that  the  linearization 
procedure  on  which  the  theoretical  fault  coefficients  are  based  is  no 
longer  valid.  For  those  systems  which  did  not  look  for  a surge  bleed 
change,  the  two  surge  bleed  cases  gave  totally  erroneous  results  as  ex- 
pected. How  an  unsought  for  independent  delta  affects  the  calculated 
deltas  is  discussed  in  Reference  A. 

For  the  cases  where  a -5%  change  in  Aj|^  was  made,  the  diagnostics 
systems  calculated  the  change  to  be  between  approximately  -8%  and  -10%. 

This  is  caused  by  the  nonlinear  nature  of  the  problem.  Comparing  the 
dependent  parameter  deltas  for  a +5%  change  in  Aj^  and  a -5%  change  in 
Ajn  showed  that  the  resulting  deltas  were  not  equal  in  magnitude  and 
opposite  in  sign  as  would  be  the  case  if  everything  were  purely  linear. 

2.  Non  Zero  Mach  No.  Conditions 

In  order  to  establish  the  proper  measurement  deltas  at  the 
different  Mach  number  conditions,  two  methods  are  available.  One  method 
uses  the  static  base  line  to  establish  a gross  delta  from  which  the 
effect  due  to  Mach  number  must  be  subtracted  to  obtain  the  net  delta 
for  use  in  the  diagnostic  routine.  Another  method  is  to  use  a base  line 
customized  for  the  particular  Mach  number.  The  latter  approach  is  used 
in  this  study.  Due  to  the  relatively  significant  differences  between 
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Che  empirical  and  the  cycle  deck  derived  Mach  effects,  only  the 
empirical  data  was  used. 

At  Mach  0.25,  the  no  fault  snapshot  cases  were  sufficiently  close 
to  the  static  base  lines  to  allow  the  use  of  the  static  base  line 
characterist.’.cs.  At  Mach  0.5  and  0.75,  however,  custom  base  lines  had 
to  be  obtained.  Using  the  no  fault  snapshot  cases,  various  engine 
parameters  were  plotted  as  a function  of  N1C2,  N2C2  and  PT7/PT2.  A 
polynomial  curve  in  Che  least  square  error  sense  was  fitted  to  the  data. 
These  curves  were  used  to  establish  the  measurement  deltas  at  Mach  0.50 
and  0. 75 . 

Eight  snapshots  were  obtained  at  the  simulated  Mach  = 0.25  condition. 
Three  had  Just  a 5%  Aj^  increase  and  five  had  no  indicated  fault. 

Figures  B-10,  B^ll,  B-12  and  B-13  show  the  diagnostic  results  for  all 
snapshots  for  systems  1,  2,  5,  and  7,  respectively.  The  static  base 
lines  and  the  mixed  fault  coefficients  were  used  to  obtain  the  calculated 
results  shown  in  Figures  B-10  through  B-13.  The  abscissa  used  to 
establish  the  measurement  deltas  and  fault  coefficients  was  engine 
pressure  ratio  (PT7/PT2).  Systems  2 and  7 are  somewhat  better  in  their 
ability  to  calculate  AAjjg,  AWBL3  and  AWBL4  than  systems  1 and  5.  The 
calculated  deltas  for  most  of  the  snapshots  agree  well  with  the  actual 
deltas. 

Using  the  Mach  No.  = 0.50  base  line  characteristics  and  the  mixed 
fault  coefficients.  Figures  B-14  through  B-17  show  the  diagnostic  re- 
sults for  Systems  1,  2,  5 and  7,  respectively.  Snapshots  175  and  176 
were  the  only  non-zero  Mach  number  conditions  which  had  multiple  faults 
introduced.  Figures  B-14  through  B-17  show  the  calculated  deltas  to  be 
very  close  to  actual  deltas  for  these  cases. 

Using  the  Mach  No.  = 0.75  base  line  characteristics  and  the  nixed 
fault  coefficients.  Figures  B-18  through  B-21  show  the  diagnostic  re- 
sults for  Systems  1,  2,  5 and  7 for  the  simulated  Mach  = 0.75  snap- 
shots. The  calculated  deltas  for  the  no  fault  conditions  agree  quite 
well  with  the  actual  faults.  Only  three  snapshots  were  obtained  with 
just  5%  area  increase.  In  two  out  of  these  three  snapshots  the  cal- 
culated change  in  WBL3  is  about  +1%  compared  to  the  true  value  of  0. 
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Of  the  different  diagnostic  systems  considered.  System  7 was  chosen 
as  the  primary  diagnostic  system.  Although  it  does  not  calculate  the 
surge  bleed  change,  it  does  calculate  the  changes  in  Ay^j,  WBL3  and  WBL4 
quite  well.  In  addition  to  these  three  primary  perturbations  introduced. 
System  7 also  calculates  the  change  in  the  gas  producer  turbine  effective 
nozzle  area  (AA5)  and  also  the  change  in  the  gas  producer  turbine 
efficiency  (Ahth) • for  most  of  the  97  snapshots  the  calculated  deltas 
in  A5  and  q-jqy  were  very  close  to  the  assumed  values  of  0. 

System  2 is  chosen  as  an  alternate  to  System  7.  This  system  is 
simpler  in  that  it  only  considers  the  actual  perturbations  Introduced. 
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This  sytem  also  calculates  the  changes  in  Aj^,  WBL3  and  WB14  quite  well. 

A comment  sometimes  encountered  concerning  Hamilton  Standard's  tech- 
nique for  multiple  fault  engine  diagnostics  is  that  excessive  computer 
storage  is  required  to  store  and  interpolate  the  fault  coefficient 
matrices  and  to  do  the  necessary  matrix  inversion  to  obtain  the  required 
diagnostic  matrix.  This  of  course  is  not  the  case.  It  is  important  to 
realize  that  all  the  matrix  inversions  necessary  to  establish  the  final 
diagnostic  matrix  are  done  offline.  Only  the  final,  inverted,  diagnostic 
matrices  need  be  stored.  Secondly,  it  is  usually  not  necessary  to  store 
more  than  one  diagnostic  matrix.  In  a practical  diagnostic  system,  the 
actual  diagnostics  would  be  done  within  a specified  power  band.  How  wide 
this  band  is  depends  on  the  particular  application.  Usually,  a single 
diagnostic  matrix  evaluated  at  the  mid  power  condition  within  this  band 
is  all  that  is  required  to  do  the  engine  diagnostics. 

The  required  memory  increment  to  include  the  gas  path  analysis  of 
System  7 in  the  computer  of  a total  engine  monitoring  system  would  be 
quite  minimal.  If  we  assume  all  base  lines  have  a common  abscissa,  then 
approximately  152  sixteen  bit  word  storage  capacity  would  be  required  to 
do  the  following: 

a.  To  store  5 base  lines  at  three  Mach  numbers.  It  is  assume 
that  the  base  lines  can  be  adequately  represented  by  5 common  abscissa 
values  for  each  base  line.  Using  a 16  bit  work  to  represent  two  points, 

40  words  are  required. 

b.  To  store  the  elements  of  a single  7x5  diagnostic  matrix  D. 

For  precision,  35  sixteen  bit  words  are  used. 

c.  To  calculate  and  store  the  fault  vector  x by  performing  the 
matrix  multiplicate  x = D m,  given  the  measurement  vector  m. 

Note  that,  as  a kind  of  bonus,  the  percentage  change  in  turbine  in- 
let temperature  and  referred  specified  fuel  consumption  measured  against 
the  assumed  abscissa  N2C2  are  also  computed  in  this  analysis.  Since 
these  two  parameters  are  dependent,  their  changes  can  be  calculated  on 
the  basis  of  the  changes  in  the  independent  variables.  This  is  an  auto- 
matic fallout  of  the  gas  path  analysis  technique.  The  calculated  change 
in  T5C2  and  SFC  will  be  the  change  that  would  be  observed  if  N2C2  were 
held  fixed.  This  does  not  imply  that  N2C2  must  be  used  as  an  abscissa 
in  establishing  the  measurement  deltas.  The  assumed  abscissa  for  T5C2 
and  SFC  need  not  be  the  same  for  both,  nor  do  they  necessarily  have  to 
be  the  same  as  the  abscissa  used  in  establishing  the  measured  parameter 
deltas. 

It  should  be  noted  that  WBL3  and  WBL4  faults  are  relatively  small 
(2-3%  range).  In  the  sense  that  bleed  is  analogous  to  compressor  or  tur- 
bine degradation  in  adversely  affecting  engine  performance,  the  bleed  levels 
detected  and  diagnosed  in  this  program  were  equivalent  to  relatively  minor 
degradation  problems  elsewhere  in  the  engine. 
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FIGURE  B-2:  POSSIBLE  GAS  PATH  ANALYSIS  SYSTEMS  FOR 

NAPTC  TF30-P-408  DATA  ANALYSIS 
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FIGURE  B-3:  DIAGNOSTIC  SYSTEM  PREDICTION  RESULTS  FOR  SNAPSHOT  114 


A'^BL4.oI 


|-”Tr  ” 

...  I:.. 

— *-H 
1 } 

!- aa 

JN 

!T:  h 

; 

r • ■■■ 

j. 

; 1 

! 

AWbu  3.4 

|- 

... 

...  t 

r •!  • 

i — Lu- 

A''''bL3:5 

-1 — 

— I 

2 L-_3 

.1 

r* 

: 

4 

. 1 

L . 

*■ 

. 



I 

r 

--  i 

i—  -- 

...r 

^ ^ 

i 

! 

* 

. j 

. t 

« 

k I 

F-  ' 

i ■ 

1 

r 

i 

i 

L 

J 

•p'*  ■ " 

t*— *1 

i 

I 

i 

- : —I 

- 1 

i 

i 

- .M 

■ 

: ■ ■■ 

■ - 1 

1 

. 1 

: ' 

. .j 

I - 

..T4 

..  ...j 

. I 

i. ... 

. 

L-  . . 

1 

• r 

■ t 

■;  t 



■ 

mm 

iM 

■ 

m 

aHnnR 

t ■' 

“ 

, 

3 



, . I . ..  i 

iR  SNAPSHi 


Ai/Uitjo 


2 A 

1 



■ ^ 

' * ! 

H 

I — - — , 

*1* — r 

F — H 

h--1 

U— • • 

_ - . 

- 1 

I ! 

' J 

I ■ 

j 

; ' 

1 ■ 

' i 

■ I 

! » , 

i 

i 

i—  J 

1 

1 

1 ■ ■ ; 

xl  \ 

i i 

[ 

r 1 

'■  ; 

■ < 

■; 

; 

1 I 

- - -j 

j 

I 

; 

r 1 

! 

0 - 

1 ■ ' . 

hi  ^ 

h*  ' 

■ -j 

1 

N ! 

i 

• J 

' ■ 1 

.1 

^ i 

L-  •] 

h 

I ! 

! I 

• i 

■ 4 

■| 

1 

• - 

i 

1 

1 

’ - 

■n  ■ 

i 

_j 

1 — 

1 

i 

! 

i ’ 

...J 

_._i 

i 

1 

* i 

..  r L . ' " 

i 

, 1 

' j 

■ i 

1 s 

} 

r* 

V - ! 

-i 

•4 

^ u 

L.. 

r* 

FT"H 

0 SI 

1 

■ ■ i 
u J 

’■*■  I 

Lfr  . . 1 

H 

! 

i 

2 A 

1 ■ ! 

..-J 

i i 

1 

1 

.... . j 

: "-J 

-J 

! 

1 

4 

j 

} 

! 

j 

FIGURE  B-6:  SYSTEM  1 TF30-P-408  MN=0  DIAGNOSTIC  RESULT 
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FIGURE  B-7:  SYSTEM  2 TF30-P-408  MN=0  DIAGNOSTIC  RESULTS 


FIGURE  B-8:  SYSTEM  5 TF30-P-408  MN=0  DIAGNOSTIC  RESULTS 
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FIGURE  B-12:  DIAGNOSTIC  RESULTS  SYSTEM  5 TF30-P-408  MN=0,25 
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FIGURE  B-13:  DIAGNOSTIC  RESULTS  SYSTEM  7 TF30-P-408  MN=0.25 
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FIGURE  B-14:  DIAGNOSTIC  RESULTS  SYSTEM  1 TF30-P-408  MN=0.5 
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FIGURE  B-15:  DIAGNOSTIC  RESULTS  SYSTEM  2 TF30-P-408  MN=0.5 
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FIGURE  B-17:  DIAGNOSTIC  RESULTS  SYSTEM  7 TF30-P-408  MN=0,5 
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FIGURE  B-21:  DIAGNOSTIC  RESULTS  SYSTEM  7 TF30-P-A08  MN=C.75 
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APPENDIX  C 

Variable  Geometry  Cycle  Deck  - Description  and  Status 

A.  GENERAI,  DESCRIPTION 

This  cycle  deck  has  been  written  to  conforc  to  Aerospace  Standard 
681-C,  which  addresses  Gas  Turbine  Engine  Steady  State  Performance 
Presentation  for  Digital  Computer  Programs,  and  Aerospace  Recommended 
Practice  755-A,  concerning  Gas  Turbine  Engine  “erformance  Station 
Identification  and  Nomenclature.  The  cpeclfJc  gas  turbine  cycle 
simulated  by  this  deck  is  a two-spool,  T-fan^,  mixed-flow  turbofan 
with  afterburning.  In  addition,  this  cycle  deck  allows  for  variable 
area  high  and  low  pressure  turbine  (HPT  and  LPT)  nozzles  and  variable 
exhaust  nozzle  (A8)  area.  Figure  C-1  is  a schematic  diagram  of  the 
engine  components  showing  the  station  identification  numbers. 

B.  OVERALL  OPERATION 


Basically,  this  cycle  deck  consists  of  the  thermodynamic  equations 
of  a two-spool,  mixed-flow  turbofan.  These  equations  account  for 
flow  continuity,  energy  balance,  speed  matching,  and  pressure  balance. 

In  addition,  the  program  is  required  to  satisfy  the  requested  HPT,  LPT, 
and  A-8  areas.  The  theromdynamlc  properties  of  air  are  obtained  from 
GASTAB,  which  is  a computerized  gas  table. 

C.  INPUTS 

1.  Design  Inputs 

This  cycle  deck  has  been  written  as  a design  or  Parametric  deck, 
so  that  any  engine  in  the  two-spool,  T-fan,  mixed-flow  turbofan  category 
can  be  simulated.  To  simulate  a particular  engine,  the  engine  design 
point  information  must  be  provided  to  the  program,  operating  in  the 
Design  Mode.  The  design  point  can  be  at  any  desired  flight  condition 
and  must  be  at  intermediate  power,  i.e.,  rating  code  = 50.  Table  C-I 
shows  the  required  design  point  inputs. 

Once  a design  run  has  been  balanced,  that  is  adjusted  to  obtain  a 
static  pressure  balance  at  the  mixer  and  adjusted  to  obtain  the 
desired  performance,  a design  deck  is  generated.  This  deck  contains 
the  "standard"  inputs  and  the  design  point  information  required  to  use 
the  deck  in  its  Normal  Mode. 

2.  User  Inputs 

Operating  in  its  Normal  Mode,  the  program  accepts  inputs,  such 
as  altitude  and  Mach  number,  horsepower  and  bleed  air  extraction. 


1 T-fan  implies  fan  and  low  pressure  compressor  on  low  spool. 
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changes  in  component  performance,  changes  in  HPT,  LPT,  and  exhaust 
nozzle  areas,  and  rating  code  (ZRC) , which  is  the  power  level  input. 
The  program  can  be  operated  to  generate  a single  data  point  or  to 
generage  a "power  series"  at  the  given  input  conditions.  The  "power 
series"  is  controlled  by  a power  code  (ZPC)  input.  In  all  cases,  a 
new  input  will  supersede  the  standard  input,  and  will  remain  until 
changed  or  until  the  standard  deck  is  reloaded. 

TABLE  C-I 

Design  Point  Inputs 


1. 

T2* 

Fan  Inlet  Temperature 

2. 

P2* 

Fan  Inlet  Total  Pressure 

3. 

TAMB* 

Ambient  Temperature 

4. 

PAMB* 

Ambient  Pressure 

5. 

ZALT* 

Altitude 

6. 

ZMN* 

Flight  Mach  Number 

7. 

VA* 

Flight  Airspeed 

8. 

SERAM 

Inlet  Recovery  Computation  Method 

9. 

ZRR 

Inlet  Recovery  (used  only  if  SERAM  = 1) 

10. 

SIM 

Inlet  Condition  Input  Selector 

11. 

W2 

Inlet  Airflow 

12. 

BPR 

By-pass  Ratio 

13. 

P13Q2 

Outer  Fan  Pressure  Ratio 

14. 

EFO 

Outer  Fan  Efficiency 

15. 

P21Q2 

Inner  Fan  Pressure  Ratio 

16. 

EFI 

Inner  Fan  Efficiency 

17. 

P23Q21 

Low  Pressure  Compressor  Pressure  Ratio 

18. 

ECL 

Low  Pressure  Compressor  Efficiency 

* Not  all  are  required  - value  of  SIM  dictates  requirements 
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19. 

P3Q23 

High  Pressure  Compressor  Pressure  Ratio 

20. 

ECH 

High  Pressure  Compressor  Efficiency 

21. 

WB22Q 

Customer  Bleed  Air  Extraction,  LPC  Dishcharge 

22. 

WB29Q 

Customer  Bleed  Air  Extraction,  HPC  Discharge 

23. 

T4 

Burner  Discharge  Temperature 

24. 

EB 

Burner  Efficiency 

25. 

FHV 

Fuel  Heating  Value 

26. 

ETH 

High  Pressure  Turbine  Efficiency 

27. 

ETL 

Low  Pressure  Turbine  Efficiency 

28. 

WBHVIQ 

Turbine  Cooling  Flow,  HPT  First  Vane 

29. 

WBHBIQ 

Turbine  Cooling  Flow,  HPT  First  Blade 

30. 

WBHV2Q 

Turbine  Cooling  Flow,  HPT  Second  Vane 

31. 

WBHB2Q 

Turbine  Cooling  Flow,  HPT  Second  Blade 

32. 

WBLVIQ 

Turbine  Cooling  Flow,  LPT  First  Vane 

33. 

^LBIQ 

Turbine  Cooling  Flow,  LPT  First  Blade 

34. 

WBLV2Q 

Turbine  Cooling  Flow,  LPT  Second  Vane 

35. 

WBLB2Q 

Turbine  Cooling  Flow,  LPT  Second  Blade 

36. 

TMHVl 

Maximum  Allowable  Turbine  Metal  Temperature, 
HPT  First  Vane 

37. 

TMHBl 

Maximum  Allowable  Turbine  Metal  Temperature, 
HPT  First  Blade 

38. 

TMHV2 

Maximum  Allowable  Turbine  Metal  Temperature, 
HPT  Second  Vane 

39. 

TMHB2 

Maximum  Allowable  Turbine  Metal  Temperature, 
HPT  Second  Blade 

40. 

TMLVl 

Maximum  Allowable  Turbine  Metal  Temperature, 
LPT  First  Vane 

41. 

TMLBl 

Maximum  Allowable  Turbine  Metal  Temperature, 
LPT  First  Blade 
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43. 

TMLB2 

Maximum  Allowable  Turbine  Metal  Temperature, 
LPT  Second  Blade 

44. 

PF4 

Burner  Discharge  Pattern  Factor 

45. 

RPFH 

High  Pressure  Turbine  Rotating  Pattern  Factor 

46. 

RPFL 

Low  Pressure  Turbine  Rotating  Pattern  Factor 

47. 

T7MAX 

Maximum  Afterburner  Discharge  Temperature 

48, 

EAB 

Afterburner  Efficiency 

49. 

XMAB 

Afterburner  Mach  Number 

50. 

CDS 

Exhaust  Nozzle  Discharge  Coefficient 

51. 

A8 

Exhaust  Nozzle  Throat  Area 

52. 

XNHT 

H.gh  Pressure  Rotor  Speed 

53. 

XNLT 

Low  Pressure  Rotor  Speed 

54. 

PWXH 

High  Pressure  Rotor  Power  Extraction 

55. 

PWXL 

Low  Pressure  Rotor  Power  Extraction 

56. 

DP16 

Fan  Duct  Pressure  Drop 

57. 

DPB 

Burner  Pressure  Drop 

58. 

DP56 

Turbine  Diffuser  Pressure  Drop 

59. 

DPAB 

Afterburner  Friction  Pressure  Drop 

60. 

P61D6 

Mixing  Pressure  Drop 

61. 

XMN16 

Fan  Duct  Mach  Number 

62. 

XMN6 

Turbine  Diffuser  Mach  Number 

3.  Variable  Geometry  Schedules 

Another  group  of  inputs  are  used  to  define  the  HPT,  LPT, 
and  A-8  area  schedules.  The  area  schedules  are  expressed  as  a function 
of  rating  code,  and  are  at  this  point  limited  to  three  (3)  linear  seg- 
ments, as  shown  in  Figure  C-2.  Once  a set  of  schedules  has  been 
determined,  these  schedules  can  be  Incorporated  into  the  "standard" 
or  design  deck  for  convenient  re-use  at  a later  time.  Any  input  area 
schedules  will  supersede  the  standard  area  schedules. 
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D.  OUTPUT 


The  data  generated  by  this  cycle  deck  is  output  to  the  line 
printer  and/or  magnetic  tape.  There  are  three  (3)  different  printer 
formats  for  a single  data  point.  Format  1 is  a complete  listing  of 
all  Inputs,  all  stored  computed  values,  and  a summary  page  of  engine 
parameters  at  each  station.  The  second  format  is  just  the  summary 
page.  The  third  format  is  just  the  Inlet  Conditions  and  Overall 
Performance  sections  of  the  summary  page  (Format  2 is  the  standard 
format) . 

In  addition  to  the  engine  performance  data,  the  summary  page 
includes  a matrix  of  tlumerical  Status  Indicators  (NSI's)  which  are 
warning  messages.  These  warning  messages  Include  such  things  as 
improper  or  inconsistent  inputs.  Iteration  non-convergence,  flight 
envelope  exceedance,  internal  function  or  map  exceedance,  compressor 
surge,  and  insufficient  or  excess  cooling  flow  messages  (see  Table  C-II) 

TABLE  C-II 

Numerical  Status  Indicator  Listing 
NS I MESSAGE 

0200  Illegal  Altitude  Input,  Altitude  reset  to  0. 

0201  Illegal  Mach  Number  Input,  Mach  number  reset  to  0. 

0202  Illegal  Flight  Speed  Input,  Flight  Speed  reset  to  0. 

0203  Illegal  Ram  Recovery  Input,  Ram  Recovery  reset  to  1. 

0204  Illegal  Ram  Recovery  Selector,  Ram  Recovery  Selector  reset  to  0 
0300  Illegal  Power  Code  Input,  Power  series  not  executed 

0500  Flight  Envelope  Exceedance 

070X  Excess  Cooling  Flow,  at  station  4,  (X-1) 

1801  Exhaust  Nozzle  Function  limit  exceedance 

1802  Fan  Efficiency  Function  limit  exceedance 

1803  Compressor  Efficiency  Function  limit  exceedance 
270X  Insufficient  Cooling  Flow,  at  station  4.  (X-1) 

7001  High  Pressure  Turbine  Map  Exceedance,  Velocity  Ratio  >4.0 
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NS I MESSAGE 

7002  High  Pressure  Turbine  Map  Exceedance,  Work  Ratio  >2.0 

7003  High  Pressure  Turbine  Map  Exceedance,  Referred  Speed 
Ratio  <.897 

7004  High  Pressure  Turbine  Map  Exceedance,  Referred  Speed 
Ratio  >1.1J 

7011  Low  Pressure  Turbine  Map  Exceedance,  Velocity  Ratio  >4.0 

7012  Low  Pressure  Turbine  Map  Exceedance,  Work  Ratio  >2.0 

7013  Low  Pressure  Turbine  Map  Exceedance,  Referred  Speed 

Ratio  <.'65 

7014  Low  Pressure  Turbine  Map  Exceedance,  Referred  Speed 
Ratio  >1.17 

7051  Fan  Map  Exceedance,  Referred  Flow  Ratio  >1.10 

7052  Fan  Map  Exceedance,  Referred  Flow  Ratio  <.35 

7053  Fan  Map  Exceedance,  Percent  Fan  Pressure  Ratio  >1.2 

7054  Fan  Map  Exceedance,  Percent  Fan  Pressure  Ratio  <.35 

7061  Low  Pressure  Compressor  Map  Exceedance,  Referred  Flow 
Ratio  >1.2 

7062  Low  Pressure  Compressor  Map  Exceedance,  Referred  Flow 
Ratio  <.45 

7063  Low  Pressure  Compressor  Map  Exceedance,  Percent  LPC 
Pressure  Ratio  >1.2 

7064  Low  Pressure  Compressor  Map  Exceedance,  Percent  LPC 
Pressure  Ratio  <.45 

7071  High  Pressure  Compressor  Map  Exceedance,  Referred 
Flow  Ratio  >1,2 

7072  High  Pressure  Compressor  Map  Exceedance,  Referred 
Flow  Ratio  <.45 

7073  High  Pressure  Compressor  Map  Exceedance,  Percent  HPC 
Pressure  Ratio  >1.2 
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NSI  MESSAGE 

7074  High  Pressure  Compressor  Map  Exceedance,  Percent  HPC 
Pressure  Ratio  <.45 

7091  Fan  Surge 

7092  Low  Pressure  Compressor  Surge 

7093  High  Pressure  Compressor  Surge 

7100  Inlet  Referred  Flow  > Maximum  limit 

7101  Fan  Discharge  Mach  Number  >.S 

7110  High  Pressure  Turbine  Pressure  Ratio  > Maximum  limit 

7111  Low  Pressure  Turbine  Pressure  Ratio  > Maximum  limit 

9101  Non-convergence,  A-8  iteration:  DATA  INVALID 

9102  Non-convergence,  FP4,  FP44  iteration:  DATA  INVALID 

9200  Illegal  Inlet  Conditions  Input,  Calculated  Ram 

Recovery  >1.:  RUN  NOT  EXECUTED 

9201  Illegal  Inlet  Mode  Selection  Input:  RUN  NOT  EXECUTED 

9300  Illegal  Rating  Code  Input:  RUN  NOT  EXECUTED 

9301  Illegal  Power  Code  Input:  RUN  NOT  EXECUTED 

E.  FEATURES 


Some  of  the  features  of  this  program  which  make  it  flexible  are 
as  follows.  Each  compressor  and  turbine  are  Individually  represented 
by  performance  maps.  The  compressor  maps  are;  (a)  referred  speed  as 
a function  of  pressure  ratio  and  referred  flow  and,  (b)  efficiency  as 
a function  of  pressure  ratio  and  referred  flow.  The  turbine  maps 
are;  (a)  flow  parameter  as  a function  of  pressure  ratio  and  referred 
speed  and,  (b)  efficiency  as  a function  of  velocity  ratio  and  cor- 
rected work.  These  are  stored  as  percentages  of  design  so  that  the 
same  maps  can  be  used  for  many  engines.  Pressure  drops  are  not 
considered  to  be  constant.  The  fan  duct,  burner,  diffuser,  mixer, 
and  afterburner  friction  pressure  drops  are  computed  as  a function  of 
the  referred  flow  and  the  design  value  of  the  pressure  drop,  providing 
more  realistic  calculations.  Also,  the  afterburner  momentum  pressure 
loss  is  computed  as  a function  of  afterburner  Mach  number  and  temper- 
ature rise.  Whenever  possible,  the  program  coding  was  organized  into 
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building  blocks,  that  is  a subroutine  for  a compressor,  a burner,  a 
turbine  vane,  a turbine  blade,  etc.  This  has  been  done  to  facilitate 
construction  of  other  engine  cycles  without  duplicate  effort. 

For  purposes  of  identification,  a program  indicator  code  (IND)  is 
included  with  the  data.  This  six  digit  code,  XXYYZZ  is  used  to 
identify  the  cycle  (XX)  (i.e.,  two-spool,  T-fan,  mixed-turbof an 
augmented  = 01),  the  specific  design  (YY),  and  the  variable  geometry 
schedules  (ZZ)  used  to  generate  this  data. 

Another  option  within  tin's  program  involves  the  number  of  turbine 
stages.  Both  the  hlgii  and  low  pressure  turbines  have  the  option  of 
being  single  or  two-stage  components  Independently. 

F.  PERIPHERAJ,  CAPABILITIES 

Once  data  has  been  written  on  magnetic  tape,  there  are  a number 
of  ways  to  retreive  desired  information.  One  option  is  to  list  all 
runs  on  tape.  The  output  of  this  option  is  one  line  printer  line 
per  run  on  tape  showing  the  Case  (or  Run)  number,  the  program  indi- 
cator code,  the  day,  date  and  time  the  run  was  generated,  and  the 
NSI  array.  It  is  also  possible  to  output  a chosen  run  or  set  of 
runs  in  any  of  the  three  data  output  formats. 

Another  tape  feature  which  is  available  is  called  Tape  Edit.  This 
allows  output  of  up  to  eight  stored  values  (Inputs  or  outputs)  for  all 
or  any  selected  group  of  runs.  The  "editing"  is  accomplished  by 
requesting  that  only  runs  which  satisfy  the  editing  criteri?  be 
output,  and  any  stored  value  may  be  used  as  an  editing  criteria.  Each 
editing  criterlal  must  be  assigned  a minimum  and  maximum  value.  For 
example,  if  there  exists  a tape  of  power  series  of  many  different 
engine  designs,  each  run  at  many  altitudes  and  Mach  numbers,  the  user 
may  request  data  within  a specified  band  of  altitudes,  and  within  a 
specified  range  of  thrusts.  Up  to  five  such  editing  criteria  are 
available.  In  addition,  editing  on  the  maximum  value  of  the  NSI 
array  and  the  program  indicator  code  is  available. 

G.  FURTHER  EFFORTS 

The  building  block  approach  to  the  coding  of  this  deck  allows  for 
the  creation  of  other  gas  turbine  cycles  with  a minimum  effort.  All 
of  the  programming  required  for  Input/output  and  program  control  will 
remain  essentially  unchanged.  The  peripheral  capabilities  are  indep- 
endent of  the  cycle  and  are  therefore  available  for  use  on  an  • cycle. 

The  cycles  which  are  planned  to  be  developed  in  the  immediate 
future  are  the  (a)  two-spool,  unmixed-turbof an  with  duct  burning  and/ 
or  afterburning,  (b)  two-spool,  mixed  turbofan  with  afterburning, 

(c)  two-spool,  regenerative  turboshaft  with  a free  power  turbine,  and, 

(d)  high  bypass  ratio  turbofan. 
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Other  features  which  are  planned  to  be  incorporated  into  this 
deck  and  all  other  cycle  decks  are  (a)  English  or  Metric  units, 

(b)  static  pressure  mismatching  (for  mixed  flow  systems)  and, 

(c)  automated  sensitivities.  Automatea  sensitivities  implies  that 
a peripheral  program  will  be  developed  which  can  compare  any  single 
run  on  tape  with  any  defined  operating  line  (whose  data  is  also  on 
tape).  This  method  of  comparison  is  chosen  by  the  user.  For 
example,  this  program  can  be  used  to  generate  the  sensitivity  of 
specific  fuel  consumption  at  constant  thrust  of  a particular  design 
to  a 1%  change  in  a component  efficiency. 

H.  VARIABLE  GEOMETRY  DATA  FOR  PARAMETER  INTERRELATIONSHIPS  STUDY 


Variable  geometry  data  from  this  cycle  deck  was  supplied  to 
Hamilton  Standard  Division,  United  Technologies  Corporation  for 
study  using  the  parameter  interrelationships  method  of  gas  path  analysis. 

This  data  represents  an  engine  whose  geometry  schedules  are  designed 
for  30,000  ft.,  Mach  0.8  Cruise  condition  because  the  cruise  condition 
is  an  appropriate  choice  for  steady  state  gas  path  analysis.  The 
geometry  schedules  were  selected  to  provide  Installed  performance 
gains.  This  is  accomplished  by  maintaining  the  inlet  design  airflow 
while  reducing  thrust  (thereby  reducing  or  eliminating  spillage  drag) 
and  by  opening  the  exhaust  nozzle  to  reduce  boattail  losses  (aft  end 
losses).  Additionally,  maintaining  overall  compression  ratio  tends  to 
minimize  uninstalled  TSFC.  The  operating  limits,  by  which  the  geometry 
schedules  were  restricted,  are  the  compressor  surge  limits  (15% 
reduction  in  surge  margin),  maximum  turbine  pressure  ratio  (P^/P44, 
P44/P5), 

maximum  inlet  referred  flow  (W2R) , maximum  burner  discharge 
temperature  (T4) , and  maximum  fan  discharge  Mach  number. 

Figures  C-8  through  C-12  are  plots  of  the  engine  performance  data 
showing  the  fixed  geometry  and  variable  geometry  operating  lines 
plotted  versus  net  thrust.  The  variable  geometry  data  covers  the 
range  of  40  to  100%  of  intermediate  net  thrust.  Figure  C— 8 shows 
that  constant  inlet  referred  flow  is  held,  accomplishing  reduced 
spillage  drag.  Figure  C-9  shows  uninstalled  TSFC  improvement  from 
70  to  100%  thrust,  while  below  70%  FN,  uninstalled  TSFC  penalties 
are  present.  However,  those  lower  power  conditions  are  associated 
with  the  largest  installation  losses,  so  a net  TSFC  Improvement  may 
be  realized.  Figure  C-13  shows  the  geometry  schedules  used  to 
generate  this  data. 
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A.  INTRODUCTION 


This  studv  was  conducted  on  the  TEDD  Phase  II  program  to  determine 
the  effects  oi  parameter  errors  on  the  diagnostic  capabilities  of  the 
system  and  to  get  an  indication  of  which  parameters  are  the  most 
critical  and  what  measurement  accuracies  are  required.  The  parameters 
stuuies  were  the  ones  used  in  the  steady-state  performance  analysis, 
plus  turbine  blade  temperature  (pyrometer  output) , whicn  was  used  in 
the  turbine  life  subroutines.  The  other  TEDD  parameters  were  used 
mostly  in  simple  limit  exceedance  checks  (such  as  oil  temperature  and 
pressure,  fuel  pressure,  etc.),  or  for  on-off  type  indications  (fuel 
control  selector  switch,  squat  switch,  12th  stage  bleed,  etc.),  where 
extrem.e  accuracy  is  not  as  important. 

B.  MEASUREMENT  ERROR 

Before  presenting  the  results  of  the  sensitivity  study,  a brief 
discussion  of  parameter  measurement  error  is  in  order.  Measurement  error 
has  two  types  of  components,  fixed  error  and  random  error.  Random  error 
(non-repeatability)  is  seen  in  repeated  measurements,  which  do  not  and 
are  not  expected  to  agree  exactly,  for  various  reasons.  The  variation 
between  repeated  measurements  is  called  "precision  error".  The  standard 
deviation  is  used  as  a measure  of  the  precision  error.  A large  standard 
deviation  means  large  scatter  in  the  measurements.  Data  smoothing  tech- 
niques can  be  used  to  minimize  the  effects  of  random  error.  A statistical 
analysis  of  TEDD  Phase  II  data  will  be  presented  later. 

The  constant  or  systematic  error  is  known  as  bias.  In  repeated 
measurements,  each  measurement  has  the  same  bias.  Bias  ran  be  cat- 
egorized into  various  classes.  Known  biases  may  be  eliminated  by 
calibration,  but  unknown  biases  are  not  correctlble.  A comprehensive 
study  of  m'asurement  error  in  gas  turbine  engines  is  provided  in 
Reference  D-1. 

The  absolute  level  of  accuracy  becomes  relatively  unimportant  when 
a differential  or  trending  technique  is  used  for  diagnosis.  In  this 
approach,  a custom  baseline  is  obtained  for  each  engine  using  the 
actual  engine  sensors  and  total  diagnostic  system,  and  the  deltas  or 
changes  from  the  original  values  are  used  for  diagnosis.  This  method 
eliminates  the  effect  of  bias  errors,  as  well  as  engine  to  engine 
variations . 

The  above  (delta)  approach  was  used  for  the  TEDD  performance 
analysis.  The  absolute  measurement  accuracy  is  important,  however,  for 
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other  applications,  such  as  in  critical  limit  exceedance  checks  (RPM, 
turbine  inlet  temperature)  and  in  computations  such  as  the  turbine 
I.  life  subroutines. 

C . TURl^lNE  LIFE  SUBROUTINES 

The  first  sensitivity  analysis  per^^orraed  was  for  turbine  blade 
temperature,  as  utilized  in  the  TEDD  turbine  erosion  life  and  turbine 
creep  life  dlacjncstlc  subroutines.  Calculations  of  turbine  life  were 
made  at  several  different  temperature  levels  and  with  plus  or  minus 
errors  (variations)  from  these  levels.  For  each  temperature  variation, 
the  change  in  calculated  turbine  life  was  expressed  as  a percentage 
(plus  or  minus)  change  from  the  original  level.  The  results  are  shown 
in  Figures  D-1  and  D-2.  It  can  be  seen  that  these  computations  are 
quite  sensitive  to  any  error  in  turbine  blade  temperature.  For  erosion 
life,  this  sensitivity  is  as  high  as  1.2%  computation  error  per  degree 
temperature  error,  or  43%  increase  in  calculated  turbine  erosion  life 
at  -40*F  temperature  error.  The  creep  life  calculations  is  even  more 
sensitive  to  temperature  error,  up  to  a maximum  of  56%  increase  in 
calculated  turbine  life  at  -15“F  temperature  error.  These  calculated 
turbine  life  values  include  a factor  of  safety  of  70%  for  creep  life 
and  90%  for  erosion  life.  Although,  they  have  not  been  evaluated 
experimentally  to  compare  with  the  actual  engine  part  conditions,  they 
Illustrate  some  of  the  difficulties  involved  in  making  accurate  tur- 
bine life  calculations.  Besides  the  problem  of  measurement  accuracy 
(absolute  accuracy),  there  are  temperature  variations  from  blade  to 
blade  and  gradients  along  each  blade.  The  creep  life  calculation  also 
has  an  rpm  input,  but  this  can  be  measured  more  accurately,  and  any 
contribution  to  the  life  computation  error  has  been  found  to  be  rela- 
tively small. 

D.  STEADY-STATE  PERFOR>b\NCE  ANALYSIS 

The  major  portion  of  this  sensitivity  study  dealt  with  the  TEDD 
Steady-State  Performance  Analysis.  This  study  was  made  using  the  TF30- 
P-408  engine  performance  deck  baselines  and  deck  fault  data  processed 
through  the  TEDD  Performance  Analysis  Program.  For  each  of  the  eleven 
diagnosable  faults,  data  was  processed  with  errors  input  to  each  of  the 
twelve  different  performance  parameters  used,  in  +1%  increments  to  +5%. 
Errors  greater  than  +5%  were  not  considered,  since  they  are  larger  than 
normal  parameter  inaccuracies,  and  in  most  cases,  can  be  detected  by  the 
instrumentation  check  portion  of  the  program.  In  fact,  this  instru- 
mentation check  had  to  be  removed  from  the  program  for  this  study,  since 
it  could  even  detect  smaller  errors. 

The  engine  operating  condition  chosen  was  sea  level  static, 
standard  day,  with  a power  setting  of  2200°R  turbine  inlet  temperature, 
which  is  approximately  midway  in  the  engine  pressure  ratio  (ERP)  band 
where  performance  analysis  is  done  (ERP  = 1.6  to  maximum  power). 
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Table  D-I  identifies  the  eleven  faults  analyzed  and  the  fault  levels 
used  for  this  study.  Table  D-Il  lists  the  measured  performance  para- 
meters used  in  the  diagnosis. 


IPCE 

HPCE 

HPTE 

LPTE 

dPTA 

LPTA 

ANZ+ 

ANZ- 

BL09 

BL12 

BL16 


TABLE  D-I 

Performance  Fault  Identlf i.cation  and  Level 

= Lew  Pressure  Compressor  Efflclencj-  (5%  Decrease) 

= high  Pressure  Compressor  Efficiency  (6%  Decrease) 
= High  Pressure  Turbine  Efficiency  (6.?  Decrease) 

= Low  Pressure  Turbine  Efficiency  (6%  Decrease) 

= High  Pressure  Turbine  (Nozzle)  Area  (3%  Increase) 
- Low  Pressure  Turbine  (Nozzle)  Area  (3%  Increase) 

= Exhause  Nozzle  Area  Increase  (5%) 

= Exhaust  Nozzle  Area  Decrease  (4%) 

= Ninth  Stage  Bleed  Flow  (4  Ib/sec) 

= Twelfth  Stage  Bleed  Flow  (5  Ib/sec) 

= Sixteenth  Stage  Bleed  Flow  (3  Ib/sec) 


TABLE  D-I I 

Measured  Performance  Parameters 

Tx2  Inlet  Air  Total  Temperature 

PX2  Inlet  Air  Total  Pressure 

Tt3  Low  Pressure  Corap.  Disch.  Total  Temperature 

Px3  Low  Pressure  Comp.  Disch.  Total  Pressure 

Tt4  High  Pressure  Comp.  Disch.  Total  Temperature 

Pt4  High  Pressure  Comp.  Disch.  Total  Pressure 

P^  Burner  Pressure 

Tx7  Low  Pressure  Turbine  Disch.  Total  Temperature 

PX7  Low  Pressure  Turbine  Disch.  Total  Pressure 

N1  Low  Pressure  Rotor  Speed 

N2  High  Pressure  Rotor  Speed 

Wf  Fuel  Flow 

The  results  of  this  study  are  expressed  in  the  series  of  plots 
shown  in  Figures  D-3  through  D-14.  Each  figure  shows  the  effect  of 
error  in  a given  parameter  on  the  diagnostic  results  for  each  of  the 
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faults  analysed.  The  diagnostic  program  provides  a measure  of  the 
accuracy  of  the  fault  signature  match  called  Confidence  Factor,  In 
the  upper  plot  of  Figures  D-3  through  D-14,  the  sensitivity  study 
results  were  expressed  in  terms  of  percent  parameter  error  versus 
percent  decrease  in  the  Confidence  Factor  from  the  zero  error  factor 
(which  in  most  cases  was  over  9Q7„) . Also  shown  in  the  lower  plot  of 
each  figure  Is  the  effect  of  parameter  error  on  percent  change  in 
calculated  fault  level,  which  is  used  in  the  calculation  of  Confidence 
Factor. 

The  results  shown  in  Figures  D-3  through  D-14  are  summarized  in 
Table  D-III.  The  parameter  error  criteria  used  was  a 10%  allowable 
decrease  in  Confidence  Factor.  Using  this  criteria,  the  smallest 
parameter  error  (+  or  -)  that  would  cause  this  10%  decrease  was  deter- 
mined from  the  plots  for  each  parameter  and  each  fault.  These  allowable 
errors  are  listed  in  Table  U-III,  with  the  average  allowable  error 
computed  for  each  parameter,  and  the  parameters  ranked  in  order  of 
average  allowable  error  (sensitivity).  It  is  seen  that,  percentage- 
wise, errors  in  N2  and  N1  will  affect  the  accuracy  of  the  diagnostic 
results  the  most,  with  average  allowable  errors  of  0.48%  and  0.65%, 
respectively.  These  average  allowable  errors  only  serve  to  provide 
an  overall  measure  and  ranking  of  parameters.  The  sensitivity  of  each 
parameter  can  vary  considerably  with  each  fault,  and  the  required 
accuracy  for  each  parameter  must  be  bases  upon  the  most  sensitive 
case.  However,  this  ranking  by  average  allowable  error  is  very 
similar  to  a ranking  by  smallest  parameter  error  to  cause  a 10%  de- 
crease in  Confidence  Factor.  These  figures,  ranging  from  0.1%  for 
N2  to  1.7%  for  T2,  are  listed  in  the  last  column  of  Table  D-III. 

The  results  show  T2  to  be  the  least  sensitive  parameter,  but  this 
is  due  to  the  basis  used,  and  the  actual  T2  error  involved  is  not 
large.  It  should  be  noted  that  the  parameter  errors  introducted  were 
Ip  terms  of  percent  of  point,  in  the  measured  engineering  units,  i.e., 
in  the  case  of  T2,  the  point  value  is  59°F,  and  a 1%  error  represents 
a 0.59°F  error.  At  an  inlet  temperature  condition  of  12°F,  a 0.59° 
error  would  be  a 4,9%  of  point.  For  the  performance  analysis  compu- 
tations, T2  is  converted  from  °F  to  the  absolute  temperature  in  °R. 

Now  the  same  0.59°  error  becomes  only  0.11%  of  point  in  terms  of  the 
absolute  temperature  (519°R),  so  the  percent  error  is  considerably 
reduced.  The  same  is  true  of  T3  and  T4,  but  these  temperatures  are 
at  higher  levels,  so  the  conversion  to  absolute  temperature  has  a 
lesser,  but  still  significant  effect. 

In  order  to  bring  these  percent  error  numbers  into  perspective. 

Table  D-IV  is  presented.  This  lists  the  actual  value  of  each  para- 
meter at  the  sea  level  static,  standard  day  engine  operating  condition 
chosen  and  the  engineering  unit  equivalent  of  a 1%  of  point  error. 
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TABLE  D-IV 

Sea  Level  Static  Standard  Day  Values 


Point  Value 


Parameter 

at  T5=2200°R 

1%  Error  = 

N2 

13,701  rpm 

137  rpm 

N1 

8,951  rpm 

90  rpm 

Pt7 

60.52"  HgA 

0.61"  Hg 

Pt2 

29.92"  HgA 

0.30"  Hg 

?T4 

435.03"  HgA 

4.35"  Hg 

Wf 

5914  pph 

59  pph 

Tt4 

761°F 

o^ 

0 

Pb(Pt5) 

201.21  psia 

2.0  psi 

PT3 

153.76"  HgA 

1.54"  Hg 

Tt3 

407°F 

4.1° 

Tx7 

967°F 

9.7° 

Tt2 

59»F 

0.59° 

The  Table  D-IV  data  is  deck  baseline  data  with  no  faults,  so  the 
point  value  will  vary  somewhat,  depending  on  the  fault,  but  the  1% 
column  is  a good  approximation  of  the  magnitude  of  the  errors  involved. 

It  is  of  interest  to  analyze  the  TEDD  Phase  II  performance  para- 
meter data  for  repeatability.  This  information  is  readily  available 
from  the  TEDD  steady-state  Calibration  and  Data  Analysis  Program  which 
statistically  examined  the  data  to  identify  problem  areas.  Table  D-V 
shows  data  taken  at  a steady-state  operating  condition.  This  data 
reading  consists  of  200  consecutive  samples  on  each  channel,  for  which 
an  average,  standard  deviation,  and  statistical  average  are  computed. 
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TABLE  D-V 


Statistical  Analysis 

of  TEUD  Phase 

II  Data 

Paramete; 

Standard 

Deviatlon(s) 

Statistical 

Average 

N2 

!.3,578  rpm 

40.4** 

13,575 

N1 

8, 7-^5 

9.8** 

8,754 

FT7 

59.47"  HgA 

0.05 

59.47 

Pt2 

29.80"  HgA 

0.17** 

29.81 

?T4 

415.1"  HgA 

0.58 

415.1 

Wf 

1110  cps  * 

5.1** 

1111 

Tt4 

708. 9°F 

4.0** 

709.1 

Pb 

187.9  psia 

0.44 

188.0 

Pt3 

148.9"  HgA 

0.13 

148.9 

Tt3 

393. 8°F 

7.3** 

394.0 

Tt7 

1035. 2°F 

3.9 

1035.3 

Tt2 

60.6“F 

2.1** 

60.6 

* Flowmeter  output  signal 

**  Repeatability  unsatisfactory  without  data  smoothing 

The  average  is  a mathematical  average  of  200  consecutive  samples. 
The  standard  deviation  is  estimated  by  S as  follows: 


S = 


E (Xi  - X)2 
i=l 

n-1 


Where  Xj  is  the  ith  sample, 
X is  the  mean  of  ^ samples 


Statistically,  if  the  sample  is  random  (if  it  has  no  periodic  noise), 
95%  of  all  samples  will  fall  within  a band,  centered  on  the  mean, 
which  is  approximately  4S  wide.  Any  samples  which  fall  outside  this 
band  may  be  conc^dered  insignificant,  and  may  be  ignored.  The  statis- 
tical average  is  a new  mean  computed  from  the  points  remaining  after 
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all  insignificant  points  are  rejected.  Ignoring  the  insignificant 
samples  (5%),  the  data  scatter  is  approximately  + 2S.  Comparing  the 
actual  scatter  shown  in  Table  D-V  with  the  maximum  allowable  % error 
shown  in  the  last  column  of  Table  D-III,  it  is  seen  that  the  para- 
meters N2,  Nl,  Pt.?j  Wf>  Tt4,  Tt3  and  Tx2  had  too  much  scatter  to  use 
the  raw  data  without  smoothing.  The  limited  moving  average  smoothing 
technique  employed  in  TEDT  Phase  II  reduced  the  scatter  in  the  per- 
formance analysis  results.  However,  further  investigation  was 
required  in  this  area,  as  described  in  Appendix  A. 

The  characteristics  of  the  sensitivity  plots  in  Figures  D-3 
through  D-14,  are  unique  to  the  TEDD  Steady-State  Performance  Analysis 
Program.  Occasional  anomalies  such  as  a large  change  in  slope  or  a 
reversal  can  be  caused  by  changeover  from  one  range  to  another  or 
other  aspects  of  the  analysis  approach.  If  a plot  flattens  out  at 
higher  error  levels,  i.e.,  no  more  decrease  in  Confidence  Factor  with 
increasing  error,  this  is  due  to  the  "saturation"  (reaching  100%)  of 
the  various  individual  Confidence  Values  computed  which  use  the  subject 
parameter,  plus  the  non-use  of  this  parameter  in  the  average  fault 
level  computation.  In  actual  use,  the  instrumentation  check  portion 
of  the  program  would  have  identified  the  problem  of  bad  instrumentation 
before  most  of  these  error  values  were  reached.  To  fully  explain  the 
plot  characteristics,  the  details  of  the  analysis  program  must  be 
understood  (Reference  I>-2),  but  this  is  not  necessary  to  use  the  data. 

The  calculated  fault  level  (lower  plot  of  Figures  D-3  through 
D-14)  is  determined  by  the  program  just  prior  to  the  Confidence  Factor, 
and  the  results  of  the  fault  level  calculation  affect  the  Confidence 
Factor.  For  most  of  the  faults  only  3 or  4 curves  are  used  for  fault 
level  calculation.  These  curves  have  been  pre-selected  for  each  fault 
as  exhibiting  large  and  linear  changes  with  increasing  fault  level. 

Thus,  some  parameters  are  not  used  for  some  fault  level  calculations; 
and  therefore,  error  in  these  parameters  will  not  affect  these  cal- 
culations, and  there  is  no  corresponding  curve  shown  for  that  fault. 

Since  an  error  in  calculated  fault  will  cause  an  error  in  Confidence 
Factor,  it  is  generally  true  that  the  larger  the  change  in  calculated 
fault  level,  the  larger  the  decrease  in  Confidence  Factor  tends  to  be. 

An  area  where  the  program  was  not  fully  refined,  was  in  faults 
No.  4,  5,  and  6;  and  this  had  an  effect  on  the  results  in  these  cases. 

For  these  faults,  all  or  substantially  all  of  the  23  curves  were  used 
to  calculate  fault  level,  instead  of  selecting  a representative  few. 

In  addition,  all  23  curves  (rather  than  a selected  group)  were  used  to 
calculate  Confidence  Factor,  As  a result,  these  three  faults  show  up 
as  being  the  most  sensitive  to  parameter  error  In  the  majority  of  cases. 
An  example  is  seen  in  Figure  D-12  (TT3) . This  parameter  is  not  used 
to  determine  the  calculated  fault  level  for  any  fault  except  these  three, 
and  with  further  program  refinement  would  not  be  used  in  these  cases 
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either.  This  parameter  is  relatively  insensitive  to  changes  in  engine 
condition  due  to  these  faults,  so  conversely,  if  it  is  utilized  in 
the  diagnosis,  an  error  or  change  in  this  parameter  would  falsely 
Indicate  a relatively  large  charge  in  the  engine.  This  shows  the 
icipor  ance  of  pr'^ner  selection  of  curves.  Making  this  correction  in 
the  program  voula  increase  many  of  the  allovable  errors  to  some 
extent . 

E.  SUMMiKRY 

This  sensitivity  study  is  an  indicati  i of  the  kinds  of  accuracies 
required  for  the  most  critical  paran'.ete's  used  on  the  rEDD/TF30-P-408 
system.  In  addressing  accuracy,  both  the  absolute  accuracy  and  the 
repeatability  were  considered.  The  absolute  accuracy  is  important  in 
calculations  such  as  turbine  life  and  in  critical  limit  exceedances. 
However,  in  the  aelta  approach  used  for  steady  state  performance  analysis, 
it  is  not  the  absolute  accuracy  but  the  data  repeatability  that  is  most 
important,  and  this  can  be  significantly  improved  by  data  smoothing 
techniques.  The  accuracies  (repeatabilities)  required  are  subject  to 
variation  depending  on  the  type  of  engine,  the  faults  to  be  analyzed, 
and  the  type  of  diagnostic  system.  For  any  given  engine/diagnostic 
system  the  required  accuracy  also  can  vary  with  engine  flight  con- 
diton,  power  setting,  level  of  fault  detection  desired,  and  multipli- 
city of  faults.  Also,  these  results  reflect  only  the  error  in  single 
parameters  and  not  a combination  of  error  in  different  parameters. 

Two  parameters  which  were  not  included  In  this  TEDD/TF30-P-408 
study,  but  will  be  important  for  advanced  engines,  for  performance 
analysis  and/or  limit  exceedance  monitoring,  are  position  measurement 
(variable  compressor,  turbine  and  nozzle;  power  lever)  and  turbine 
inlet  temperature. 
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